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Let W be a rank n irreducible finite reflection group and let pi{x),. .. ,p n {x), x G R n , be a basis of 
algebraically independent W-invariant real homogeneous polynomials. The orbit map p : R n —> R n : 
x — > (pi(x),... ,Pn{x)) induces a diffeomorphism between the orbit space R 71 /W and the set S = 
p(R") C R". The border of S is the p image of the set of reflecting hyperplanes of W. With a given 
basic set of invariant polynomials it is possible to build an n x n polynomial matrix, P(p), p G R", 
sometimes called P-matrix, such that P a b(p(x)) = Vp a (x) • Vpt(x), Vo, b = 1,..., n. The border 
of 5 is contained in the algebraic surface det(P(p)) = 0, sometimes called discriminant, and the 
polynomial det(P(p)) satisfies a system of differential equations that depends on an n-dimensional 
polynomial vector A(p). Possible applications concern phase transitions and singularities. If the 
rank n is large, the matrix P(p) is in general difficult to calculate. In this article I suggest a choice 
of the basic invariant polynomials for all the reflection groups of type S n , A n , B n , D n , VnGN, for 
which I give generating formulas for the corresponding P-matrices and A-vectors. These P-matrices 
can be written, almost completely, as sums of block Hankel matrices. Transformation formulas allow 
to determine easily both the P-matrix and the A-vector in any other basis of invariant polynomials. 

Examples of transformations into flat bases, a-bases, and canonical bases, are considered. 
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I. INTRODUCTION 

When a rank n irreducible finite reflection group W acts in R", there exist a basis of n algebraically inde¬ 
pendent W-invariant real homogeneous polynomial functions pi(x ),... ,p n (x), x G R n , such that all W-invariant 
polynomial (or C°°) functions q(x) can be uniquely written as polynomial (or C°°) functions of the basic set: 
q(x) = q{pi{x),... ,p n {x)), with q(p) a polynomial (or C°°) function of p = (p±,... ,p n ) G R™ [H, |45|. By saying 
that q(x) is a W-invariant function , one means that q(gx) = q{x), Vg G W, Vi G R. n . The orbit (of the W action) 
containing x G M n is the set Wx = {gx | g G W}. Any two orbits are either coincident or disjoint. The orbit space of 
the W-action in is the quotient space W l /W, whose points are the orbits. 

The reflecting hyperplanes of W divide the space R n in open subsets, called chambers. The reflecting hyperplanes 
in the border of a given chamber C are called the walls of C and their subsets belonging to the topological closure C 
of C are called the faces of C. The closure of a chamber is a fundamental domain of the W-action in R", and W acts 
transitively on the chambers (13: Ch. V, §3.1, Lemma 2, and §3.3). This means that the closure C of a chamber C 
contains exactly one point for each orbit, that there are ord(TT) chambers, and that, given any two chambers C± and 
6 * 2 , there exists a unique element g £ W that maps in a one to one manner C\ in C'a. Given a chamber C, if x G C, 
the orbit Wx has exactly ord(IT) points, one for each chamber, and is called a regular or principal orbit , while, if x 
is a point in a face of C , the orbit Wx contains less than ord(W) points, all lying in reflecting hyperplanes, and is 
called a singular orbit. 

The W-invariant polynomial functions are constant on the orbits of R ra , so it is well defined the map 
p : R” — > R™ : x —> (pi(x),... ,p n (x)), called the orbit map. The basic invariant polynomials separate the orbits, 
that is, the points of R™ p = p{x ) and p' = p{x') are different if and only if the orbits Wx and Wx' are different: 
p ^ p' O Wx Wx'. The orbit map p then represents in a one to one manner the orbits as points of R n and the 


vittorino.talamini@uniud.it 



2 


orbit space R n jW as a subset S = pfW 1 ) of 1". Actually, S = p(C), where C is the closure of any chamber of W, 
and p is a diffeomorphism between C and S. This implies that S has the geometrical shape of a simplicial cone of 
R n , as that one of C, with its vertex at the origin, because of the homogeneity and the positive degrees of the basic 
invariant polynomials. The interior of S contains the images through the orbit map of the principal orbits, while the 
border of S contains the images through the orbit map of the singular orbits. From a different point of view, the set 
S is a simple example of orbifold (& Prop. 13.2.1). 

With a given set of basic invariant polynomials p a {x ), a = l,...,n, it is possible to build an n x n matrix 
function, P(p ), with matrix elements that are polynomial functions of p = (pi,... ,p n ) £ R n , that has the property 
that when the variables p a are substituted with the basic invariant polynomials p a (x ), Va = 1 its matrix 

elements P a b(p(x )) coincide with the invariant polynomials P a b{x), Vo,6 = 1 ,... ,n, obtained by making the scalar 
products among the gradients of the basic invariant polynomials, that is: P a b{p{x )) = P a b{x) = Vp a {x) ■ \7pb(x ), 
V a, b = 1,..., n. Following Ref. lj, I call the matrix P(p) with the name P-matrix, but the matrix P(p) has also been 
used with other notations and called “convolution of invariants” matrix [§], “complete displacement of invariants” 
matrix J24| or it was used without giving it a specific name in Refs. 0, Hil, [HT| . The P-matrix is strongly related 
with the set S: S coincides with the subset of R™ in which the P-matrix P(p) is positive semi-definite 0, H| , the 
border of S is contained in the algebraic surface determined by the equation det(P(p)) = 0, called discriminant , and 
at the interior of S one has det(P(p)) > 0. In addition, det(P(p)) satisfies a system of differential equations [40l [5lj. 
Eq. m below, called boundary equation , that depends on the matrix P(p) and on an n-dimensional polynomial 

vector A(p). The first appearance of the P-matrices has been in Ref._0, just for the groups A n , and their use 

for other finite reflection group has been first made in Refs. [HI.[36l [dill.[HUl52j| . The A-vectors first appeared in Ref. j5lj . 

Both the P-matrix P(p) and the A-vector A (p) depend on the choice of the basic invariant polynomials 
Pi(x),... ,p n (x). The explicit determination of the P-matrix from its definition requires to calculate the matrix 
elements P a b(x) = S7p a (x) ■ Vpb{x), and to express them as polynomials of the basic invariant polynomials. When 
the rank n of the group is large, these calculations require an enormous effort, even with a computer, and are often 
impossible to take to completion in a reasonable time. The possibility to use generating formulas to determine the 
matrix elements of the P-matrices are then particularly welcome. In Ref. Q, Remark 3.9, Arnold writes a generating 
formula for the P-matrices corresponding to the groups of type A n , saying that that formula was communicated to 
him by D. B. Fuks, and gives a hint for its proof (with a typographical error, a definition <7o = 0 that should be 
a o = 1). The formula is correct but the proof of that formula, to my knowledge, has never been published. In Ref. 
[24j], Givental writes generating formulas for the P-matrices corresponding to groups of type B n and D n . His results 
are correct but his proofs are only sketched and based on differential geometry and singularity theory techniques that 
are difficult to follow by a non-specialist. Without taking into account the generating formulas given by Arnold [1] 
and Givental j24j, I was looking for bases of invariant polynomials in which the P-matrices look particularly simple 
and with many similarities when comparing the P-matrices corresponding to groups of the same type and with rank 
differing by one unit. This allowed me do discover bases of invariant polynomials for which one can use generating 
formulas to calculate the matrix elements of the corresponding P-matrices. These generating formulas are Eq. HO), 
for the groups of type S n . Eq. (l32l) . for the groups of type A n , Eq. for the groups of type B n , and Eq. OT, 

for the groups of type D n . They are valid for all values of the ranks n £ N. Unexpectedly, the bases I chose for the 
groups of type B n and D n are the same chosen by Givental while those I chose for the groups of type A n differ only 
for a scalar multiples from those chosen by Arnold. I give proofs of these generating formulas that are based only on 
elementary calculus, algebraic manipulations and induction principle, so the proofs I give for the groups A n , B n and 
D n are very different from those suggested by Givental in |24| . although somewhat longer. These four generating 
formulas, with their algebraic proofs, are the main result of this article. In correspondence with the chosen sets of 
basic invariant polynomials, generating formulas for the A-vectors that appear in the boundary equation m are 
also determined, for all the groups of type S n , A n , B n , D n , and for all ranks n G N. These results are given in 
Eqs. (|26|) , (J33J) , (J30j) (|38j) and USD- The proofs are based on Theorems [T] and [2] below, that are new, and use sums 
over the positive roots of the reflection group W. Almost the same A-vectors, for the groups of type A n , B n and 
D n , were obtained in Ref. (52j by Yano and Sekiguchi using the induction principle on the rank n of the reflection 
group W, and the explicit expressions of the fundamental anti-invariant, product of the linear forms vanishing on the 
reflection hyperplanes of W. These different proofs for the same generating functions, both for the P-matrices and 
for the A-vectors, are another sign of the linking role that the finite reflection groups have among different branches 
of mathematics and physics. 
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To keep simpler the reading of the article, the generating formulas of the P-matrices and of the A-vectors for the 
groups of type S n , A n , B n , D n , are reported in Sections HTTlIVIl respectively, without proofs. The proofs, some of 
them quite long, are collected in a separate section, Section Ixl 

As a byproduct, the P-matrices obtained from the generating formulas for the groups of type A n , B n , D n , show a 
high level of symmetry and can be written, almost completely, as sums of block Hankel matrices. These results are 
reported in Section [VIII without proofs and their proofs are reported in Section III 

Several authors determined explicit bases of invariant polynomials for the irreducible finite reflection groups, 
for example Coxeter [18|, Ignatenko [28], Mehta 31], and many others that I do not cite. A few authors tried to 
select, among the infinitely many possibilities, some distinguished bases of invariant polynomials, by requiring some 
supplementary conditions to be satisfied. For example, Saito, Yano and Sekiguchi [39| defined the flat bases, Sartori 
and Talamini [H the a-bases and Flatto [2l| and Iwasaki [29j the canonical bases. Their results will be recalled at the 
end of next section. The bases of invariant polynomials that are here suggested are mainly justified by the possibility 
to use generating functions to determine the corresponding P-matrices and A-vectors. It turns out, however, that the 
matrix elements of the P-matrices obtained from the generating functions here reported (Eqs. (l25l) . (1321) . (l35l) . and 
(140|) ) are polynomials (with only integer coefficients) with a small number of terms, generally much less than those 
obtained using other bases of invariant polynomials. At the end of next section I report the number of terms of the 
matrix elements of the P-matrices of Bg in 4 different bases of invariant polynomials, precisely, the basis in Eq. (1341) . 
a flat basis, an a-basis and a canonical basis. 

The generating formulas for the P-matrices and the A-vectors here reported for the groups of type S n , A n , B n 
and D n , correspond to a particular choice of the basic invariant polynomials. However, one may easily obtain 
the P-matrices and the A-vectors corresponding to any other basis of invariant polynomials of these groups, 
using the transformation formulas between bases of invariant polynomials and the transformation formulas for the 
P-matrices and the A-vectors, given in Eqs. (ED and ED- Some examples, reported in Section HXl show how one can 
determine in this way the P-matrices and the A-vectors corresponding to a general flat basis, a-basis or canonical basis. 

The following Section [TT] contains more or less known results, but also some new ones. The proofs of the statements 
reported in Section El are generally omitted, because they can be found in the literature, except those of the new 
Theorems |T] and [5] and of a few other results, collected in Theorems El Eland [5J that will be widely used in the sequel. 
The proofs written in Section El are generally short and simple. 

Section fVIIIl deals very shortly with the exceptional and the non-crystallographic reflection groups. 

Orbit spaces and P-matrices have interes ting ap plications in mathematics and physics, especially in the study of 
phase transitions and symmetry breaking flL I26U321. l43j] . in the study of singularities H, H3, [36], and in constructive 
invariant theory j4Ej. In the rest of this introduction, I recall some of these applications of orbit spaces and P-matrices. 
This part is not necessary to follow the rest of the article. 

• Phase transitions and symmetry breakings. 

The isotropy subgroup, or stabilizer, of a point x € K™ is the subgroup W x of W such that W x = {g £ W | gx = 
x}. Points in a same orbit have conjugated isotropy subgroups. The orbit type of the orbit Wx is the conjugacy 
class [W x ] of subgroups of W. An orbit type stratum is a set of points x G K n that have conjugated isotropy 
subgroups. The number of orbit type strata is finite and coincides with the number of conjugacy classes of 
isotropy subgroups of W. The images of the orbit type strata through the orbit map form the orbit type strata 
of S. An orbit type stratum of S of dimension r, with r = 0,1,..., n, is the union of a finite positive number 
of connected components of dimension r, called primary strata. If r > 0, the boundary of an r-dimensional 
primary stratum is the union of primary strata of dimension lower than r. Each primary stratum of dimension 
r is the image through the orbit map of the set of fixed points of an isotropy subgroup of W of rank n — r. The 
decomposition of obit type strata in primary strata de pen ds on the reflection group W and is well known for all 
irreducible reflection groups (see for example Refs. (91. Il4l] for Weyl groups of simple Lie algebras, and (22|, Ch. 
3, for all irreducible finite reflection groups, including the non-crystallographic ones). 

If IT is a general compact linear group with a basis of n algebraically independent homogeneous polynomials, 
the union of the primary strata of dimension r, Vr = 0,1 ,..., n, coincides with the subset of K" in which the 
P-matrix P(p) has rank r and is positive semi-definite, if r < n, or positive definite, if r = n 0 , 111 . These 
conditions on the rank and on the (semi-)positivity of the matrix P{p), can be translated, using the minors of 
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P(p), in a (generally complicated and overdetermined) system of algebraic equations and inequalities on the 
variables Pi,. ■ ■ ,p n of R". (However, in the case of an irreducible finite reflection group W, the easiest way to 
find the algebraic equations and inequalities defining a primary stratum a of S of dimension r = 0,..., n — 1, 
is to consider the linear equations and inequalities defining the intersection of n — r faces of a chamber G, 
that has a for image in S through the orbit map p. Using these linear equations and inequalities as linear 
equations and inequalities among parameters in the orbit map, one obtains a set of n parametric equations 
and some inequalities, that depend on r free parameters, and define the primary stratum a). Moreover, in any 
point of the discriminant surface, the rows of the P-nratrix are a set of n vectors tangent to the discriminant 
surface. With some more complications, one can also determine, in a similar way, starting from the P-matrix, 
also the stratification of the orbit space of a general compact linear group acting in R ra , having a basis of 
to > n algebraically dependent polynomials. In this case one has to take into account the syzygies, that is 
the polynomials equations in the variables pi,...,p m , that are identically satisfied when these variables are 
substituted with the basic invariant polynomials pi(x),... ,p m (x), x £ R”. The set S C R m , representing the 
orbit space, is contained, in this case, in the intersection Z C R m of all the independent syzygy surfaces [34}. 


Let G be a rank n simple Lie group with Weyl group W. W is a rank n irreducible finite reflection group. 
The adjoint representation of G, indicated with Ad(G), is a representation of G in its Lie algebra g, that is a 
real vector space of dimension dim(g) = Card(P) + n, where 1Z is the set of roots of W. In Ref. [lE| Chevalley 
proved that it is always possible to obtain a basis of W-invariant polynomials from a basis of Ad(G)-invariant 
polynomials_by taking the restrictions of the basic Ad(G)-invariant polynomials to the Cartan subalgebra t) of 
g (see also [ll|, Ch. VIII, §8, no. 3). The bases of homogeneous W-invariant polynomials and of homogeneous 
Ad(G)-invariant polynomials contain then the same number n of elements with the same degrees, but they differ 
for the number of independent variables: the basic W-invariant polynomials depend on n variables, while the 
basic Ad(G)-invariant polynomials depend on dim(g) variables. One proves that the P-matrix corresponding to 
the W-invariant basis obtained by taking the restrictions to the Cartan subalgebra of the basic polynomials in an 
Ad(G)-invariant basis is identical to the P-matrix corresponding to the Ad(G)-invariant basis ([42}, Section III). 
This implies that W and Ad(G) have orbit spaces represented by the same set <S, stratified in orbit type strata 
and primary strata in exactly the same way. The isotropy subgroups corresponding to the different primary 
strata are, in the two cases, subgroups of W or of Ad(G), and the ones are the Weyl groups of the others. The 
P-matrices here obtained for the irreducible finite reflection groups of type A n , B n and D n can then be used 
also for the adjoint representations of the compact simple Lie groups of type A n . B n , C n and D n (the basic 
invariant polynomials of B n and C n are the same, and so also their orbit spaces). 

The W-invariant G°°-functions f(x) can be studied as functions f(p) defined in S , if f(x) = f(p(x)), Vx £ 
Dom(/) C R", however, f(p) is, in general, defined also outside S. The study of extrema of invariant functions 
can also be performed in the set S as a study of constrained extrema in the primary strata of the set S. If 
V(x) is a W-invariant potential function, and xq £ R" is a point in which V(x) takes its minimum value, then 
V(p) is minimum in the point p(x o), because V(x) = V(p(x)). If [Wo] is the orbit type of the orbit through 
Xo , then Wo is the residual symmetry group of the minimum point Xq , and the orbit type [Wo] is the residual 
symmetry type of the orbit type stratum hosting the minimum point of V(x). If V{x) depends on some external 
physical parameters (like pressure, temperature, external magnetic or electric field intensity, etc., in the case of 
thermodynamic potentials, or coupling constants, mixing angles, etc., in the case of Higgs potentials), when these 
external physical parameters change with continuity, the minimum point changes its position with continuity. 
If the new position x' 0 lies in the same primary stratum of xo , the residual symmetry group is at most changed 
by conjugation, the same that happens if one moves xo along its orbit, so minimum points lying in the same 
primary stratum represent the same residual symmetry type of the physical theory. If instead the new position 
x' 0 lies in another primary stratum, the two primary strata containing xq and x' 0 , respectively, have different 
dimensions and must necessarily be bordering. Then the residual symmetry type [Wo] is changed to [Wq], with 
the groups in [Wo] proper subgroups of groups in [Wq], (or vice versa, with Wo and W( exchanged). In this 
situation one has a phase transition (if V is a thermodynamic potential), or a (spontaneous) symmetry breaking 
(if V is a Higgs potential). (When the minimum point changes with continuity it cannot shift from a primary 
stratum of type [Wo] to another primary stratum of type [Wo] without passing through an orbit type stratum 
of different type [Wq] , this is the reason to speak only about primary strata and not about orbit type strata in 
the previous sentences). The set of orbit type strata then represent the set of different residual symmetry types 
possible for the physical system, and the bordering relations between primary strata represent possible phase 
transitions allowed for the physical system. It is easy to build graphs in which the nodes represent the possible 
primary strata (or orbit type strata) and the lines connecting the nodes the possible phase transitions. These 
graphs represent in a clear way all the possible phase transitions allowed to a given physical system. 
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• Singularity Theory. 

One of first goals of Singularity Theory was to study the degenerate critical points of smooth functions of one 
or several variables. It is better and simpler to work in complex spaces with complex holomorphic functions. 
A critical (or singular) point is a point zo £ C k in which the gradient Vf(z)\ z=zo of a function / : C k —> C 

vanishes and a degenerate critical point is a critical point zo in which also the Hessian det || <9 2 /(z)/ dzidzj |_ || 

vanishes. A critical point Zq £ C k is isolated if there is a sufficiently small neighbourhood U C C k such that Zq in 
the only critical point in U. One can suppose that zo = O, and that f(0) = 0, because one can always perform 
a translation of the origin of the system of coordinates of C k , and add a constant to /(z), without changing the 
behaviour of the given function at the given critical point. Near the origin O £ C fc , all germs of holomorphic 
functions, having O as a critical point and such that f(0) = 0, can be divided in equivalence classes, two germs 
being in the same class if there is a biholomorphic coordinate change in C fe , sending O in itself, that takes one 
germ into the other one. For each equivalence class one picks a polynomial representative of the germs in the 
class, that is called their normal form. This is always possible, if the critical point is isolated. Let /(z) be a 
normal form having a critical point in O in which f(0) = 0. If O is non-degenerate, /(z) is a quadratic form, 
while, if 0 is degenerate, the normal form /(z) has some terms with degrees greater than two. Two normal 
forms, that have O as a degenerate critical point, and differ only for the number of quadratic terms exhibit in 
O the same type of singularity. One can then study and classify the singularities by considering only normal 
forms without constant terms, with the smallest possible number of quadratic terms, and differing only for the 
terms of degrees higher than two. 

Families of functions depending on parameters necessarily exhibit degenerate critical points. One of such a 
families is a deformation (or universal unfolding) F(z, A) of the normal form /(z), with A £ C J a set of parameters 
controlling the deformation, such that F(z, 0) = f(z). The parameter space G ? is called the base space of the 
deformation. If O is an isolated degenerate critical point for the normal form /(z), with a generic deformation 
of f(z), the degeneracy in O disappears and bifurcates into a finite number /i(/) of non-degenerate critical 
points. The positive integer fi(f) is called the multiplicity , or Milnor number , of the singularity O of /(z). The 
Milnor number /i(/) can be calculated algebraically. Let Ok be the ring of holomorphic function germs in O, 
g : C k —> C, and If C Ok the ideal generated by the k partial derivatives of f(z). The quotient Af = Ok/Ff 
is called the local algebra of / in O. Its dimension, as a complex vector space is equal to the Milnor number 
of f(z) in O: n(f) = dime Af. One can take n = n(f) monomials <f>i(z), ..., 4> n (z ), as a vector space basis of 
Af (one of the basic monomials is always equal to 1), and write the deformation F(z, A), in the following way: 
F(z, A) = f(z) + A 0 ^> a (z). One obtains a germ deformation dependent on the least possible number of 

parameters, such that all other germ deformations can be reduced to it after a suitable biholomorphic coordinate 
change. This is a so-called miniversal deformation of /(z). For a miniversal deformation of /(z) one then has 
j = n = p,(f), and A £ C". One can suppose f(z) a weighted homogeneous polynomials, this implies that 
one has to assign a weight w% to each of the coordinates Zi, i = 1 ,,k, of z £ C fc , in such a way that all 
monomials in /(z) have the same weight w. One has to take for w and Wi the smallest positive integers for 
which /(z) is weighted homogeneous. A miniversal deformation F(z, A) of /(z) must then also be a weighted 
homogeneous polynomial of weight w, and this leads to assign weights d a also to the deformation parameters 
A 0 , Va — 1,..., n. 

Let f(z) = F(z, 0) have a degenerate critical point in O with critical value f(0) = 0 and Milnor number n. 
One is interested to determine all different singularities that can be realized with a deformation of the normal 
form f(z). The set of parameters A £ C™ in the base space for which F(z, A) has zero as a critical value, that 
is, for which there exists a point Za £ such that F(z\, A) = 0 and X7 z F(z, A)| z _ = O, forms a hypersurface 

£/ C C", called level bifurcation set of /(z) (or bifurcation set of zeroes of f(z)). This set can be found from 
the conditions for which the equation F(z, A) = 0 has multiple z roots, at least one, and this happens if and 
only if the discriminant Af of the polynomial £(z, A) is vanishing (one here considers z±,...,Zk as variables and 
Ai,..., X n as coefficients. The discriminants of polynomials in one or several variables are well described, for 
example, in Chapters 12-13 of [23|). The polynomial A/ depends only on the coefficients A a , a = 1,..., n, and 
is a weighted homogeneous polynomial, if F(z, A) is. The set of points A £ C n , satisfying the equation A/ = 0, 
forms an algebraic hypersurface of C™ that coincides with the bifurcation set of zeroes £/. For this reason E f 
is also called the discriminant set. We have then that if A £ £/, the equation in z, F(z,X) = 0 has multiple 
roots, and if A ^ £/, the equation in z. F(z, A) = 0 has all different roots. 

In [3|, Arnol’d discovers an unexpected link between the isolated singularities and the irreducible finite reflection 
groups of types A n . D n , E 6 , £ 7 , and E$ (these ADE groups are all the crystallographic finite reflection groups 
with only simple edges in their Coxeter-Dynkin diagrams). To understand this link, one has first to consider 
the action of a finite irreducible reflection group W on C". This action can be obtained with the same real 
orthogonal matrices of the action of W in R™, and this implies that the real and the imaginary part of a point 
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z £ C" transform without mixing, because, V g £ W and Mz £ C n , one has gz = g{x + iy) = gx + igy , with 
x, y £ R” and i = y/—l. The complex reflection hyperplanes have the same equations of the real reflection 
hyperplanes, but the variables involved are now complex instead of real. The basic invariant polynomials of 
this actions are the same as those of the action of W in R™, so, there is a unique corresponding P-matrix P(p) 
for both the real and the complex actions of W. In particular, the coefficients of the matrix elements of P(p) 
are real numbers, and det(P(p)) is a polynomial with real coefficients. In contrast, the independent variables 
on which P(p) depends are now complex numbers, because p £ C n . Some of the orbits in C n are formed by 
only real points, and these have image in S C C n , where S is now a subset of the set of points of C™ with 
vanishing imaginary part. Orbits formed by only complex numbers have images, through the orbit map, both 
in points with vanishing imaginary part (lying outside S in the real subspace of C n ), and in points with not 
vanishing imaginary part. In any case, the orbit space of the action of W in C" is the whole C": p( C n ) = C". 
A point 2 £ C" belonging to a complex reflecting hyperplane have image, through the orbit map, in a point 
p = p{z) £ C" in which det(P(p)) = 0. This because the polynomial det(P(p)) can be written as a product of 
squares of linear forms whose vanishing determine the reflection hyperplanes (Eq. ([!J]) below, in which x £ 1" 
has to be substituted by z £ C"). Then, when 2 belongs to a complex reflecting hyperplane det(P(p(z))) = 0. 
The set of the complex reflecting hyperplanes is then mapped to an algebraic surface E of C" with equation 
det(P(p)) = 0 (with p £ C n ). This surface is no more a boundary surface of C n , but just a set of points of C" 
that contains the images through the orbit map of the singular orbits passing through the reflecting hyperplanes. 

For all the simplest types of normal forms f(z), for which O is an isolated singularity with finite Milnor number 
n, and f(0) = 0, Arnold [2J found that there is one and only one rank n irreducible finite reflection group 
W, of one of the types A n , D n , E e , E 7 , and E s , such that there exist a biholomorphic map between the 
orbit space C n of the action of W in C" and the base space C ra of the miniversal deformation F(z, A) of f(z), 
keeping the origin O £ C n fixed, and such that the surface ScC" (containing the image through the orbit 
map of the reflection hyperplanes of the IF-action in C n ) is mapped to the surface E/ C C” (forming the 
discriminant set of the singularity). For this reason the image of the reflection hyperplanes in the complex 
orbit space is called discriminant. One often can choose, without much effort, the basic invariant polynomials 
of W in such a way to obtain E = E/. The equation of E/ can be obtained in this case using the P-matrix: 
E/ = {p £ C n | det(P(p)) = 0}. The weights d a that were assigned quite arbitrarily to the parameters A a of 
the deformation are then equal to the degrees of the basic invariant polynomials of the group W associated to 
the singularity: d a = w( A 0 ) = deg(p a ( 2 )), Va = 1,... ,n. As a singularity determines a finite reflection group, 
a finite reflection group determines a singularity. The simplest isolated singularities can then be classified using 
the same symbols of the classification of the irreducible finite reflection groups. There exist then singularities of 
types A n . D n , Eq , E 7 , and Eq. The existence of a biholomorphic map between the orbit space of the IF-action 
in C" and the base space of a miniversal deformation of f(z), such that the origin is mapped to the origin, 
and E C C" (the p-image of the complex reflection hyperplanes of W) is mapped to the discriminant set of the 
singularity E f C C n , was actually proved in general by Slodowy [43], using previous results by Brieskorn [11 ]. 

Other singularities can be associated in a similar manner to the crystallographic finite reflection groups with 
a non-simple edge in their Coxeter-Dynkin diagrams |4}. The groups of type B n , C n and F 4 correspond to 
singularities in O of normal forms invariant for the reflection about a hyperplane containing O, that is exhibiting 
a Z 2 symmetry (this reflecting hyperplane is called boundary and these singularities are thus called boundary 
singularities ), while the group Gi correspond to a singularity exhibiting a Z 3 symmetry. In all these cases there 
exists a biholomorphic map between the orbit space of the IF-action in C" and the base space of a miniversal 
deformation of the normal form of the singularity f(z), such that E C C n (the p-image of the complex reflection 
hyperplanes of IF) is mapped to the discriminant set E/ C C n , and the boundary and the origin are mapped to 
themselves. The non-crystallographic groups H 2 = T 2 (5), H 3 and H 4 are related to singularities that are formed 
by the time evolution of wave fronts avoiding obstacles (see Ref. 0]). In this case there exist a biholomorphic 
map between the p-image E of the complex hyperplanes of the non-crystallographic reflection group and a subset 
(the closure of a primary stratum) of the discriminant set E/ of the singularity (see Theorem 4 of Ref. 0). 

The relations between singularities and finite reflection groups are much deeper than how is here written. 
For example, one can prove that the finite reflection group that is associated to an isolated singularity is the 
monodromy group of the singularity that has its simple roots forming a basis of the vanishing cycles. There are 
many books and articles treating in many depth levels singularity theory, so I stop here my short review (for 
details see for example @, 0, SH tH3 ) ■ 

In Refs. EH Arnol’d used the P-matrices to obtain some other results in singularity theory. In particular, in 
1, he used a linearized form of the P-matrices, called linearized convolution or linearized displacement matrix, 
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in the cotangent space of the orbit space C n /W at the origin, to determine certain integral invariants, called 
indices of the singularities. This kind of calculations have been also done by Givental [3|, who first determined 
explicit formulas for the P-matrices of the groups £ n , D n , £4 and £ 6 . Saito [36] proved that on the base space 
of the universal unfolding of an isolated hypersurface singularity it is possible to define a flat structure, obtained 
through the period map associated to a primitive form. This flat structure implies the existence on the base 
space of the universal unfolding of a system of coordinates with a flat metric, i.e. with an everywhere vanishing 
curvature. The coordinates of the base space are consequently said flat coordinates. The original references 
are Ref. 27 of Ref. [36| (unpublished) and Ref. ]35| (all reviewed in Ref. [33)). Saito, Yano and Sekiguchi 
proved [36, 39|, just using standard properties of finite reflection groups, that in the orbit space C n /W ~ C™ 
of the effective action of an irreducible finite reflection group W onC", it is possible to choose a system of 
coordinates (that is, a basis of W-invariant polynomials that determines a system of coordinates in the orbit 
space, obtained by identifying the basic polynomials with the coordinates) with a real and constant flat metric. 
It is said consequently that the orbit space acquires a flat structure and that the corresponding basic invariant 
polynomials form a flat basis. To obtain a flat basis of invariant polynomials in Ref. j.3f| , Saito, Yano and 


Sekiguchi used the P-matrices. Even if they did not know the complete form of the P-matrices, they succeeded 
to determine the flat bases of all irreducible finite reflection groups, except £7 and Eg, just studying the terms 
of the P-matrices containing the highest degree basic invariant polynomial. The flat bases of the groups £7 
and £g, together with their P-matrices, were found in Ref. [49| using the method suggested in Ref. [§9j. In 
1991 Dubrovin, trying to give a coordinate independent formulation of the Witten Dijkstra-Verlinde-Verlinde 
associativity equations that appeared those years in the setting up of the two dimensional topological field 
theory, was lead to define a geometric object that he called Frobenius manifold. He also proved that the 
(complexified) orbit space of an irreducible finite reflection group, endowed with the flat structure introduced 
by Saito, satisfies all the axioms of a Frobenius manifold [Hj, |20| (in effect this is the most elementary example 
of Frobenius manifold). Recent articles describing the set up of the flat structure and the Frobenius manifold 
structure on the (complexified) orbit space of a finite reflection group and on the base space of the universal 
unfolding of an isolated hypersurface singularity are Refs. [37] and 


• Constructive invariant theory. 

The P-matrices can be defined for a general real compact matrix group G, not necessarily a finite group. All 
such groups admit a finite basis of real G-invariant homogeneous polynomials. The basic invariant polynomials 
are in general not algebraic independent and sometimes it may be difficult to determine all the polynomials in 
a basis. The construction of the P-matrix can help to overcome this difficulty. The degrees of some of the basic 
invariant polynomials, as well as the degrees of some of the syzygies among the basic invariant polynomials can 
be inferred from the Molien function of G, and this is the first calculation to do. Then, one has to determine 
all the basic invariant polynomials of the lowest degrees, for example by calculating the averages in the group 
G of some proper monomials. Let these invariant polynomials be pi(x),... ,Pk{x ), with the labeling increasing 
according to their degrees. These are supposed to be only some, and not all, of the basic invariant polynomials 
of G. Then one starts to determine the P-matrix elements with these basic polynomials, from the lowest degree 
element to the highest degree element, using the definition, that is by expressing the invariant scalar products 
X7p a {x) • \7pb(x) as polynomials of the known basic invariant polynomials pi(x),... ,pk(x). If one is not able 
do do this, then one sets pk+i(x) = ’S7p a {x) ■ Vpbix), and one discovers so a new basic invariant polynomial 
(and if deg (pk(x)) > deg(pfc+i(x)) one has to relabel the basic invariant polynomials). When one finishes the 
construction of the P-matrix and the degrees suggested by the Molien function are in accord with those of 
the basic polynomials so found, one has determined the full basis of invariant polynomials of G. There are 
shortcuts and particular cases. The interested reader can find a wider exposition and some examples in Ref. [48[ . 


II. BACKGROUND MATERIALS AND PRELIMINARY RESULTS 

In this review section I shall report mostly known results, but also some new ones, like Theorems [l] and [2] The 
arguments here reviewed concern the theory of finite reflection groups, the theory of linear actions of compact groups 
and of the use of invariant theory to describe their orbit spaces. The proofs are generally omitted, except for some 
theorems that are widely used in the sequel. For a deeper understanding of arguments related to linear actions of a 
general compact group or to the use of invariant theory in this context one may read Refs. 0 , m, h, while among 
the many good books on the theory of finite reflection groups and their relations with Lie groups I mainly used 



Refs. EaUa- There one can find the proofs or the references to many of the results that are here just stated. 

I consider here group actions in R n , in which it is given a set of n orthonormal canonical basis vectors ei,..., e„. A 
reflection R of R" is a linear transformation of R™ that leaves pointwise fixed a (n — l)-dimensional hyperplane -kr of 
R", called the reflecting hyperplane of R, and sends each x £ R" to its symmetric x' £ R" with respect to ttr. This 
means that £ ttr, and that x — x' 1. ttr. The line passing through the origin of R" orthogonal to ttr is mapped 
into itself by the reflection R but is reflected about the origin. It is called the reflecting line of R. 

Given any non-null vector a £ R n , the reflection R a that turns a into its opposite, can be defined by the following 
formula: 

R a x = x — ^ 01 X a Vx £ R" , (1) 

a ■ a 

where • indicates the canonical scalar product in R". Any non null vector proportional to a would define the same 
reflection. The reflection hyperplane of R a has equation a ■ x = 0, unique if one excludes possible multiplicative 
factors. We will indicate with 


l a (x) = a ■ x 


the linear form that vanishes in the reflection hyperplane of a. We will identify the reflection R a , with the orthogonal 
matrix that realizes the reflection R a in R". 

By multiplying a set of reflections in all possible ways one obtains a group of orthogonal matrices, called a reflection 
group. In general this group is infinite. The irreducible finite reflection groups require particular angles between the 
reflection lines and were classified by Coxeter na in the following types: A n . n > 1, B n , n > 2, D n , n > 4, Eq, 
E-j, Eg, F 4 , H 3 , H 4 , him), m > 5. The number appearing as a subscript is the rank n of the group and is equal to 
the dimension of the space relative to which the group action has no fixed points besides the origin. The limitations 
exclude group isomorphisms (and for this reason the groups of type C n are not in the list because, as reflection 
groups, they coincide with the groups B n , Vn > 3, and the group G 2 is not in the list because of the isomorphism 
12 ( 6 ) ~ G 2 ). A reflection has determinant equal to —1, and hence it is an improper rotation. The product of two 
reflections has determinant equal to +1, and hence it is a proper rotation. A finite reflection group has an equal 
number of proper and of improper rotations. Only a few number of the improper rotations are reflections, the others 
do not fix pointwise an hyperplane. 

Let n be a positive number and W C 0(n) a rank n irreducible finite reflection group. If we choose a vector r £ R n 
parallel to a reflecting line, the group W maps this vector into a finite set <J> r = {gr, Vg £ W} of vectors of R", all 
parallel to reflecting lines, and all of the same length, because W is orthogonal. For all vector v £ $ r , also — v £ <fv, 
because R v £ W. The set <fv is stable with respect to W transformations. The vectors in $ r are called roots and 
forms a root system (according to the definition of root system given in 27], Section 1.2). With certain reflection 
groups W ( B n , Vn > 2, F 4 , h(m), Vto > 6 , m even), there are reflecting lines that are not parallel to any root in 
<fv. Then it is necessary to define a second vector r' parallel to one of these reflecting lines left over and repeat the 
construction of a new root system <!>+. The root systems <J> r and <lv are disjoint, and at most two root systems are 
necessary to cover in this way all reflecting lines of W. The lengths of the roots r and r' at this point are arbitrary. 
One can determine, not in a unique way, a set of n linearly independent roots, called simple roots, such that all roots 
are linear combinations of the simple roots with real coefficients, all of the same sign. Those with positive coefficients 
are called positive roots. Let’s use the symbols 77 and 77+ for the set of roots and the set of positive roots of W. 
Clearly, 77 = {±r, Vr £ 77+}. There is a one to one correspondence between the set of reflections of W and the set 
77+, of positive roots of W. The reflections obtained with Eq. © from the simple roots are called simple reflections 
and are sufficient to generate the whole group W. The Z-span of the simple roots is called the root lattice. A finite 
reflection group W for which there exist a root lattice that is stable with respect to W transformations is called 
crystallographic. A finite reflection group with two different root systems that is crystallographic has necessarily the 
roots in the two root systems with lengths in the ratio \/2 or v/3- When the ratio is \/2, and n > 3, there exist 
two different W-stable root lattices, that correspond to two different crystallographic reflection groups: B n and C n . 
There is a one to one correspondence between the crystallographic irreducible finite reflection groups (A n , n > 1, 
B n , n > 2, C n , n > 3, D n , n > 4, Eq, E 7 , E$, F 4 , G 2 — h( 6 )) and the Weyl groups of the simple Lie algebras. 
The non crystallographic finite reflection groups are H 3 , H 4 and / 2 (m), m > 5, m ^ 6. The relative length of the 
roots of the two root systems of the groups him), m > 5, m / 6 is not fixed. If W is crystallographic, all roots 
belong to the root lattice, and this implies that they are linear combinations of the simple roots with integer coefficients. 


Consider now the algebra of polynomials of the n variables x \,..., x n . The W-invariant polynomials are such that 
p(gx) = p{x), V g £ W, Vx £ R ra . There exist n W- invariant real homogeneous polynomials p\{x),... ,p n (x), such 
that for every invariant polynomial p{x), there exists a unique polynomial p{p) in n indeterminates p = (pi,... ,p n ) £ 
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R”, such that p(x) = p(p\(x),... ,p n (x)), \/x £ R n . The polynomials pi(x),... ,p n (x) are called basic invariant 
polynomials. There are infinitely many possible choices of a set of basic invariant polynomials. For the algebraic 
independence of the basic invariant polynomials, an equation like p(p\(x),...,p n (x)) = 0 cannot be identically 
satisfied VieR". 

The degrees d \,..., d n , of the basic invariant polynomials of all the irreducible finite reflection groups, were deter¬ 
mined by Coxeter [IS], and for a given irreducible finite reflection group they are all different, except in the case of the 
groups D ni with even n, in which there are two basic invariant polynomials of degree n. Usually, the basic invariant 
polynomials are labeled according to their degrees, for example by taking 

da ^ d a - j_i , V fl — 1,...,7T 1. (2) 

If the group has no fixed points besides the origin, di = 2, and one may define 

n 

p ± (x) = X ■ x = \\x\\ 2 =^2xi 2 , ( 3 ) 

2—1 


a natural choice for real orthogonal actions. 

Coxeter proved [HI, that if the finite reflection group W is irreducible, the sum and the product of the degrees 
d a , a = 1,... ,n, are related to the number N of reflections of W and to the order, ord(IU), of W, by the following 
formulas: 


n 

^ d a = N + n , N = card(77. + ), (4) 

a =1 


n d a = ordW • (5) 

a =1 

Moreover, for W irreducible, given the usual degree labeling @ of the basic invariant polynomials, one has the 
identities: d a + d n - a +i = d n + 2, Vo = 1,..., > an d, using Eq. fU, these imply N = nd n /2. ( d n is the greatest 

degree of the basic invariant polynomials). 


With a basis {p±(x),... ,p n (x)} of invariant polynomials, one may calculate an n x n real symmetric matrix P(x) 
which has its elements P a b(x ) that are the scalar products of the gradients of the basic invariants SH: 


n 

Pab(x ) = Vpa(x) ■ S/pbix) = ^ 

i= 1 


dpajx) dp b {x) 
dxi dxi 


V o, b = 1 ,..., n. 


( 6 ) 


The matrix elements P a b{x) are homogeneous polynomials of degrees d a + db — 2. By convention I shall use the Latin 
indices a, &, c,..., to label the basic invariant polynomials, and the Latin indices i,j , fc,..., to label the coordinates of 
R” in their domain. From the homogeneity of the basic invariant polynomials and the standard form f3j of pi(x), 
one obtains, using Euler’s theorem on homogeneous polynomials, that the first row and column of P(x) have fixed 
form: P\ a {x) = P a i(x) = 2d a p a (x), Va = 1,... ,n, where no sum over a is understood. 

The matrix P(x) can also be written in the following way: 

P(x)=j(x)j T {x), (7) 


in which the exponent T means transposition and j(x ) is the jacobian matrix 


jai(x) 


dpg(x) 

dxi 


Va = 1,..., n, 


V? = 1 ,..., n . 


Then, P(x) is a positive semi-definite matrix, in fact x T P(x)x = > 0, Vx £ R". Consequently, all 

principal minors of P(x ) are non-negative, in any point x £ R". 

A classic result by Coxeter [l 8 j |. Theorem 6.2, claims that det(j(z)) is proportional to the (anti-invariant) product 
of the linear forms l r (x) whose vanishing determine the set of the reflecting hyperplanes, that is: 


det(j(a:)) = c l r (x), 

r€iz+ 


( 8 ) 
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in which 72+ is the set of positive roots of W, l r (x) = r ■ x, and c is a non-zero constant. Then, one has: 

det(P(x)) = c 2 %(x), (9) 

r£lZ + 


that implies det(P(x)) = 0 if x belongs to the set of the reflecting hyperplanes and det(P(x)) > 0 otherwise. The 
definition of l r (x ) as it is given above depends on the normalization of r, and the constant c then depends on the 
normalization of both the roots and the basic invariant polynomials. 

The matrix elements P a b{x ) are W- invariant polynomials. This follows from the orthogonality of W and the 
covariance of the gradients of IT-invariant functions. In fact, if p(x) is an invariant polynomial, p(x ) = p{x 9 ), 
V g £ W, V x £ R n , where x 9 = gx is the g- transformed point of x, and one has, V g £ W and V x £ R n : 


dp( x) 
dxi 

that implies: Vp(x) = 


dp(x 9 ) ^ dp{x 9 ) dx 9 dp(x 9 ) T dp(x 9 ) 

dxi ^ dx 9 dxi dx 9 ^ ^ dx 9 ’ 

3=1 3 3=1 3 3=1 3 

g T S7 9 p(x 9 ) or S7 9 p{x 9 ) = gS7p{x ), where V 9 = (^,..., gfr). Then: 
V 9 p Q (x 9 ) • \7 9 p b {x 9 ) = gX/p a {x) ■ g\7pb(x) = Vp a (x) • Vp 6 (x), 


that is, Pab(gx) = P a b(x), that tells us that the matrix elements of P(x), calculated with Eq. (0, using the variables 
x, or the transformed variables x 9 = gx, V g € W, are the same. Hence, P 0 b(x) is a IT-invariant polynomial. 

As a consequence, we have that all the matrix elements P a b(x) can be expressed as polynomials of the basic invariant 
polynomials: 

Pab(x) = Pab(pi{x),...,p n {x)), Vx € K™ . (10) 

One can then define a matrix function of p £ M™, P(p), having P a b{p) = Pab{pi, ■ ■ ■ ,Pn ) for elements, such that 
Eq. m is satisfied. The real symmetric matrix P(p) is called the P-matrix. Its explicit form depends strongly on 
the given basis Pi(x),... ,p n (x) of IT-invariant polynomials. 


It is useful to define an operator P that applied to a given basis of homogeneous invariant polynomials p a {x), 
a = 1,... ,n, gives the corresponding P-matrix P(p), obtained through the application of Eqs. 0 and ([TO)) . This 
will be understood in the following. 


The map p : R ra —> R™ : x —> p{x) = (p i(x),... ,p n (x)) is called the orbit map and maps M" into an n-dimensional 
semi-algebraic connected proper subset S of R". There is a one to one correspondence between the orbits in R" 
and the points in S , because the basic polynomials separate the orbits, that is, given any two orbits in R", at least 
one of the basic invariant polynomials p a {x), a = 1,.. . ,n, takes a different value in the two orbits. The orbit map 
p actually induces a diffeomorphism between the orbit space W 1 /W of the W action in R™ and the set S, so S can 
be identified with the orbit space W 1 /W. Points lying in a same orbit of R" have conjugate stabilizers, or isotropy 
subgroups. The set of points in R ra with conjugated isotropy subgroups is called a stratum of the W action in R". 
Its image in the set S through the orbit map is called a stratum of S. Clearly there is a one to one correspondence 
among strata of R n and strata of S. The number of strata is finite and equals the number of different conjugacy 
classes of isotropy subgroups of W. 

The P-matrix P(p) contains all information needed to characterize geometrically the connected set S C R™ 
representing the orbit space M n /W. For example, the interior of the set S is the only subset of R" where the 
P-matrix is positive definite EM. The equation detP(p) = 0 determines the discriminant, an (n — 1 [-dimensional 
algebraic surface in R™ that contains the image, through the orbit map, of the set of the reflection hyperplanes of R", 
and this image coincides with the boundary of the set S. The fc-dimensional strata of S are found from the equations 
and inequalities expressing the conditions: rank(P(p)) = k and P(p) > 0 (this last condition means semi-positivity). 
S is a semi-algebraic set because all these defining conditions can be obtained through polynomial equations and 
inequalities. This method to determine the algebraic equations and inequalities of the strata of the orbit space is 
however not very easy to perforin, and for the finite irreducible reflection groups there are better methods, shortly 
recalled in the introduction. 


As the basic invariant polynomials pilx),... ,p n (x) are homogeneous, it is natural to consider the coordinates 
pi,...,p n as a set of n graded, or weighted, coordinates whose degrees, or weights, coincide with the degrees of 
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the basic invariant polynomials, that is: w(p a ) = deg(p a (x)) = d a , Va = 1 The variables p and x are used 

to label points in different spaces 1" the variables p are weighted, while the variables x are not. As homogeneous 
polynomials are important in what follows, it is convenient to recall two definitions. A polynomial p(p) is said of 
weight d if the IT-invariant polynomial p(x) = p(p\(x),... ,p n (x)) is of degree d in x. If p(x) is homogeneous then 
p(p) is said weighted homogeneous , or w-homogeneous. Each coordinate p a € R” is itself a ui-homogeneous polynomial 
with weight d a = w(p a ), the matrix elements P a b(p) of the P-matrix P{p) are ic-homogeneous polynomials of 
weights d a + db — 2, and the determinant of the P-matrix, det(P(p)), is a w-homogeneous polynomial, of weight 
J2a= i(2da — 2) = 2N = nd n , if W is irreducible, where Eq. ((H) has been used, N is the number of roots of W, and 
d n is the highest degree among the basic invariant polynomials. 

A w- homogeneous irreducible factor a(p) of the determinant det(P(p)) of the P-matrix P(p) satisfies the following 
system of differential equations 0 , [5l|, called the boundary equation: 


C—1 


da(p) 

dp c 


[ a) (p)a(p ), 


\/b = 1 ,... ,n , 


( 11 ) 


where X^ a \p) = ^A^(p),..., A l?\p)^J is a polynomial vector function of p € R", dependent on a(p ), whose elements 
A [ a \p) are ui-homogeneous polynomial in the weighted variables pi,... ,p n , of weights db — 2, Mb = 1,..., n: 

w(X^\p)) = d b - 2 , V6=l,...,n. 

The relation w(X[ a \p)) = db — 2 easily follows from Eq. (fill) and the w- homogeneity of a(p) and of Pbc{p), and 
in particular from w(Pb c (p)) = db + d c — 2. The first element A ^(p) is a constant, because d\ = 2. From Eq. 
(HU), with 6=1, using Pi c (p) = 2 d c p c (no sum over c) and Euler’s theorem on homogeneous polynomials, that for 
re-homogeneous polynomials, like a(p), has the following form: 


n 

y, d c p c 

C=1 


da(p) 

dp c 


w(a) a(p ), 


one obtains the following general result, that fixes the form of A^(p): 

A i Q) (p) =2 w(a(p)). 

Eq. (fill) was first given in [5lJ , and further studied in (52|, for the case of an irreducible finite reflection group, and 
in which for a(p(x)) it was considered just the invariant polynomial det(P(x)). Eq. (fill) was rediscovered in [40|, and 
extended to the case of a general effective linear group having a basis of n algebraically independent polynomials, 
and in which for a(p) it was considered a general irreducible factor of det(P(p)). In (4l| there were further examined 
the general properties of arbitrary w -homogeneous polynomials a(p), satisfying Eq. (fill) together with a proper n- 
dimensional vector X^ a \p), whose components X < j ) a \p) are w-homogeneous polynomials of degree db — 2, b = 1,..., n. 
Two definitions are convenient. A re-homogeneous polynomial a(p), that satisfies the boundary equation (lllll . together 
with a proper vector A ( a \p), is called an active polynomial (of P(p), or corresponding to P(p)), and the vector X^ a \p), 
without imagination, is called the A -vector (of a(p), or relative to a(p)). Active polynomials differing only for constant 
multiplicative factors should be identified because they satisfy the boundary equation with the same A-vector. In [4l| . 
among other things, it was proved that: 

1 . all irreducible (in the complexes, as well as in the reals) active polynomials are irreducible factors of det(P(p)) ; 

2 . all products a(p) = ni=i( a j(p)) mi > €= N, of active irreducible polynomials ai(p), i = 1 ,...,r, are active 

polynomials too, for which A ( a \p) = X)l=i A^ 0i )(p), in which A( 0i ) (p), Vi = 1 ,..., r, are the A-vectors of the 
active irreducible polynomials a,i(p), i = 1 ,... ,r. 

The determinants of the P-matrices of the irreducible finite reflection groups are generally irreducible polynomials, 
except for the groups of type B n , Vn > 2, P 4 , Vm > 6 , even to, in which cases there are two irreducible active 

factors of det(P(p)). These are exactly the cases in which there are two different root systems. 
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Given a P -matrix P(p) of an irreducible finite reflection group, the boundary equation, Eq. (1111) . is satisfied by all 
the irreducible factors of det(P(p)), and in any case by det(P(p)) itself: 


Y Pab(p) 


6=1 


d det(P(p)) 
dpb 


= Ai det{P)) (ri det(P(p)), 


a = 1 ,..., n. 


( 12 ) 


There is a general formula giving the A-vector \( det ( p ^(p) corresponding to the active polynomial det(P(p)) that 
appears in the boundary equation m of a general irreducible reflection group. It is given by the following theorem. 


Theorem 1 Let W be an irreducible finite reflection group, and 1Z+ a set of positive roots of W . For all r G 
1Z + let l r (x) = r • x, be the linear form vanishing on the reflection hyperplane of r. The \-vector A( det ( p ))(p) = 
(A^ det ( p ))(p),..., An det * P ' ) ' ) (p)), that appears in the boundary equation 112 d . can be calculated with the following formula: 


Ai det(P »(p(*))=2 ^ 

r6R + 


S7p a (x) ■ r 
l r (x) 


Va = 1,..., n , 


(13) 


valid V x G R n . 

Proof. Eqs. (IToT) and 0 imply det(P(p(x))) = det(P(x)) = (det(j(x))) 2 . Using this result in Eq. (fl2l) . written in 
terms of the original variables x G R n , one has, at the first member: 


Y Pabip) 


d det (P(p)) 


6=1 


dpb 


= Y ^Pa( x ) ’ Vp b (a 

p=p(x) b—1 


d(det(j(x))' 2 

dpb{x) 


= Vp a (x) • V(det(j(^))) 2 = 


= 2 det(j(x)) Vp a (x) ■ V det (j(x)). 


Using Eq. (0), one obtains: 


/ 


V det(j(x)) = c V l r (x) = c Y J 


re n. 


ren+ 


\ 


S7l r (x) l r flx) = det(j(x)) Y VU(a;) = det{j{x)) Y 


V 


r' * i 


/ 


re n. 


l r (x) 


rG 7£+ 


l r (x) 


because V/ r ( x) = r. Using this result in the preceding equation, one obtains: 


^Pabip) 


<9det(P(p)) 


6=1 


dpb 


= 2 (det(j(x))) 2 Vp a (x) • Y JJffj = 2 det(P(p(x))) Y ■ 


p=p{x) 




ren+ 


l r (x) 


Comparing with the second member of Eq. G2D, written in terms of the original variables x G 


Ai det(P)) (p) det(P(p))(p) 


p=p(x) 


, one finds that 


Ai det(p)) (p(x))=2 Yi 


ren+ 


Vp a {x) • r 
l r (x) ’ 


a = 1 ,..., n. 


As the first member of Eq. in is an invariant polynomial in x G R”, also the second member is an invariant 
polynomial in x G R ra , and can be expressed as a polynomial of the basic invariant polynomials. 

A straightforward consequence is the following. 
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Theorem 2 Let W be an irreducible finite reflection group with two different root systems, that is W is one of the 
groups B n , n > 2, F 4 , J 2 (m), m > 6 , m even, (B n , n > 2, G n , n > 3, F 4 and G 2 — / 2 ( 6 ), if one considers Weyl 
groups of simple Lie algebras). Then det (P(p)) has two different irreducible active factors, s(p) and l(p), such that 
det (P(p)) = c 2 s(p) l(p), with s(p(x )) = JX rgK+ K{ x ), an d Kp( x )) = rireTC + , R{x) 2 , respectively, where P +)S and 
Lt +,1 are the sets of positive short roots and of positive long roots of W (more in general, the sets of roots of the 
two root systems, respectively), c 2 is the constant appearing in Eq. (9\). and for any root r, l r (x) is the linear form 
l r (x ) = r ■ x. The X-vectors X^(p) and X^(p), corresponding to the active polynomials s(p) and l(p), respectively, 
that appear in the boundary equation m can be obtained from the following formulas: 


Xi s \p(x)) = 2 

r£K+, s 


S7p a {x) ■ r 
l r (x) 


^a ) iP(x))=2 ^ 

r£K+,l 


VPapc) 

l r (x) 


r 


Va = 1,..., n , 


(14) 


valid Va; £ R n , and one has: 

Ai s) (p) + A«(rt=Ai det(p)) (p), Va = l,...,n. (15) 

Proof. A group of orthogonal matrices, as is W, acts on a vector space by maintaining the vector lengths. Then the 
long roots and the short roots are transformed into themselves by W transformations (each root system is stable with 
respect to W-transformations). This implies that the set of points contained in the union of the reflection hyperplanes 
orthogonal to the long roots is invariant under W transformations, and so does the set of points contained in the 
union of the reflection hyperplanes orthogonal to the short roots. We can then rewrite Eq. ©. using Eq. m, in the 
following way: 


where 


det(P(p(a;))) = c 2 






c 2 s(p(x)) l(p(x )), 


s{p{x)) = l r {x ) 2 , l(p(x))= l r (x ) 2 

r&Tl+,s r£H +t i 


are invariant polynomials, and so they can be written in terms of the basic invariant polynomials p\{x),... ,p n [x). 
The proof of Theorem [T] can then be repeated with s(p(x)) or l(jp(x)) in place of det (P(p(x))), and it gives Eqs. (fill) . 
By comparing Eq. m and Eqs. m, one immediately finds Eq. m- o 


In the hypothesis of Theorem [2] the equations s(p) =0 and l(p) = 0 contain different orbit type strata of S because 
points in the set of reflecting hyperplanes with images in s(p) = 0 have isotropy subgroups that cannot be conjugate 
to any isotropy subgroup of points in the set of reflecting hyperplanes with images in l(p) = 0. 

It is interesting to see how the P-matrices transform when the basic invariant polynomials are transformed. We 
will examine two classes of basis transformations. The first one is related to a change of the canonical basis in R", 
that naturally implies, among others, a change of the reflection hyperplanes, of the matrix representatives of the 
reflection group, and of the basic invariant polynomials. The second one is related to a change of the basic invariant 
polynomials performed without changing the basis of R n , that is, obtained through algebraic transformations of the 
basic invariant polynomials such to maintain their degrees and algebraic independence. 

A change of the (orthonormal) coordinate system used in R”: x —> x' = Rx , R £ 0(n), determines a change of the 
matrix representatives of any matrix group acting in R”. Let’s call W and W' the matrix group before and after the 
transformation. It is well known that W' = {g' \ g' = RgR r , V g £ W}. In fact, if y = gx , y' = Ry, x' = Rx, one has 
y' = Ry = Rgx = RgR T x' = g'x' and by comparison g' = RgR T . The basic invariant polynomials change accordingly 
from the expressions p a (x) to other expressions p' a {x'), so that, Va: £ R", one has: p' a (x') = p a {x), p a {gx) = p a {x), 
V g £ W, and p' a {g'x') = p' a (x'), V g' £ W'. This happens in particular to the matrix elements of P(x), as they are 
invariant polynomials in x: P a b(x) = P' ab (x'). It is convenient to state and prove the following well known theorem, 
because we will use it very often in the following. 


Theorem 3 Let x —> x' = Rx, with R £ 0(n), be an orthogonal coordinate transformation o/R", W a reflection 
group, defined in the old coordinates x, and W' = {g' \ g' = RgR T , Mg £ Wj the corresponding reflection group in 
the new coordinates x'. One has then: 
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1. The W'-invariant polynomials p'{x') are related to the W-invariant polynomials p(x ) by the following equation: 

p'{x') = p(R T x') = p{x). 


2. If pi(x),... ,p n (x), is a basis of W-invariant homogeneous polynomials, then p[(x') = pi(R T x'),... ,p' n (x') = 
p n (R T x'), is a basis of W'-invariant homogeneous polynomials, and if p(x) = p(pi(x ),... ,p n (x)) is a W- 
invariant polynomial, then the IP' -invariant polynomial p'(x') = p(R T x') has the same functional dependence p 
on the new basis, as that one of p(x) on the old basis, that is: 

P{x)=pipi(x),...,p n (x)) <^=> p'{x')=p(p' 1 (x'p' n (x')). 

3. The P-matrices P(p) and P(p'), constructed with the bases p a (x) and p' a {x'), a = 1 respectively, are 

identical, except for the substitution of the coordinates p a with the coordinates p' a , Va = 1 ,,n, that is: 

P(p') = P(p) # . 


Proof. Q] If p(x) is a IP-invariant polynomial, p(gx) = p(x), V g £ W and Vx £ K™. One defines the polynomial 
p'(x') in such a way that p'(x') = p(x), that also implies p'(x') = p(R T x'), because x = P T x'. It is then just to prove 
the bP'-invariance of p'(x'). One has: p(x) = p(gx) = p((R T g'R)(R r x')) = p(R T g'x'), but also p(x) = p(R T x '), so 
that: p(R T g'x') = p(R T x'), Vg' £ IP' and Va/ £ M", that is p'(g'x') = p'(x'), Mg’ £ IP' and Vs' £ R", and the IP' 
invariance of p'{x') is proved. 

[2] With the basis of IP-invariant polynomials p a {x), a = 1 one can form n IP'-invariant polynomials 

Pa(x') = p a (R T x'), Va = 1 These n polynomials are such that p a {x) = p' a (x'), Va = 1 and form 

a basis for the IP'-invariant polynomials. In fact, if p'(x') is a IP'-invariant polynomial, p'(x') = p(x), with p(x) a 
bP-invariant polynomial. But then p'{x’) = p{x) = p(pi(x),... ,p n (x)) = pip'^x'),... ,p' n (x')), and this proves that 
p'{x') can be written as a polynomial of the basic invariant polynomials p' a {x'), a = 1 ,..., n, with the same functional 
dependence p on the new basis, as that one of p(x) on the old basis. 

[3] By substituting in the P-matrix P(p) the variables p a , a = 1 ,...,n, with the basic bP-invariant polyno¬ 
mials p a {x), a = 1 one obtains the matrix P(pi(x),... ,p n (x)) = P{x), whose elements are determined 

by Eq. © and are homogeneous bP-invariant polynomials of x. From item 2, using the variables x ', one has: 
P'(x') = P(R T x') = P(pi(R T x'),... ,p n (R T x')) = P(p' 1 (x r ), ... ,p' n (x')), whose matrix elements are homogeneous 
bP'-invariant polynomials of a/, with the same functional dependence on the new basic invariant polynomials as those 
of P(x) on the old basic invariant polynomials. It remains to show that the matrix P'(x') coincides with the matrix 
P(x') obtained from Eq. ©, using the basic polynomials p'i(x'),... ,p' n (x'). One has, Va, b = 1,..., n: 


Pab{p'l(x'), ■ ■ ■ ,p' n (x')) = Pab{x') = Vp^(x') • S7p' b {x') = ^ 

i= 1 


dp' a {x') dp' b {x') 
dx\ dx\ 


E 


dpajx) dp b (x) 
dx\ dx[ 


^ dp a (x) dRjix'i dp b (x) dRj m x' m 
. dxj dx\ dxk dx[ 


A dp a (x) T dp b {x) T _ dpgjx) dp b (x) _ 

dxj jl dxk kl ^rf dxj jk dxk 

) —-L 3i &—-L 


= E = ^Pa(x) ■ 'Vpb(x) = Pab(x) = Pab(pi(x), . ■ . ,Pn(x)) ■ 

3 = 1 0 3 

The P-matrix P(p') is then obtained from the P-matrix P(p) just by replacing the variables p a with the variables p ' a , 
Va = 1,..., n. o 

The basic invariant polynomials, as well as the matrices representing the group of linear transformations, are 
dependent on the system of coordinates used in M”, but the P-matrix is independent from this system of coordinates. 
The P-matrix depends only on how a given basis of bP-invariant polynomials generates the ring of bP-invariant 
polynomials pR"] w '. For this reason the P-matrix can be considered a basic tool of constructive invariant theory ITR . 
Theorem [3] holds for a general compact linear group G , because the action of G can always be obtained with 
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orthogonal matrices and G possesses a finite basis of homogeneous invariant polynomials. 


From Theorem [3l item[2l we know that under an orthogonal coordinate transformation of M n , x —> x', the basic 
invariant polynomials change: p a (x) -» p' a (x'), Va = 1 so the weighted coordinates of M n , that are then- 

representatives in the space 1" hosting the set S , also change: p a —»• p' a . However, for Theorem [31 item [31 the P-matrix 


changes only for the substitutions p a —> p' a - 


P(p') = 


Pip) 


, and this implies that the set S (with its strati- 

p=p' 


fication), that is determined by the P-matrix, only, does not change under orthogonal coordinate transformation of K". 


Let’s now examine transformations of the basic invariant polynomials obtained without changing the coordinate 
basis in K". One changes the basic set of invariant polynomials: p a (x) —> p' a {x), a = 1 in such a way to 

maintain the homogeneity, the degrees d a , a = 1 ,n, and the algebraic independence of the p' a (x), a = 1 ,... ,n. 
As Va = 1,. ..,n, the p' a (x) are hF-invariant homogeneous polynomials of degrees d a , it is possible to write them as 
polynomials of the basic invariant polynomials p a (x ), Va = 1 ,..., n, in such a way that: 

Pa( x ) =p'a(Pi( x ),---,Pn{x)), Va=l,...,n, V X g IT . 


Removing the dependence on x, this basis transformation can be seen as a particular coordinate transformation of 
p £ M", that is: 

p'a =Pa(Pl>--->Pn) = p'a{p) > V a = 1, . . . , Tl , 

in which p' a (p) is a w-homogeneous polynomial of weight d a of the weighted variables pi,... ,p n . 

The jacobian matrix of this transformation, J(p), has elements: 


J a b{p ) 


dPa(p) 

dpb 


Va, b = 1,..., n , 


(16) 


that are w- homogeneous polynomial functions of weights d a — db of the weighted variables p\,... ,p n , so they are 
vanishing if d a < db, and constant if d a = db- 

The requirement that the transformed invariant polynomials be algebraic independent gives some constraints. One 
can require, for example, that the expression of p' a {p) contains necessarily the variable p a , Va = 1 ,...,n. For 
simplicity, let’s suppose to adopt the convention (O on the labeling the basic invariant polynomials according to their 
degrees. If all the degrees d a are different, that is, in the case of all irreducible finite reflection group different from D n , 
with even n , the diagonal elements J aa {p), Va = 1 ,..., n, must be non-zero constants, and the non-diagonal elements 
Jabip )> with a < b, must vanish, and this implies that J(p) is a lower triangular matrix with non-zero elements along 
its diagonal. In the case of D n , with even n, two degrees are equal and it is possible to have a constant non-singular 
non-diagonal 2 x 2 sub-matrix along the diagonal of J{p), in correspondence with the rows and columns corresponding 
to the two coordinates of weight n, but for the rest J( p ) is a lower triangular matrix, with non zero constant elements 
along its diagonal. The only jacobian matrices that can be accepted as giving rise to transformations between basis 
of (algebraically independent) invariant polynomials must meet all the above requirements. 

This implies that, in any case then, det(J(p)) is a non-zero real constant. The inverse transformation p' —> p is 
then everywhere well defined. The inverse transformation p’ —> p is a IT-homogeneous polynomial map, like the 
transformation p —>■ p ’, that is, p a = p a (p'), Va = 1,..., n, are w-homogeneous polynomial functions of the weighted 
variables p[,... ,p' n , whose weights are w(p a (p')) = d a , and such that p a (p') contains necessarily p' a , Va = 1 ,..., n. 

The transformation p —> p' causes a transformation of the P-matrix and this implies also a transformation of the 
geometric shape of the set S describing the orbit space. Because the maps p —> p’ and p’ —» p are polynomial, this 
transformation of S is a diffeomorphism, and this is enough to know that S maintains in any case its topological 
shape and stratification. 


It is convenient to give the following definition. An allowed basis transformation is a transformation p' = p'{p) with 
Jacobian matrix J(p) that satisfies the following requirements: 

1 . J a b(p) is a re-homogeneous polynomial of weight d a — db] 

2 . det(J(p)) is a non-zero constant. 

Let P a b(p') and P a b(p), Va,b = 1 ,..., n, be the matrix elements of the P-matrices P{p') and P{p), corresponding 
to the bases p[(x),... ,p' n {x), and p\(x),... ,p n (x ), respectively. The relation between the matrices P(p') and P(p) is 
given by the following theorem: 
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Theorem 4 If J(p) = \\dp' a {p) / dpb\\ is the jacobian matrix of an allowed basis transformation p■ = p'(p), whose 
inverse transformation is p = p(p'), the relation between the P-matrices P{p') and P(p) corresponding to the bases 
Pi(x),... ,p’ n (x), and p\(x ),... ,p n (x), is given by the following P-matrix transformation formula: 


Proof. 


Pip') 


J(p) P{p) J T (p) 


p=p(p') 


Pcdip') 


p'=p'(p(*)) 


(Vp'c(x) ■ Vp' d (x)) 


P'(x)=p'(p(x)) 


Vp'c(p(x)) ■ Vp' d (p(x)) = 


(17) 


that is: 


n 

= E 

a,b= 1 


dp'cip) 

dp a 


V Pa(x ) 

p=p(x) 


dp diP) 

dpb 


n 

V Pb{x) = ^2 J ca(p)Pab{p)Jbd(p ) 

P=P( a: ) a,fe=l 


p=p(x) 



=p'(p) 


J{p) P(p) J T (p) 


(18) 


Using the inverse transformation p 


p(p'), this formula can be written in the form of Eq. (1171) . 


o 


It is interesting to see how an active polynomial a(p) (possibly coincident with det(P(p))), and its corresponding 
A-vector A^ (p) transform under transformations of the basis of invariant polynomials. As we will use these results 
in the following we will state and prove the following theorem. 


Theorem 5 Let W be a rank-n irreducible finite reflection group, pi(x),... ,p n (x), a basic set of W -invariant poly¬ 
nomials, and P{p) the corresponding P-matrix. Let a(p) be an active polynomial of P(p) and X^ a \p) its \-vector. 

1. If one performs an orthogonal coordinate transformation x —> x' = Rx, and p' a {x') = p a (R T x'), a = 1,... ,n, 
are the basic invariant polynomials in the new coordinates x', the active polynomial a'(p'), corresponding to a(p) 
in the new coordinates x', and its X-vector A^“ \p'), are the following: 

a'(p') = a (p)\ P =p' > A (a ' } (p') = A (a) (p) 

p=p' 

2. If one performs an allowed basis transformation p' = p'(p), with Jacobian matrix J(p) = \\dp' a (p)/dpb\\, the 
active polynomial a'(p'), corresponding to a{p) in the new basis of invariant polynomials p'(x), and its X-vector 
X (o 'V), are the following: 


a '(p') = a{p)\ p =p(p') , 


X b a \p') = E J bc{p)^i\p) 

C—l 


p=p(p') 


Mb = 1 ,n. 


(19) 


Proof. [H Let a(jp) be irreducible (in the reals or in the complexes). Then a{p) is a factor of det(P(p)), and satisfies 
the boundary equation ED- Theorem [3] item [3] establishes that in the new coordinates x' = Rx, the P-matrix 
is unchanged in form, only one has to replace the old variables p with the new ones p'. Then det(P(p')) has the 
irreducible factor a'{p') = a(jp) | ,, that satisfies the boundary equation (fill) , in which all variables p a have been 

replaced by the variables p' a , Va = 1,... ,n. Then a'[p') is active and its A-vector A^ a \p') is obtained from X^ a \p) 
just by replacing the variables p a with the variables p' a , Va = 1 ,n, that is A^ a \p') = A^ a ^(p)| p p/ . If a(p) is a 

reducible active polynomial, it is the product of irreducible active factors of det(P(p)), and its A-vector is the sum of 
the A-vectors of its irreducible factors. Then, a'(jp') = a{p) | , is a reducible polynomial in the variables p' , that is 

the product of irreducible active factors of det(P(p')), and so it is an active polynomial, and its A-vector is the sum 
of the A-vectors of its irreducible active factors. 

[2] Let a(p) be irreducible (in the reals or in the complexes). Theorem [I] establishes that P{p') = 
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J(p) P(p) J T (p) 


p=p(p') 


Then det(P(p')) = j 2 clet(P(p)) 


p-p(p') 


where j = det(J(p)) is a non-zero real constant, 


and a'(p') = a(jp )| p= ^( p /) must be an irreducible factor of det(P(p')). It results that a'(p') is an active polynomial 
because it satisfies the boundary equation CD- In fact, V a = 1,..., n, we have, at the first member of Eq. CD, 
written in the new coordinates p', and understanding the summation symbols: 


p , /\ da\p') 

Pbc{p } ~d^r 


( JbdiP ) PdeiP) Jjjp)) 


( da (p)\ 

p=p(p') \ dpf J 


p=p(p') 


dp/(p ') 
dp' c 


( Jbdip ) Pde{p )) 


dp f (p') (dfi c (p) 


p=p(p') 


dp’c 


dp e 


fda(p)\ 

p=p(p') V dpf J 


p=p(p') 


( Jbd{.P) Pde(p)) 


P=P(P') 


s„ (AM 


\ dpf 


— Jbdijp) L = ; 


P=P(P') 


P=P(P') 


^) 


P-P(p') 


J M| p=f(p0 (\ { ;\p)a(p)) | p=?(p() = (^(p)A^(p))[ = ^ o a'{p'). 


p-p(p') 


By comparing with Eq. CD, written for the variables p ', and a'{p') 1 we see that a'(p') is an active polynomial and 

Mb = 1 ,..., n. 


that its A-vector has the elements A^“ \p') given by the equation: A^“ \p') = ^ Jbd(p ) A ^(p)') 


p=p(p') 


Many of the definitions and results given in this Section can be extended to a general compact linear group, not 
necessarily a finite reflection group, because its algebra of invariant polynomials admits a basis made by a finite 
number of homogeneous invariant polynomials. In particular, the P-matrix can be defined for all these groups. 

As already recalled in the introduction, the construction of the P-matrix from the definition, that is using Eqs. © 
and m, when the rank n of the group is large, requires an enormous effort, computer memory and computing time, 
and is often impossible to take to completion. The generating formulas (l25ll . ([32]) , (l35l) and CUll . obtained in Sections 
mm for a given choice of the basic invariant polynomials, are then very useful, especially because, using Eq. dm 
they allow to determine the explicit form of the P-matrices in all bases of invariant polynomials of the groups of type 
S n , A n , B n and D n . Examples [5] SI in Section m show the use of the P-matrix-transformation formula, Eq. CD, 
and the A -vector transformation formula, Eq. C9- 

It is now possible to recall the results of Refs. [39j], [U| and j2l| with some precision. All of them are concerned 
with the determination of distinguished bases of invariant polynomials satisfying some supplementary conditions. 

In Ref. [H Saito, Yano and Sekiguchi proved that, for all irreducible finite reflection groups, it is always possible 
to choose the basis of homogeneous invariant polynomials p\{x),... ,p n {x), in such a way that the derivative of the 
P-matrix P(p), made with respect to the highest weight variable ph , is a constant matrix, that is, it is always possible 
to have 

A a b = = a real number, V a, b = 1 ,n. ( 20 ) 

dph 

Moreover, given two bases for which this property holds, their basic invariant polynomials may differ only for scale 
constant factors. This means that, if all the basic polynomials and the corresponding variables pi,... ,p n are ordered 
according to their degrees, or weights, (so ph = p n ), and because the degrees of the basic invariant polynomials satisfy 
the relations d a + d n _ 0+ i = d n + 2 , Va = 1 ,..., n, it is possible to choose the basic invariant polynomials in such a 
way that p n appears in P{p) only along the auxiliary diagonal connecting the elements Pi n (p) and P n i(p)- In |39] 
these bases were called flat, because they establish a flat structure on the orbit space, that is a system of coordinates 
with a flat metric. The authors of Ref. [3f| also determined explicitly the flat bases of invariant polynomials of 
all irreducible finite reflection groups, except E 7 and E s (and they also calculated explicitly the P-matrices of the 
groups E e , Pj, G 2 , H 3 , H±, I 2 (to), m > 2). The flat bases of invariant polynomials for E 7 and E%, together with the 
corresponding P-matrices were found in [49| . Ex. [2] in Section IIXI shows how one can determine the P-matrix and 
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the A-vector of its determinant in a flat basis of A 3 , starting from the P-matrix and the A-vector obtained with the 
generating formulas here given for the groups of type A„, and using Theorems 0] and [5] 

In Ref. [4d| Sartori and Talamini proved that, in the case of a compact linear orthogonal group having a basis of 
algebraically independent invariant polynomials, with pi(x) = ||a ;|| 2 (and this excludes groups of type S n ), it is always 
possible to choose the basis of invariant polynomials in such a way that the only non zero element of the A-vector 
A ( a \p), of any active polynomial a(p), is the first one, that is: 

X[ a) {p) =2w(a), X ( b a \ p ) = 0, V6 = 2, (21) 

Moreover, given two bases for which this property holds, their basic invariant polynomials are related by an allowed 
basis transformation not dependent on pi. In Ref. [4lJ these bases were called a-bases. Ex. [3] in Section HXl shows 
how one can determine the P-matrix P(p) and the A-vector of its determinant ,\( det ( p ))(p) in a det(P)-basis (that is 
in an a-basis in which the active polynomial a(p) is det(P(p))) of A 3 , starting from the P-matrix and the A-vector 
obtained with the generating formulas here given for the groups of type A„, and using Theorems 0] and 0J 

In Ref. plj | Flatto defined a bilinear map that maps two polynomials p{x) and q(x) into a polynomial (p,q)(x), 
obtained in the following way: 


{p,q){x) =p(V)q(x), 

where p(V) means that in p(x) one has to replace all occurrences of x £ with ^ 77 , Vi = l,...,n, VfcgN, and p(V) q{x) 

is the polynomial obtained by the application of the differential operator p(V) to the polynomial q(x). If p{x) and 
q{x) are homogeneous of degree d p and d q , respectively, then ( p,q)(x ) is 0 if d p > d q , a real number if d p = d qi a 
homogeneous polynomial of degree d q — d p if d p < d q . This also shows that the map (•, •) is non-symmetric in the two 
arguments. In Ref. [ 2 ]] Flatto proved that, for all effective irreducible finite reflection groups W, besides multiples, 
there exist a unique basis of homogeneous invariant polynomials, pi{x),... ,p n (x ), if one requires that p\(x) = ||a:|| 2 
and that all basic invariant polynomials satisfy the conditions: 

{Pa,Pb)(x) = 0 , Va, b = 1,..., n, such that d a < db , (22) 

(p a ,Pb){x) = (p b ,p a )(x) = 0, if d a = db = n , a^b, (23) 

where the second one concerns the case W = D n , even n, only. In addition to the conditions (l22l) and (l23l) . Iwasaki 
f2fj ] proposed to fix the normalization of the basic polynomials, by requiring ( p a ,p a ){x ) = 1, Va = 1 ,...,n, and 
called canonical such a basis. However, both the basic invariant polynomials and the P-matrices of the canonical 
bases (according to Iwasaki’s definition) contain irrational coefficients, and pi(x) has not the standard form d3j). For 
this reason, I decided to leave undetermined the normalization of the basic invariant polynomials, and call canonical 
any basis of invariant polynomial satisfying Eqs. (l22l) and (l23l) . Two canonical bases differ then only for scalar 
multiples of the basic polynomials. Because of the special form of pi(x), all polynomials of degrees greater than 2 in 
a canonical basis satisfy the differential equation {pi,p a )(x) = V 2 p 0 (x) = 0, Va = 2,..., n, and are thus harmonic 
polynomials. The canonical bases of invariant polynomials of the reflection groups of types A„, B n , D n and / 2 (m), 
were determined in Ref. jUj and those of the reflection groups iJ 3 , P 4 and P 4 were determined in Ref. [3(|. Ex. 
0] in Section HXl shows how one can determine the P-matrices and the A-vectors in a canonical basis of A 3 , starting 
from the P-matrix and the A-vector obtained with the generating formulas here given for the groups of type A n , and 
using Theorems 0] and [5] 

I calculated the P-matrices of the group Bs for the following bases: the basis given in Section |Vl a flat basis, a 
det(P)-basis and a canonical basis. The form and complexity of the P-matrix depends strongly on the basis. This can 
be seen from the matrices reported here under, that, for each of these 4 choices of the basis, they have in their matrix 
elements the number of terms of the w-homogeneous polynomials forming the corresponding matrix elements of the 
P-matrix. All flat bases differ only for scalar multiples of the basic invariant polynomials, and so do all canonical 
bases. This implies that the number of terms in the matrix elements of the P-matrices of all flat bases are the same 
and that the number of terms in the matrix elements of the P-matrices of all canonical bases are the same. On the 
contrary, the number of terms in the matrix elements of the P-matrices of the a-bases depend on the a-basis chosen. 
For example, there are det(P)-bases of B$ in which the element Pgs{p) has 146 terms, the maximum number of terms 
that is possible for a homogeneous invariant polynomial of degree d 3 + d 3 — 2 = 30, when expressed in terms of basic 
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invariant polynomials of B$. The det(P)-basis used to obtain the matrix written above, was obtained from the basis 
defined in Section |Vl by avoiding, in the basis transformation, all terms independent from p-j. Anyway, this calculation 
seems to suggest that one of the simplest form of the P-matrix is obtained with the bases here proposed. 


Number of terms in the matrix elements of the P-matrices of B% 
in correspondence with different bases of invariant polynomials. 


basis of Section lYl 


flat basis 
11111 


1 

2 

2 

2 

2 

2 

2 

1 


1 

3 

4 

5 

7 

8 

10 

10 

1 

2 

3 

3 

3 

3 

2 

1 


1 

4 

7 

10 

13 

18 

21 

25 

1 

2 

3 

4 

4 

3 

2 

1 


1 

5 

10 

15 

21 

26 

33 

39 

1 

2 

3 

4 

4 

3 

2 

1 


1 

7 

13 

21 

28 

36 

44 

54 

1 

2 

3 

3 

3 

3 

2 

1 


1 

8 

18 

26 

36 

46 

57 

66 

1 

2 

2 

2 

2 

2 

2 

1 


1 

10 

21 

33 

44 

57 

68 

79 

\1 

1 

1 

1 

1 

1 

1 

1 J 



10 

25 

39 

54 

66 

79 

88 / 


det(P)-basis 


canonical basis 


( 

1 

1 

1 

1 

1 

1 

1 

1 \ 


/I 

1 

1 

1 

1 

1 

1 

1 \ 


1 

3 

4 

5 

5 

4 

5 

4 


1 

3 

5 

7 

11 

15 

22 

29 


1 

4 

7 

7 

9 

9 

8 

7 


1 

5 

7 

11 

15 

22 

29 

40 


1 

5 

7 

12 

14 

14 

13 

11 


1 

7 

11 

15 

22 

29 

40 

52 


1 

5 

9 

14 

19 

19 

18 

16 


1 

11 

15 

22 

29 

40 

52 

70 


1 

4 

9 

14 

19 

26 

25 

22 


1 

15 

22 

29 

40 

52 

70 

89 


1 

5 

8 

13 

18 

25 

31 

27 


1 

22 

29 

40 

52 

70 

89 

116 

V 

1 

4 

7 

11 

16 

22 

27 

33 / 


u 

29 

40 

52 

70 

05 

00 

116 

146 / 


III. GENERATING FORMULAS FOR THE SYMMETRIC GROUPS S„ 

The symmetric groups S n are not irreducible finite reflection group but they are strictly connected with the 
irreducible finite reflection groups of type A n , and we will use them to obtain results for the groups A n . 

Let ei,...,e n be the unit vectors of a canonical basis of R n , and xi,...,x n be the corresponding coordinates of 
a vector x G R", so that: x = Y^i=i x i e i- Sn acts on R™ by permuting in all possible ways the n coordinates 
of the vectors x G R". Its order is n\. S n can be generated by the Q) transpositions Uj, 1 < i < j < n, 
that exchange two coordinates Xi and Xj leaving all the others unchanged. The transpositions iy are reflections 
of the space R": Uj corresponds to the reflection about the hyperplane of equation Xi — Xj = 0. Defining 
the vector a = e; — ej (or any non zero multiple of it), the action of t,j on any vector x G R" is given by Eq. 
©. Hence, S n is a finite reflection group. The positive roots are for example the Q) vectors e* — ej, V i < j = 1,..., n. 

A possible choice of the basis of invariant polynomials of S n , is obtained by using the n elementary symmetric 
polynomials: 


pi(x) = si(xi,. 
p 2 {x) = S 2 (x 1 , ... ,x n ) 
p 3 (x) = s 3 (x 1 , .. .,x„) 


Xi + ... + x n , 

Xix 2 + Xix 3 + . . . + X n -iX n , 

XiX 2 X 3 + XiX 2 X4 + . . . + X n - 2 X n -iX n , 


(24) 


Pn{x) = S n (x l,...,x n ) = X\X 2 ■ ■ ■ X n . 


The degrees of the basic invariant polynomials are 1, 2,..., n. It is clear that the elementary symmetric polynomials 
s a (x). a = 1 ,..., n, are invariant for all permutations of the coordinates x ±,..., x n of x. 

Theorem 6 1. The P-matrix corresponding to the basis of invariant polynomials of S n given in Eq. has 
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matrix elements that can be obtained from the following generating formula: 

min (a,b) 

Pabip) = [n+l-m.m(a,b)]p a -ip b -i- ^ {a + b- 2i) pt-i p a+& _i_;, Vo, b = 1,..., n, (25) 

i=max(l,a+6-n-l) 

in which one has to consider po = 1 . 


,3. T/ie determinant of the P-matrix P(p) given by Eq. 
X-vector: 


satisfies the boundary equation m with the following 


A (de t (P)) (p ) = ( A ( det ( P »(p), . . . , A ^\ P )) , 

^i det(i 3 )) (p) = ~{n- a + 2)(n - a + l)p a -2 , Va=l,...,n, (26) 

m which one has to consider p-\ = 0 and po = 1. 

Proof. In Section IX A I o 

Using Eq. ©, we have a formula to express the discriminant A (a:), that is the lowest degree invariant polynomial 
vanishing in the set of the reflecting hyperplanes of S n , in terms of the elementary symmetric polynomials: 

n 

A (a:) = {x t - Xj) 2 = det (P(p(x ))), 
i<j =i 

in which P(p) is the P-matrix obtained with the generating formula ( 1251 ) . and in which the variables p a , a = 1 ,..., n, 
have been substituted with the elementary symmetric polynomials in Eq. ( 1241 ) . A formula of this kind was, for 
example, invoked by Artin in Ref. Q, p. 548. 


IV. GENERATING FORMULAS FOR THE GROUPS OF TYPE A n 

The symmetric group S n +i acts in R ra+1 by permuting the n + 1 variables x±, ..., a; n +i> in all possible ways. This 
action is completely reducible. The hyperplane H having equation x\ + X 2 + ■ ■ ■ + x n +i = 0 is invariant under this 
action (and so is the line orthogonal to it). The group A n is the group acting in the n-dimensional subspace H of 
R ” +1 in the same way as S n +\. As is a finite reflection group, A n is a finite reflection group. The positive 
roots of A n , as vectors of R n+1 , have the same expressions as the positive roots of S n + 1 , and are, for example, the 
("J 1 ) vec t° rs e i~ e ji V i < j = 1,... ,n + 1. The order of A n is (n + 1)!, the same as that one of S n +i- The basic 
invariant polynomials of A n for this action are the basic invariant polynomials of S n + 1 , of degrees 2 , 3,..., n + 1, 
for example those obtained from Eq. (1241) . with n replaced by n + 1. One can verify that Eqs. © and m are 
satisfied by these degrees. These invariant polynomials, are however expressed in terms of n + 1 variables, and 
one wants to write them in terms of only the n variables of the n-dimensional subspace H C R ™ +1 where A n 
acts effectively. To this goal one can consider the unit vector u n+ \ = (1,1,..., 1 )/y/n + 1 € R” +1 , orthogonal to 
the invariant hyperplane H and a rotation matrix R n + i € 0(n + 1) that transforms u n+ \ into the canonical unit 
vector e n +i. As a consequence of this rotation, the invariant subspace H is mapped to the invariant subspace of 
R " +1 with equation x n +i — 0, and the invariant polynomials forming a basis of S n +1 (and of A n ) change form, 
according to Theorem [3] item [ 2 ] At this point if one puts x n+ i = 0 in the expressions of the “rotated” basic 
invariant polynomials of S n +i, the linear invariant becomes identically zero in R" +1 , and the others become the n 
basic invariant polynomials of A n expressed in terms of the n variables xi,... ,x n , only. A renaming of the basic 
invariant polynomials of A n obtained in this way is also convenient because one usually wants the index i, that 
labels the basic invariant polynomials of A n , start at 1 and end at n. As the basic invariant polynomials of A n are 
obtained from those of S n + 1 , one expects that also the P-matrix of these basic invariant polynomials of A n can 
be obtained from that one of S n +±, but this is not at all straightforward, as it can be seen from the proof of Theorem[7] 

To put all this in a more formal way, also convenient for the proofs, let’s introduce the following notation. A point 
of Rn+i w jq indicated with x, and its projection in R”, obtained by taking = 0, with x. The basic invariant 
polynomials of S n +i, obtained from Eq. (IM1) . with n replaced by n + 1, are indicated with s a (x), a = 1,..., n + 1. In 
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the rotated coordinates x' = R n +±x, because of Theorem^ item[ 2 ], the basic polynomials s a (x), a = 1,..., n + 1, of 
become s' a {x') = s a (R / l+1 x'), and, in particular, s^fV) oc x' n+l . Then, by taking x ' n+1 = 0, that corresponds to 
take s' 1 (x l ) = 0, one finds n basic invariant polynomials of A n : s' a+ 1 (x'), Vo = l,...,n, whose degrees are d a = a + 1, 
Va = 1,... ,n. There are however infinitely many ways to define a set of basic invariant polynomials for the group 
A n . We require two main things. The first one is that we want the standard form of the quadratic invariant, that 
is, we want that pi(x') = x 'i ■ The second one is that we want that the P- matrix, that is originating from the 

basic invariant polynomials of A n , has only integer coefficients. This will strongly reduce the possible choices of the 
basic invariant polynomials of A n . One of these choices is the one presented in the following, in items 1., 2., 3. and 
4.. 

Let’s first define the matrix R n +i £ 0(n + 1) that transforms u n + 1 into e n +i, that is used in the proof of Theorem 
[7] In the definition below, it is written R in place of R n + i to leave free space for the indices. 


Rn +1 — R • 


Rk,i 

— 1 / y/k{k + 1 ) , 

1 < k < n , 

1 < i < k, 

Rk,k -\-1 

= -k/y/k{k + l), 

1 < k < n , 


Rk,i 

= 0 , 

1 < k < n , 

k + l<i<n + l, 


= 1 /y/n + 1 , 


1 < i < n + 1 . 


(27) 


To say it in words, the /c-th row of R n +i (with 1 < k < n) is the vector that has the first k elements equal to 1, 
the next element equal to — k and all other elements equal to 0 , all multiplied by a normalization factor equal to 
1 /y/k + fc 2 3 , and the last row of R n +i, the (n + l)-th row, has all its n + 1 elements equal to 1/y/n + 1. 

By construction all row vectors are mutually orthogonal and normalized to 1, and this is sufficient to say that the 
matrix R n +i is orthogonal. An easy calculation shows that its determinant is 1, so it is a proper rotation matrix. To 
see this substitute column 2 with the sum of column 2 and column 1, substitute column 3 with the sum of column 3 
and the new column 2, and so on. The determinant does not change for these substitutions and the resulting matrix 
is a lower triangular matrix with its diagonal elements equal to d *. = k/y/k 2 + k = y/k/{k + 1), V/c = 1 ,,n, and 
d n+ 1 = (n + 1 )/y/n + 1 = y/n + 1 , so its determinant is equal to nL-=i = V n + 1 TIfc=i \fk/(k + 1 ) = 1 . 


We can now define the basic invariant polynomials of A n we shall use in Theorem [7] below. 


1. Start from the n + 1 elementary symmetric polynomials si(x),..., s„+i(x), obtained from Eq. (l24| . with n 
replaced by n + 1, that form a basis for the action of the group S n +i in R" +1 : 

Sl(x) = X1 + ... +X n+ 1, 

S 2 (x) = X1X2 + XyX 3 + . . . + X n X n+ i, 

S 3 (x) = X 3 X 2 X 3 + X1X2X4 + . . . + X n -lX n X n+ i, (28) 


Sn+l(x) 


X1X2 ‘ ’ ' X n -\-\. 


2 . Rotate the system of coordinates with the matrix R n +i, that is x' = R n +1 x. Theorem[3j item [71 states that the 
basic invariant polynomials of S n +1 in the rotated system of reference have the following expressions: 

s' a (x') = s a {Rl + 1 x'), a = 1 ,..., n + 1. (29) 

With this transformation the polynomial .$( (x r ) becomes proportional to x' n+1 . Any orthogonal matrix R n +1 
transforming u n +1 into e ra +i can be used. In the proofs we shall use the matrix R n +\ defined in Eq. (12711 . All the 
other orthogonal matrices transforming u n+ 1 into e n+ i differ from the given R n +i by orthogonal transformations 
mixing the first n variables, only. Theorem 0 item [3], ensures that the P-matrix corresponding to the basic invariant 
polynomials s' a (x'), a = 1,... ,n + 1, has the same form of the P-matrix (given by Eq. (l25l) f corresponding to the 
basic polynomials s a (x), a = 1,..., n + 1, and differs only for the substitutions s a —> s' a , a = 1,..., n + 1. It is clear 
from Theorem [3] that any other orthogonal matrix of order n + 1 transforming u n +i into e n +i would give different 
expressions for the basic invariant polynomials but the same form of the P-matrix. 

3. Perform a scale transformation on the n- 1-1 basic invariant polynomials (x 1 ),..., s' n+l {x') in the following way: 

t a (x') = -2 y/{n+ l ) a_2 s' a (x '), a = 1,... ,n + 1. (30) 
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The transformations (l29l) and m can be done in the reverse order without changing the resulting polynomials 
t a {x'), Va = 1,..., n. The scale transformation (1301) is necessary only if one wants the quadratic invariant polynomial 
of the standard form © and a P-matrix with only integer coefficients. 

4. Put x ' n+1 = 0 in the arguments of the basic polynomials t a (x'), Va = 1,..., n + 1, that is, take the restrictions 
of the basic polynomials t a (x'), Va = 1,... ,n + 1, to the hyperplane x' n+l = 0. One finds t\{x') = t\{x')\ x , ^_ 0 = 0, 

and the remaining n basic invariant polynomials, t^ix'),...,tn+ i(x , ) ) that are expressed in terms of the n variables 
in x ', only, form the n basic invariant polynomials of A n . Rename and renumber these basic invariant polynomials of 
A n in the following way: 


Pa(x') = t a +l(x') = t a +l(x')\ x ^ i=0 , Va = l,..., 7 l. (31) 

In Example 1, in Section llXl it is shown how to write the explicit form of the matrix P 4 , and how one can determine 
the basic invariant polynomials of A 3 from those of 1 S 4 , following the rules written above. 


Theorem 7 1. The P-matrix corresponding to the basis of invariant polynomials of A n obtained by Eqs. 128 1) \31\) . 

has matrix elements that can be obtained from the following generating formula: 

Pab{p ) = [(n + 1 ) min (a,b) - ab) p a - 1 p b - 1 + 2 (a + b)p a+b _ 1 + 

min(a,6) — 1 

-(n + l) {a + b - 2i) pi_x pa+b-i-i, Va, b = 1,..., n , (32) 

i=max(2,a+6- n— 1) 

where po = 0 and p a = 0, V a > n. The matrix P{p) obtained with Eq. \32\) is symmetric. 


2. The determinant of the P-matrix P(p), given by Eq. 132)) . satisfies the boundary equation U2\) with the following 
X-vector: 


A( det ( p ))(p) = (2 n (n + 1), X^ et{P) \p), ..., X^\p)), 


A { a det(P) \p) 

in which one has to consider po 
Proof. In Section IX B I 


= — (n + l)(n - a + 2 ){n - a + l)p a -2 , 
= 0 (that implies / \( det l p ))(p) =Qj. 


Va = 2,...,n, 


(33) 


o 


Using Eq. ©, we have a formula to express the discriminant A(x), that is the lowest degree invariant polynomial 
vanishing in the set of the reflection hyperplanes of A n , in terms of the basic invariant polynomials of A n : 

A(x) = det(P(p(x))), 

in which P(p) is the P-matrix obtained with the generating formula (1321) . and in which the variables p a , a = 1,..., n, 
have been substituted with the basic invariant polynomials pi(x),... ,p n (x), of A n , obtained in Eqs. (l28l) (l3TT) . The 
polynomial A(x), as calculated above, is unique, modulo constant factors. 

A generating formula like Eq. (15121) was written without proof in Remark 3.9 of Ref. Q, for a basis q a (x), 
a = 1 ,... ,n, that differs from the basis p a {x ), a = 1 ,..., n, defined above, only for scalar factors, precisely: 


q a (x) = (-1)“ 2 (n + 1) a 2 p a (x), Va=l,...,n. 


An equation that differ very little from Eq. (l33l) . was written in Eq. (2.1.5.2) of Ref. J52|, and proved in a different 
manner. 
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V. GENERATING FORMULAS FOR THE GROUPS OF TYPE B n 

Let ei,..., e ra be the unit vectors of a canonical basis of 1". A set of simple roots for B n is the following one: 

Oii = Ci 6i-\-1 , V f = 1, . . . , 77. 1, Oin = Cn . 

A set of simple roots for C n differs from that one for B n only for the length of a n : instead to be a factor \/2 
shorter than the others, is a factor \/2 longer than the others. One can then take the last root of C n equal to 
a n = 2e n . The n simple roots determine n reflections with the formula ©• The n reflections obtained from the 
simple roots of C n are the same as those obtained from the simple roots of B n . By multiplying in all possible 
ways these n reflections one obtains the reflection group B n . Then, the simple roots of B n and the simple 
roots of C n determine the same finite reflection group B n . The reflections determined by the simple roots on, 
V* = l,...,n — 1, when applied to a general vector x € M" exchange the coordinates Xi and Xi+i of x, leaving 

all the other coordinates unchanged, while the reflection determined by a n changes the sign of the coordinate 

x n of x. B n then acts on ffi™ by permuting in all possible ways the n variables x\,... ,x n and by changing 
in all possible ways their signs. The order of B n is then ord(i? n ) = n\ 2". One of the possibilities is to consider 
the vectors e*, i = 1 ,..., n, as the positive short roots, and the vectors e j ±6 j, 1 < i < j < n as the positive long roots. 

A possible choice of the basis of invariant polynomials of B n , already proposed by Coxeter in [18], is obtained by 
using the n elementary symmetric polynomials in the variables x |, that is: 

Pi (x) = si x 2 n ) = x\ + ...+x 2 n , 

p 2 (x) = s 2 {x\,...,x 2 n ) = xlxl + x\x\ + ... + x 2 n _ lX 2 n , 

( \ / 2 2\ 2 2 2 , 2 2 2, i2 2 2 (o a\ 

P3\%) , % n ) H” • • • H - ^‘n—2'^n—l'^n'i 

/\ (2 2\ 22 2 
Pn{X) = S n ( X 1 ,...,X n ) = X X X 2 ■ ■ ■ X n . 

The degrees of the basic invariant polynomials are then 2,4,..., 2 n, and one can verify they satisfy Eqs. © and ©. 
It is also clear that the polynomials p a (x), a = 1,... ,n, defined in Eq. (l34l) . are invariant for permutations and sign 
changes of the variables. 

Theorem 8 1. The P-matrix corresponding to the basis of invariant polynomials of B n given in Eq. has 

matrix elements that can be obtained from the following generating formula: 

min(a,6) — 1 

Pabip) = ^2 4 (a + b - 1 - 2z) pi p a+6 _i_i, Va, b = 1,... ,n, (35) 

i=max(0,a+6— n— 1) 

in which one has to consider po = 1. The matrix P(p) obtained with Eq. (,9,51) is symmetric. 

2. The determinant det (P(p)) of the P-matrix P(p) obtained from Eq. fgdl) has two active factors. One active 
factor is the polynomial s(p) = p n , representing the highest degree basic invariant polynomial p n (x), as well as the 
invariant polynomial square of the product of the n linear forms defining the reflection hyperplanes corresponding 
to the short roots. The other active factor is the polynomial l(p), representing the invariant polynomial square 
of the product of the 2 (”) linear forms defining the reflection hyperplanes corresponding to the long roots. 
s(p) (= p n ) satisfies the boundary equation (ill) with the following X-vector: 

A (a) (p) = (A^ } (p),...,AW(p)), X ( a s) (p) =4(n-a + l)p a -i, Vo = l,...,n, (36) 

l(p) satisfies the boundary equation (111) with the following X-vector: 

X W {p) = (X { 2(p),...,X ( n\p)), AW(p) = 4(n- a+ l)(n-a)p a -i, Va = l,...,n, (37) 

det(P(p)) satisfies the boundary equation 11121 ) with the following X-vector: 

A (det(P)) (p ) = ( A ( det ( p »(p),..., X^\p)), X i det ( p ))(p) = 4 (n - a + 1 ) 2 p a -r, 

where, in Eqs. P?| ) and 138 1) . one has to consider po = 1. 


Va = 1,..., n , (38) 
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Proof. In Section IX Cl o 

Using Eq. ©, we have a formula to express the discriminant A (a"), that is the lowest degree invariant polynomial 
vanishing in the set of the reflection hyperplanes of B n , in terms of the basic invariant polynomials of B n : 

n n 

A (a;) = II *i II “ x f) 2 = s (^(*)) Kp(x)) = 4~ n det (P(p(x))), 

i=l i<j =1 

in which P{p) is the matrix obtained with the generating formula (1351) . and in which the variables p a , a = 1 ,n 
have been substituted with the basic invariant polynomials in Eq. (ESP - Here, s(p(x)) (= p n (x)) and l(p(x )) are the 
two irreducible factors of det (P(p(x))). 

Eq. (1351) was first written in Eq. (**) of Ref. [24j and proved in different manner, using techniques common in 
singularity theory. 

An equation that differs very little from Eq. (l38l) was first written in Eq. (2.2.5.2) of Ref. [52] and proved in a 
different manner. 


VI. GENERATING FORMULAS FOR THE GROUPS OF TYPE D n 


Let ei, ..., e n be the unit vectors of a canonical basis of R”. A set of simple roots for D n is the following one: 


(x% — ei ej+i, Vi — 1 ,..., n 1 , tx n — e n —i T e n . 


The given set of simple roots of D n differs from that one of B n only for the last root a n . The n simple roots of D n 
determine n reflections using Eq. ©. It is clear that the n—1 reflections obtained from the simple roots aq,..., a n -i 
are the same as those of B ni and only the reflection relative to a n differs from the corresponding one of B n . The 
reflection determined by the simple root aj,Vi = l,...,n — 1, when applied to a general vector x £ R” exchanges the 
coordinates Xi and Xj+i of x, leaving all other coordinates unchanged, while the reflection determined by a n exchanges 
the coordinates x n -\ and x n of x and changes their signs, leaving all other coordinates unchanged. By multiplying in 
all possible ways these n reflections one obtains the reflection group D n . It acts on R™ by permuting in all possible 
ways the n variables xi,...,x n and by changing in all possible ways an even number of their signs. The order 
of D n is then ord(Z) n ) = n\ 2 n ~ 1 . The 2(™) positive roots can be taken, for example, to be e,;±e 7 , Vi,j = l,...,n,i<j. 


A possible choice of the basis of invariant polynomials of D n , already proposed by Coxeter in []j§|, is obtained by 
using the n—1 lowest degrees elementary symmetric polynomials in the variables x?, and the elementary symmetric 
polynomial of degree n in the variables Xi , that is: 


Pi{x) = si(xl,...,xl) 
p 2 {x) =s 2 (x l,...,xl) 
P3{x) =s 3 {x\,...,x 2 n ) 


xl + ...+xl, 

x\x\ + x\xl + . . . + X 2 n _ lX 2 n , 


222 , 222 , ,2 2 2 

A *1*2*4 + • • • + x n _ 2 X n _ 1 X n , 


(39) 


Pn — 1(*) — &n— l(*l? • ■ ■ : *n) 
Pn{x) = S n (x 1 , ...,X n ) 


2 2 2 , , 2 2 
*1*2 ' ' ‘ *n—1 T ■ • • T ^2*3 

2 : 12:2 • • • X n . 


X 


2 

n' 


The degrees of the basic invariant polynomials are then 2,4,..., 2(n — 1), n, and one can verify they satisfy Eqs. (01) 
and 0]). If n is even there are then 2 different basic polynomials of degree n: p n / 2 {x) and p n {x), using the labeling in 
Eq. (1391) . It is clear that the polynomials p a {x), a = 1,... ,n, are invariant for permutations of the variables and for 
an even number of sign changes of the variables. The labeling of the basic invariant polynomials given in Eq. ([39D is 
not the standard one based on the degrees, but it turns out to be more convenient for what we are going to do. 
These basic invariant polynomials differ from those of B n , given in Eq. (EH1) . only for the last one, but not very much: 
the square of p n {x) of D n is the basic invariant p n (x) of B n , that is: Pn(x) = ( p® ( x )) 2 , using the notation, also used 
in the proof, with the indices B and D telling us the group to which p n {x) refers. 
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Theorem 9 1. The P-matrix corresponding to the basis of invariant polynomials of D n , given in Eq. m, has 

matrix elements that can be obtained from the following generating formula: 


Pab(p) = £™max(0,a+b-n-l) 4 (« + b ~ 1 “ 2 *) Pi Pa+b-l-i 
Pan(p) = Pna(p) = 2 (n - a + 1) £> a _i p n 
Pnn(jP ) — Pn— 1 i 


Pn—^Pn 


V a, 6 = 1,..., n — 1, 
Va = — 1, 


(40) 


in which one has to consider po = 1. The matrix P(p) obtained with Eq. is symmetric. 


2. The determinant of the P-matrix P(p), given by Eq. m satisfies the boundary equation 112 \) with the following 
X-vector: 


x^ p) \p) = (X™\p),...,X^ P) \p),0), 


A (d et (^))(p) _ 4 ( n _ a + i)( n _ a)p a -! , Va = l, — 1. 
in which one has to consider po = 1. 


(41) 


Proof. In Section IX Dl o 

Using Eq. we have a formula to express the discriminant A (a;), that is the lowest degree invariant polynomial 
vanishing in the set of the reflection hyperplanes of D n , in terms of the basic invariant polynomials of D n : 

n 

A (x) = Rtf- x ^ 2 = 41 n det (P(p( x ))) > 

i<j=1 

in which P(p) is the matrix obtained with the generating formula (1401) . and in which the variables p a , a = 1 ,n 
have been substituted with the basic invariant polynomials in Eq. (EU). 

Eq. (HUl) was first written in Ref. |2J|, at the end of p. 88. It was obtained from the analogous expression for the 
groups of type B n , Eq. (**) of Ref. (|2J|), in essentially the same way of how Eq. (l40l) is here obtained from Eq. 

©• 

An equation that differs very little from Eq. (l4lT) was first written in Eq. (2.3.5.2) of Ref. [52:] and proved in a 
different manner. 


VII. HANKEL MATRIX DECOMPOSITIONS OF THE P-MATRICES OF A n , B n , D n 

The P-matrices obtained from the generating formulas El, ESI and (BOD , for the groups of type A n , B n and D n , 
respectively, have an interesting Hankel matrix decomposition, that will be pointed out in this section. 

Let’s recall the definition of a Hankel matrix. A Hankel matrix is a matrix that has equal elements along a same 
upgoing diagonal. For a Hankel matrix H one then has the following equality among matrix elements: 

Hi,j = Hi—i,j+x 

In other words, the elements Hij of a Hankel matrix H depend only on the sum i + j, and the matrix elements with 
equal sum of the indices are equal. A square Hankel matrix of order n is completely determined by the 2n — 1 elements 
located in the first column and in the last row of H. We will need only Hankel matrices with zero elements above the 
auxiliary diagonal (the auxiliary diagonal of a square matrix H is the diagonal going from the element H n i to the 
element Hi n ). Let’s use the symbol H(y) to define a square Hankel matrix that has its last row equal to the vector 
v, and is zero above the auxiliary diagonal. If n is the dimension of v , then n is the order of H(y). 

With the groups B n and D n we will need the matrices H n (p ) and Ha l \p)- H n (p) is defined as follows: 


H n (p) = H(v(p )), 


v{p) = (v 1 (p),v 2 {p), ■ ■ -,v n (p )), 


v a (p) = {n- a+ 1) Pa—i, 


Va = 1,..., n , (42) 
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in which one has to consider po = 1, and Hi n \p), with a < n, is the n x n block Hankel matrix that has the Hankel 
matrix H a (p) in its upper left corner, and is zero elsewhere. For example, 


h 5 (p) = 


/ 

0 

0 

0 

0 

5 > 


/ 

0 

0 

3 

0 

°\ 

0 

0 

0 

5 

4pi 

, H { ^\p) = 

0 

3 

2pi 

0 

0 


0 

0 

5 

4pi 

3p 2 


3 

2pi 

P2 

0 

0 

V 

0 

5 

4pi 

3^2 

2p 3 


l 

0 

0 

0 

0 

0 

5 

4pi 

3^2 

2p 3 

Pi / 


0 

0 

0 

0 

0/ 


For the groups B n , one proves the following theorem. 

Theorem 10 If the basic invariant polynomials of B n are those given in Eq. the corresponding P-matrix, that 

can be determined using Eq. w, can be written in the following way: 


P(p) = Aj2PaH^( P ) 


(43) 


a= 1 


Proof. In Section IX El 


For the groups of type D n there is a similar theorem expressing their P-matrices in terms of Hankel matrices, but 
some definitions have to be first given. The P-matrix P(p), given in Eq. (1401) . corresponding to the basic invariant 
polynomials of D n given in Eq. distinguishes the elements in the last row and column of P(p) from the others. 

For this reason, it is convenient to define the symbol I( n ' ) for the n x n diagonal matrix with the unit matrix I n -i , 
of order n — 1, in its upper left corner and is zero in all its last row and column, and the symbol I^' 1 for the matrix 
J< n) = I n — /("), where I n is the unit matrix of order n. The diagonal of the matrix I has n — 1 l’s and a 0 in its 
last element, vice versa, the diagonal of the matrix 1^' 1 has n — 1 0’s and a 1 in its last element. For example: 


/( 5 ) = 


/ 

1 

0 

0 

0 

0 

\ 


/ 

0 

0 

0 

0 

0 

\ 

0 

1 

0 

0 

0 


r(5)_ 

’ i 0 — 

0 

0 

0 

0 

0 



0 

0 

1 

0 

0 



0 

0 

0 

0 

0 



0 

0 

0 

1 

0 



0 

0 

0 

0 

0 


V 

0 

0 

0 

0 

0 

) 


l 

0 

0 

0 

0 

1 

/ 


The matrices I^ and I^ 1 ' 1 can be used to define the matrices H[ n \p) and H^ n \p) in the following way: 


jr< n) (p)=J<"> H n (p) /<">. 


H, 


(n) 


(p) = H n (p) - H^ n \p) - p n I 0 


(«) 


The matrix H^ n \p) has the same elements of the matrix ff„(p) but with zero’s in the last row and column, and 
Hq n \p) contains the non-diagonal elements of the last row and column of H n (jp) and is zero elsewhere. For example: 


H i t\p) = 


/° 

0 

0 

0 

°V 


/ 

0 

0 

0 

0 

5 \ 

0 

0 

0 

5 

0 

, H^\p) = 

0 

0 

0 

0 

4 pi 

0 

0 

5 

4 pi 

0 


0 

0 

0 

0 

3P2 

0 

5 

4 Pi 

3P2 

0 


l 

0 

0 

0 

0 

2p 3 

Vo 

0 

0 

0 

o / 


5 

4pi 

3P2 

2p 3 

0 / 


For the groups D n , we can now give the following theorem. 


Theorem 11 If the basic invariant polynomials of D n are those given in Eq. \39\) . the corresponding P-matrix, that 
can be calculated from Eq. w> can be written in the following way: 

n —1 

P(p) =4 J2p a Hi n \p) + 4 pi H[ n \p) + 2 p n H^ n) +p n _ 1 4 n) . 

a= 1 


( 44 ) 
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Proof. In Section IX FI o 

The P-matrix of A n , corresponding to the basic invariant polynomials defined by Eqs. (EH) m, that can be 
constructed using Eq. (15^1) . can be written in terms of block Hankel matrices, too, but it is first convenient to give 
the following definitions: 


Y n =H(e i), ei = (1,0,... ,0), (n - 1 zero’s) 

K n (p) = H{u{p )), u(p) = (ui(p),...,u n (p)), U a (p) = V a -i(p) = (n - a + 2) p a -2 , Va = 1,..., n, 

where one has, by definition, po = p _i = 0, (so u\(p) = u 2 (p) = 0). Note that K n {p) can be obtained from H n {p ) 
by considering po = 0 instead of po = 1 (and this means a matrix equal to H n (p ) but with 0’s along its auxiliary 
diagonal), and by shifting all columns of the resulting matrix one position to the right or, by shifting all rows of the 
resulting matrix one position towards the bottom, if one prefers). 

Both Y n and K n (p ) are n x n Hankel matrices, with Y n that has l’s along the auxiliary diagonal and zero elsewhere. 

Let’s define Y'i"' 1 and K < a "' > (p). with a < n, be the nx n block Hankel matrices that have the Hankel matrices Y a and 
K a (p) in their upper left corner, and are zero elsewhere. 

For example, 


K${p) = 


Let’s also define the nx n symmetric (non-Hankel) matrix R^ n \p) in the following way: 

R< ab(p) = t( n +!) min(a,6) - a&]p Q _ip b _i, a,b=l,...,n, (45) 

in which one has to consider po = 0. For example, 


/° 

0 

0 

0 

0 

\ 


/° 

0 

0 

0 

0 

\ 


/° 

0 

0 

1 

0 

\ 

0 

0 

0 

0 

0 


, Kf\p) = 

0 

0 

0 

0 

0 


V ( 5 ) - 
5 x 4 — 

0 

0 

1 

0 

0 


0 

0 

0 

0 

4pi 


0 

0 

0 

3pi 

0 


0 

1 

0 

0 

0 


0 

0 

0 

4pi 

3p 2 



0 

0 

3pi 

2p2 

0 

) 


1 

0 

0 

0 

0 

) 

\0 

0 

4pi 

3p2 

2ps 


Vo 

0 

0 

0 

0 


Vo 

0 

0 

0 

0 


R {5) (p) = 


/0 0 0 0 0 \ 

0 8 p\ &P1P2 4pip 3 2pip 4 

0 6pip 2 9 pi 6p 2 P3 %P2P4 

0 4pip 3 6p 2 P3 8 pi 4:p 3 p 4 

\ 0 2pip4 3p 2 P4 4p 3 p 4 5 pi ) 


For the groups A n one proves the following theorem. 


Theorem 12 If the basic invariant polynomials of A n are those defined by Eqs. \28\) \31\) . the corresponding P- 
matrix, that can be calculated using Eq. t32 1) . can be written in the following way: 

n n 

P(p) = RW (p) + 2Y / (a + l) Pa F a (n) - (n + 1) ^ p a K W (p) (46) 

a= 1 a—1 

Proof. In Section IX G I o 

Regarding the groups of type S n , I could not find any basis of invariant polynomials such that the corresponding 
P-matrices have a simple decomposition in terms of Hankel matrices. The problem arises with the variable pi , of 
weight 1, that cannot be written in the P-matrix of S n in a simple “symmetric” way. Obviously, there exist a basis of 
invariant polynomials of S n in which the P-matrix (pi excluded) has the same simple form of the P-matrix of H ra _i, 
that can be written as specified by Theorem Q2] (but with the labels of the basic invariant polynomials of A n _i that 
have to be increased by 1, and with one more column and row in front of the n—l others), but it contains the variable 
Pi in a complicated way. 
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VIII. EXCEPTIONAL GROUPS AND NON-CRYSTALLOGRAPHIC GROUPS 

The main results presented in this article are the generating formulas given in Theorems [6j [7] [ 8 ] and [9] These are 
useful (and sometimes necessary) to determine the P-matrices and the A-vectors in the case of the finite reflection 
groups in the infinite series S n , A ni B n , and D n , respectively, for which there exist a group for each value of the 
rank n. On the contrary, the exceptional groups Eq, E 7 , E$, P 4 , G -2 and the non-crystallographic groups H 3 , are 
defined only for specific values of the rank and a generating formula for them is not necessary. Their P-matrices can 
be found in some published articles, for example in Refs. [39j,[52| (groups Eq , P 4 , G 2 , P 3 , P 4 ), [M] (groups E e , P 4 ), 
M (groups E 7 , E$). 

There is another infinite series of finite reflection groups, that one of the dihedral groups I 2 {m), m > 2, the groups 
of symmetries of the regular m-sided polygon (that includes the crystallographic groups A\ x A± ~ P(2), A 2 — P(3), 
B2 — ^( 4 ), and G2 — ^2(6)). The roots, the basic invariant polynomials, the P-matrices and the A-vectors of the 
groups of type I 2 (m) are well known, see for instance Refs. [13, lH, H3, 0, IHj. For completeness, I report briefly, 
for the groups of type him), m > 2 , one of the possible choices for the simple roots, the roots, the basis of invariant 
polynomials, the corresponding P-matrix, the irreducible active factors and their A-vectors. 

• Simple roots. The 2 simple roots can be defined as the vector a\ = (0,1) and the vector 0:2 obtained from oq 
with a rotation about the origin (0, 0) by an angle <p m = — (n — —). The two simple roots obtained in this way 
have the same length. This is necessary if m is odd, but not necessary if m is even, because if m is even two 
different root lengths are possible. To be crystallographic, the groups I 2 {A) — P 2 and / 2 (b) — G2 require the 
lengths of the two simple roots in a ratio equal to s/2 and \/3, respectively. 


• Roots. The 2m roots are obtained by rotating a\ and 0.2 by the angles =^, V* = 0,1,..., m — 1. The positive 
roots are the roots with a positive number as their first non-zero coordinate. 


• Basic invariant polynomials. 


Pi(x) = x\ +xl , 


p 2 {x) = Re[(zi +ix 2 ) m ] = ^ [{xi +ix 2 ) m + (xi -ix 2 ) m ] , i = \/-f. 

The polynomials pi(x) and p 2 (x) are invariant for the reflections corresponding to the two simple roots a± and 
02 , and then for all the group I 2 (m) generated by them. 




P-matrix. The P-matrix corresponding to the basic invariant polynomials p\[x) and p 2 (x) defined above is the 
following one: 


Pip) 


4 pi 2mp 2 \ 
2 mp 2 m 2 p™- 1 ) 


One has det (P(p)) = 4m 2 (p™ — p 2 )- If m is even, det(P(p)) = a+(p)a_(p), with a±(p) = p± P 2 two 
distinct irreducible active polynomials. 


• A-vectors. 


Ai det(P)) (p) = (4m, 0). 


If m is even one also has: 


Ai a±) (p) = (2 m, ±toVi" /2 V 
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The basis of invariant polynomials for the group / 2 (to) just described is at the same time a flat basis, a det(P)-basis 
and a canonical basis. If to is even, the two irreducible active factors of det(P(p)) define two algebraic curves a±(p) = 0 
that contain the boundary of the set S representing the orbit space of 4 (to). The algebraic curve a+(p) = 0 contains 
the images through the orbit map of the to /2 reflecting lines orthogonal to the m roots obtained by rotating the simple 
root by the angles —, Vi = 0,1,..., m — 1. The algebraic curve a_(p) = 0 has the same meaning, but is related 
to the simple root c^- This can be seen by picking a point in a reflecting line and by seeing which algebraic curve 
contains its image. For example the point (aq, aq) = (aq, 0) lies in the reflecting line orthogonal to on, and its image 
through the orbit map is the point (pi, P 2 ) = (xf, x™), that satisfies the equation a_(p) = 0, but not a+(p ) = 0. Vice 
versa, the point (aq, X 2 ) = (aq cos aq sin //-) lies in the reflecting line orthogonal to a 2 , and its image through the 
orbit map is the point (pi,P 2 ) = (x 2 , ~ x ™)i that satisfies the equation a+(p) = 0, but not a_(p) = 0. If to is even 
there are two different conjugacy classes in 4 (to) and two different singular orbit type strata of dimension 1 that are 
represented in S by two primary strata lying in the two curves, respectively. If to is odd there is a single conjugacy 
class in 4(wi) and a single singular orbit type stratum of dimension 1 . The image in S of this orbit type stratum is 
formed by two different primary, lying in the two curves, respectively, that are separated by the origin of M 2 . 


IX. EXAMPLES 


Example 1. 

This example shows the construction of the rotation matrix R 4 , with Eq. (EH), of the basic invariant polynomials of 
A 3 , with Eqs. Ell) ED- and of the P -matrix of A 3 , with Eq. (1321) . 

The matrix R n + 1 , defined in Eq. (1271) . and its transpose, in the case n + 1 = 4, that is n = 3, are the following: 


/ 1/42 

-1/42 

0 

0 ^ 


( 1/42 

1/46 

1/412 

1/2 \ 

1/46 

1/46 

-2/46 

0 

, Rj = 

-1/42 

1/46 

1/412 

1/2 

1/412 

1/412 

1/412 

-3/412 

0 

-2/46 

1/412 

1/2 

V 1/2 

1/2 

1/2 

1/2 ) 


V 0 

0 

-3/V12 

1 / 2 / 


One may easily verify that R 4 R/I = RjR 4 = I 4 , that R 4 transforms the column unit vector 114 = (^, 4 1) into the 

column unit vector a = (0,0,0,1), and that Rj makes the opposite transformation. 

To obtain the basic polynomial invariants s' a (x'), a = 1,..., 4, in Eq. (1291) . one has to make the variable substitutions 
x —7 Rjx 1 in the basic polynomial invariants s a {x), a = 1,... ,4, given by Eq. (l28l) . with n = 3. One finds so the 
following basic invariant polynomials of £ 4 : 

sj(z') 

4 (A) 

4 (*') 

44 ') 


The basic invariant polynomials of A3 are obtained from the basic invariant polynomials s' a (x'), a = 1 ,... , 4 , of £4 
in the way specified in Eqs. EH and ED- One so finds: 

Pi{x ) = aq +x 2 + x 3 , 

P2(x') = — ( 3 V 2 x[ 2 x' 2 — V2x' 2 3 + 3x^ 2 x' 3 + 3x' 2 2 x' 3 — 2x' 3 3 ^ , 

P 3 ,{x') = - x ' 3 (l 2\/2 x [ 2 x 2 — W 2 x ' 2 — 6 x[ 2 x ' 3 — 6 x 2 2 x ' 3 + x' 3 3 ^j. 


-2 , 1,1, 

- X2 + — X3 + - Xi 


1 , 1 , 1 , 1 ,\ /-1 , 1 , 1 , 1 ,\ /-: 

vT 1 + 7T 2 + Til" 3 + 2 X V + V7T 1 + 7T 2 + w Xs + 2 X V + b 

: 4 

+ — 7=4 -4 7= x 3 + —£ 4 ') ( — 7=4 + —7=4 4 7=4 + —* 4 ') 4- ... 4- ( —j=x' 2 4 7=4 4“ -4 

y/2 VE 412 d 2 V VV2 VE \/l2 J 2 4 ) WE 2 x/12 J 2 

1 ( , 2 . , 2 . , 2 \ , 3 ,2 

- + x 2 + x 3 J + -x 4 

- ( 3VEx\ 2 x' 2 — V6.4 3 4- 3V3i( 2 x( + 3V3x 2 2 x 3 — 2y/3x 3 3 — 9a;( 2 x\ — 9x' 2 2 x' 4 — 9x 3 2 x' 4 + 94 3 ) 


- (-36v / 2 4 2 44 + i2v / 2 4 3 4 + 18 4 2 4 2 1 1Q ~ ,2 ~' 2 

/3x'n 2 X 3 x' A — 8\/3 4 3 4 - 


vh x ' 3+1 2 x ' 4] = 


-3 , 1 ,\ 

—=X 3 + -Xa = 

x/l 2 d 2 / 


- 18 — Sx'^ + 12\/6 — Ay/ftx^x'^ + ^Vdtx'-^x^'. 


.0 .0 


( 47 ) 
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The corresponding P-matrix and A-vector are obtained with Eqs. (l32l) and (1551) and are the following: 

A ( d et(P))( p) = (24 ; 0j _ gpi ). 


4 Pi 

6 p 2 

8 p 3 \ 

6 P2 

4 Pi + 8 p 3 

2 piP2 

8P3 

2 piP2 

3p| - 8 pip 3 ) 


(48) 


Example 2. 

This example shows the construction of a flat basis of invariant polynomials (recall the discussion with Eq. (12U1) '). 
together with its P-matrix P(p) and A-vector of det(P(p)). Let’s consider the case of A 3 , so we can use the basic 
invariant polynomials, the P-matrix and the A-vector obtained in Example 1. 

In a flat basis, the condition expressed by Eq. dH is true. Then one has to find the basis transformation that 
transforms the P-matrix in Eq. (1-151) to a P-matrix satisfying Eq. (15U1) . In the present case one has to eliminate the 
variable p 3 from the element P 33 (p), only. One then writes the most general allowed basis transformation p' = p' ( p ) 
that transforms the variable p 3 , only, that preserves homogeneity and such that J 33 (p) = 1. One then writes: 


Pi(p)=Pi> P2(P)=P2, P 3 (p) = P3 +cipj, 

with Ci a real constant. This transformation implies the following jacobian matrix: 

m = 

that gives, using Eq. m, the following expression: 



£33 (P ; ) 


&=& ( v ) 


= J 3 afr) Pab(p) Jb3(P ) = 4 ClPl Pllip) + 4 Cl Pi Pl 3 (p) + ^33 (p) = 


(49) 


(50) 


= 4ciPi4pi +4 ci pi 8 p 3 + (3pi - 8 pip 3 ) = 16 cfp? + 3 pi + 8 (4ci - l)pip 3 


In the above calculation, it is not necessary to write the P-matrix with the new variables p', using the inverse trans¬ 
formation p(p') and Eq. (1171) . because the condition (l20l) can be obtained either by deriving the matrix P(p ')\ p '—p ( p ) 
with respect to p 3 or by deriving the matrix P(p') with respect to p 3 . Then condition (l20l) is fulfilled if 4ci — 1 = 0, 
that is, if ci = |, and Eq. (l49l) gives the following transformation from the basis (1771) to a flat basis: 


Pi(p) =Pi 

With the transformation found, Eqs. m and m give: 
P(p') = 


Pi(p)=P2, P^(p) =P3 + ^P?- 


4 Pi 

6 p ' 2 

8 P 3 

6 pi 

2 Pi 2 + 8 p^ 

5 P 1 P 2 

8 P 3 

5 P 1 P 2 

Pi 3 + 3p! 


A (det(P))(p/) = (24, 0, 4pi). 


When the rank n is large it may be somewhat difficult to find the transformation to a flat basis and a systematic 
method to follow is necessary. First of all, in place of the P-matrix P(p) in Eqs. (1171) or (Hhd, one can use the part 
of it (let’s call it Q{p)) containing (linearly) the highest weight variable p n , only. A quick look at the Hankel matrix 
decomposition in Eqs. (1751) . (1771) or (1751) . is sufficient to see what this matrix Q{p), containing p„, is. For example, in 
the case of the groups of type A n this matrix is Q{p) = 2 (n + 1) p n — (n + 1) p n K^ (p). We then use Eq. (1151) 

= J(p) Q(p) J T (p), and take the derivative of the resulting matrix with 


with the matrix Q(p), that is: Q{p') 


p'=p'(p) 


respect to p n . Let’s define A(p) = yJ{p) Q{p ) J T (p)J = J(p) d g^ ^ 7 T (p). If p' is a flat basis, A(p) is a constant 
matrix, and this means that all its elements of positive weight are zero. By requiring this, one finds the basis trans¬ 
formation. When n is large, it is convenient to proceed in small steps. Let i be the smallest index such that A-^-i (p) is 
of positive weight. Write the most general allowed basis transformation p' a = p' a {jp ), Va = 2with = 1, 

and determine all the free coefficients of this transformation formula necessary to obtain A a ^(p) = 0, Va = 1,..., i. 
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Then, for each b = i + 1,..., n, write also the most general allowed basis transformation p' b = p' b (p), with = 1 

and determine all the free coefficients (in p' b {p) as well as in p' a (p), a = 2 ,..., b — 1 , that are still undetermined) 
necessary to have A a ^(p) = 0, Va = 1,..., b. In this way one arrives at a flat basis. It is convenient to program 
a computer to do the transformation of Eqs. o or m- When applying this method with the groups of type 
D n , one has to remember that the highest weight variable is p n -\ and not p n , in particular in Eq. (1441) and in Eq. (1201) . 

Example 3. 

This example shows the construction of a det(P)-basis of invariant polynomials (recall the discussion with Eq. (1211) 1. 
together with its P-matrix P{p) and A-vector ,\( det C p )) (p). Let’s consider the most general allowed basis transformation 
p' = p'{p), with J aa {p ) = 1, Va = 1 ,n. Eq. (fl9l) states how a A-vector transforms. To obtain a det(P))-basis, 
one has to require that all components, except the first one, of the transformed A-vector of \( det ( p ^ (p) be equal to 
zero. Using the formula: 


Ai det(P) V) ^ Jab ^ A i det(P)) (p) - Va = 1 ,...,n, 

P=P (P) b 


equivalent to that one in Eq. m, one has to require: 


(51) 


Ai det(P) V) 


p'=p'(p) 


= 0 . 


Va = 2,..., n. 


(52) 


Let’s consider the case of A 3 , so we can use the basic invariant polynomials, the P-matrix and the A-vector obtained 
in Example 1. The most general allowed basis transformation with the required properties is written in Eq. (1491) . and 
its jacobian matrix is written in Eq. (USD- The P-matrix P{p) and the A-vector \( det ( p )) (p) ; corresponding to the 
basic invariant polynomials in Eq. (El), are written in Eq. (1481) . The transformed A-vector, obtained with Eq. m, 
is then the following one: 


A (det(P )) ( ^) 


= (24, 0, 8(6 c\ — l)pi), 

p'-p (p) 


and one obtains a det(P)-basis if c\ = A. Using C\ = | in Eq. (TTO . one finds the following transformation from the 
basis El) to a det(P)-basis: 


Pi(p)=Pi, p'i{p)=Pi, p' 3 {p) =P3 + ^p\- 
With the transformation found, Eqs. ED and m give: 


/ 4 pi 

Qp '2 

OO 

CO " 


6 P 2 

Ip'i 2 + &P3 

4 P 1 P 2 

, A (det(p) )(p') = (24, 0, 0) 

{&P3 

4 P 1 P 2 

Ip'i 3 -Ip'iP3 + ZP2 2 ) 



If the rank n is large it may be convenient to transform one by one the variables from that one of the greatest weight 
d n to that one of the lowest weight greater than 2. When writing the most general transformation, one may ignore 
the terms not depending on pi, because they do not transform a non-a-basis into an a-basis or vice versa. When 
one is transforming the variable indexed with a, one has to find the free coefficients of this transformation in such a 

way for which Ai dct ^ P '*' ) (p) = 0, only. What is here said for a det(P)-basis can be performed for a general a-basis, as well. 
Example 4. 

This example shows the construction of a canonical basis of invariant polynomials, qi(x),... ,q n (x), (recall the dis¬ 
cussion with Eqs. (l22l) and (l23ll ). together with its P-matrix P(q) and A-vector A( det ( p ))(g) of det(P(g)). Let’s write 
P\{x), ... ,p n {x) for the basic invariant polynomials of the groups of type A r , B n , D n , given in Sections IIVIIVI1 with 
the same labeling conventions, and pi 1 (x),... ,pi n (x) for the same basic invariant polynomials in which the order 
is according to their degrees: d,. < d, . +1 , Vj = 1,... ,n — 1. (Except for the groups of type D n , even n, one has 
ij = j, and the inequality is always strict). One has then qi(x) = pi(x) = ||x|| 2 . Suppose that the k basic invariant 
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polynomials of lowest degree of a canonical basis are already known, with k = 1 , 2 ,..., n — 1 , that is, suppose that 
(^(x),... ,qi k {x ) are already known, and that one wishes to find the next canonical invariant q ik+ 1 (x). One has then 
to write the most general homogeneous polynomial of degree di k+1 , r(x), written in terms of pi 1 (x), ... ,Pi k+1 {%), such 
that the coefficient of Pi k+1 ( x ) is 1 (except in the case of D n , even n, that when di k+1 = n and Pi k+1 (x) is the second 
basic invariant of degree n, one has to keep the coefficient of Pi k+1 (x ) arbitrary). By requiring the conditions expressed 
by Eq. (j22]l . that is ( qi a ,r } = 0,Va=l,...,fc, one finds the arbitrary coefficients that render r(x) canonical. Then 
one defines qk+i(x) = r(x) and proceed in the same way with the next value of k. In the case of D n , even n, one has 
the two more conditions expressed by Eq. (1231) . 

Let’s consider the case of A 3 , so we can use the basic invariant polynomials, the P-matrix and the A-vector of 
defyPfy)) obtained in Example 1. The most general allowed basis transformation is expressed by Eq. (l49l) . For 
k = 1 one has r(x) = p 2 (x), without arbitrary coefficients, so one defines q 2 (x) = p 2 (a?). For k = 2 one has 
r(x) = P 3 (x) + Cipi(x) 2 , with ci an arbitrary coefficient. The condition {qi,r) = 0 gives 0 = {qi,P 3 + C\p\) = 
(< 71 , 753 ) + c 1 (qi,Pi) = 2(— 1 + 10ci) qi(x). The condition (q 2 ,r) = 0 is identically satisfied and gives no conditions. 
One then has a canonical basis if c-[ = Using ci = ^ in Eq. (l49l) . one finds the following transformation from the 
basis (S3 to a canonical basis: 


<li(p)=Pi, q 2 (p)=P2-, 

With the transformation found, Eqs. (El) and ED give: 


P(q) = 


q 3 (p) =P3 + -^Pi- 


/ 4<7i 

6(72 

8q 3 

6 C/2 

fqi 2 + Sq 3 

f <7192 

\ 8 <7 3 


M <7i 3 - f <7 i<73 + 3<7 2 


A (det(P)) (g)= ( 24 , 0,~qi 


This method can be followed also if the rank n is large, but the calculations may become very long. 


X. PROOFS 

A. Proof of Theorem [6] for S„ 

Proof of the generating formula (1251) . 

The matrix obtained from the generating formula (1251) is symmetric, so we can simplify the formula and the proof by 
supposing a < b. If we prove that Eq. (l25l) is true for a < b then it is true also for b < a, hence for all allowed values 
of a, b. A further simplification is possible if we define, in addition to Po(x) = 1, that p a (x) = 0, Va < 0 and Vo > n. 
We will then prove the validity of the following formula for S n> equivalent to the generating formula (l?5l) : 

a 

Pabip) = Pbaip) = (n + 1 - a) p a -i Pb-i - (a + b - 2 i) pi-\ p a+6 _i_i, Va<6 = l, ...,n, (53) 

i=a-\-b—n— 1 

in which one has to consider po = 1 and p a = 0, V a < 0 and V a > n. 

The proof can be obtained by induction on the rank n. The lowest possible value for n is n = 1, and with 
n = 1 Eq. (1531) is trivially verified: the only basic invariant polynomial is pi(x) = X\, with which Eq. © gives 
Pn(x) = 1, and Eq. (15^1) gives Pu{p) = 1. Let’s verify Eq. (HkH) also for the less trivial value n = 2. From 
Eq. (l24l) . for n = 2, the basic invariant polynomials are p\{x) = X\ + x 2 and p 2 (x) = X\X 2 , their gradients are 
then Vpi(x) = (1,1) and Vp 2 (x) = (x 2 ,xi), and Eq. © gives: Pn(x) = 2, Pi 2 (x) = P 21 {x) = aq + x 2 = Pi(x), 
P 22 (x ) = x\ + x\ = pfyx) 2 — 2 p 2 (x). Using Eq. (1551) one finds the same form of the corresponding P-matrix: 
Pi i (p) = 2, P 12 (p) = P 2 i(p) = Pi, P 22 (p) = p\ — 2 p 2 . We have then verified that Eq. (l53l) correctly gives the matrix 
elements of the P-matrix P(p) for the basic invariant polynomials E51) of S n , for n = 1, 2. 

Using the induction principle, we suppose, for a general n > 2, that Eq. (1551) is true for S n , and prove in all 
generality that Eq. (l53l) is true also for n substituted by n + 1, that is for S„+i. This implies then that Eq. (1531) is 
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true for S n for all n > 2 . 

We use here the symbols x and y for generic vectors of R" and R n+1 , that is: x = (aq, X 2 , ■ ■ ■, x n ), and 
y = [x\, X 2 , ■ ■ ■, x n , x n +i). S n acts in R™ and its basic invariant polynomials are given by Eq. (EhO) : p Q (x), 
a = 1,..., n. S n +1 acts in R ” +1 and its basic invariant polynomials are also given by Eq. (El, but with n substituted 
by n + 1 , and we write them as p a {y), a = 1 ,..., n + 1 . 

The basic invariant polynomials of SV 1+1 can be expressed in terms of those of S n in the following way: 

Pi{y ) = Pi{x) + x n+1 , 

P 2 (y) = P 2 {x) +pi{x)x n+1 , 

P 3 {y) = P 3 {x)+p 2 {x)x n +l, 

Pn(y) = Pn{x) + p n -l{x) X n+ i , 

Pn+l{y) = Pn{x)x n+ 1 . 

Using the definitions po(x) = 1 and p n +i(x) = 0, one has the following general formula: 

Pa{y) =Pa{x) + p a -i(x) X n+ 1 V a = 1, 2, . . . , 71 + 1. (54) 

Let’s define the following “nabla” operators that map a scalar function into a vector of R" +1 : 

7 - (JL JL 9 

y V 9x ± ’ dx 2 ’ ’ " ’ dx n 1 dx n+ i 

One has then, for Eqs. (1M1) and (l55l) . Va = 1,..., n + 1: 

Vy Pa (y) = V y [p a (x) +Pa-l(x) X n +\] = V xPa(x ) + X n +1 V xPa-l(x ) + p a -l{x) e n +l , (56) 

where e n +i is the canonical unit vector of R ” +1 (with 1 in the last, (n + l)-th, coordinate and 0 elsewhere). 

We can now calculate a general element of the matrix P(y) from Eqs. EJ and (15611 . using the basic invariant 
polynomials p a (y), a = 1,..., n + 1, of S'n+i. For all a, b = 1,..., n + 1, we have: 

Pab{y ) = V.yp Q (y) ■ Vyp b (y) = 

= [v x p a (x) + x n+ i V x p a -]_(x) + p a -i{x) e n+ i] ■ [S7 x p b (x) + x n+1 v x p b -i{x) + Pb-i(x) e n+ i] = 



V x = 


d d 


d 


dxi ’ dx 2 ’ dx r 


,0 


— VxPo (xj • V xPbi.x') T X n -^-i ^xPa(x) • V ' x Pb—l{x ) T 0 T X n -\~i \7 x Pa— 1 (*^) * ^xPb(x) ~\~ 


+ X 2 n+1 X7 x p a -l(x) ■ \7 x p b -l{x) + 0 + 0 + 0 + Pa-l(x) Pb-l(x) = 


= Pab(x) + X n+ 1 [Pa,b- l(z) + Pa- 1 ,b{x)\ + xl +1 P a - l,b-l(x) + p a -l (x) p b -l (x) . (57) 

Of course, in the preceding expression, P c d(x ) = 0, if one of the indices c, d, (that in the above expression can be equal 
to a, b, a — 1 , b — 1 ) is equal to 0 or to n + 1 , because for S n the allowed indices run from 1 to n, and the gradients of 
Po(x) = 1, and p n+ ±(x) = 0 are null vectors. By the induction hypothesis, Eqs. (l53ll and (l25l) are true for S n , so we 
can use them in the preceding equation. We obtain then, Va<6=l,...,n+1: 

a 

P ab {y) = {n + 1 - a) Pa-i(x) Pb-i(x) - ^2 (a + b- 2i) p l - 1 {x) p a+b -i-i(x) + 

i=a-\-b—n— 1 


+ ^n+l 



— min(a, b 


1)) Pa —1 {x) Pb- 2 {x) 


min(a,6— 1) 

F, (a + b- 1 - 2i) Pi-i(x) p a+b - 2 -i(x) + 

i=a-\-b—n —2 
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a— 1 


+ (n + 2 - a) p a - 2 {x) pb-i{x) - y (a + b - 1 - 2z) Pi-i(x) p a+b - 2 -i{x) > + 


i=a-\-b—n —2 


a—1 


+ :r ^ +1 < (n + 2 — a) Pa-20) Pb- 2 (x) - V (a +&- 2 - 2*)pi_i(a;) p a+ f,_ 3 _i(a;) > + 


i=a+b—n—3 


+ Pa-lO) Pb-l(x) . (58) 

The first four lines of the above formula correspond to the terms of Eq. m containing the matrix elements P ab (x), 
Pa,b-i(x), P a -i t b(x) and P 0 _i i f,_i(x), respectively. In the second line instead of Eq. (15^1) we used the original formula 
(1251) and could not simplify further min(a, b — 1), because in the present case a < b (our hypothesis) we have to 
distinguish the case a < b from the case a = b. We shall do that later. The limitations on the indices regarding these 
four lines are as follows: 1st line: a < b < n, 2nd line: a < b and a < n and 2 < 6 , 3rd line: 2 < a < b < n, 4th 
line: 2 < a and 2 < b and a < b, respectively. These limitations are automatically satisfied, that is, if b = n + 1, the 
first and third line vanish, ifa = 6 = n + l the second line vanishes, if a = 1 the third and fourth line vanish, and if 
a = b = 1 the second, in addition to the third and fourth line vanish, using p-\ = 0. We can then avoid to take into 
account in the following the limitations on the indices. We can then write, Va< 6 = 1,..., n + 1: 

a 

Pab(y) = {n + 2 - a) p a -i(x) Pb-i(x) - ^2 (a + b - 2i) Pi-^x) p a+ b-i-i{x) + 

i=a-\-b—n— 1 


+ ^n+l 



— min(a, b 


1 )) p a -i(x) Pb- 2 {x) 


min(a,b— 1) 

F, (a + b- 1 - 2i) Pi-i(x) Pa+b- 2 -ii.x) + 
i=a-\-b—n —2 


+ (n + 2 - a) p a - 2 (x) p b -i(x) - 


a —1 

E 

i=a-\-b—n —2 


(a + b-l-2i) Pi-i(x) p a+b - 2 -i(x) > + 


+ X n+l S (, n + 2 - «) Pa- 2 (x) Pb- 2 {x) - ^ ( a + b - 2i') p z , _ 2 (x) p a+b _ 2 -i' ( x ) 

v i'—a-\-b—n—2 

The term inside the first curly brackets, that one multiplying x n +i, can be further simplified. 

One notes that both sums that appear inside the curly brackets contain equal terms up to the value i = a — 1, because 
a < b. Moreover the first sum contains one more term if a — 1 7 ^ min(a, b — 1), that is if a < b (hence not if a = b), 
that corresponds to i = a. One has then: 

min(a,6— 1) a—1 

y, (a + b- 1 — 2i) Pi-i(x) p a+b - 2 -i(x) - y (a + b- 1 — 2i) Pi-i{x) p a+b -2-i(x) = 

i=a-\-b—n —2 i=a-\-b—n —2 



a—1 

= - y 2 (a + b - 1 - 2i) Pi-i(x) p a+b - 2 -i(x) - (1 - S ab )(a + b - 1 - 2a) p a -i(x) p a+b - 2 -a(x) = 

i=a-\-b—n —2 

a— 1 

= (a - b + 1 - Sab) Pa-i(x) Pb- 2 {x) - y 2 (a + b- 1 - 2i) Pi-i(x) p a+b - 2 -i{x ), 

i=a-\-b—n —2 


where S ab is the Kronecker delta. 

The two terms without the sum inside the curly brackets multiplying can be transformed in the following way, 
that can be easily verified by distinguishing the different cases a < b and a = b: 

(n + 1 - min(o, b - 1 )) p a -i(x) p b - 2 (x) + (n + 2 - a) p a - 2 (x) p b -i(x) = 
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= (n + 1 - a + Sab ) Pa-i(x) p b - 2 (x) + (n + 2 - a) p a - 2 {x) p b -i{x) . 

One can then write the term inside the curly brackets multiplying as follows: 

x n +i {■•■}= x n+ i j(n + 1 - a + 5 ab ) p a -\{x) p b - 2 {x) + (n + 2 - a) p a - 2 (x) p b -i(x) + 

a-1 ^ 

+ (a - b + 1 - Sab) Pa-i(x) p b - 2 (x) - Y, 2 (a + b - 1 - 2i) pt-i(x) p a +b- 2 -i(x) > 

i=a-\-b—n—2 ) 

= x n+ i l (n + 2 - b) p a -i(x) p b - 2 (x ) + (n + 2 - a) p a - 2 {x) p b -i(x) + 


a— 1 


Y 2 {a + b - 1 - 2i) Pi-^x) p a+b - 2 -i(x) > . 


i=a-\-b—n —2 / 

We have then, finally, the following expression, Va, 5=l,...,n+l: 

a 

Pab(y) = {n + 2- a) p a -i(x) p b -i(x) - Y {a + b - 2i) Pi_ x {x) p a + b -i-i(x) + 

i=a-\-b—n— 1 

+ x n+ i < (n + 2 — 6) Pa-i{x) p b - 2 (x) + (n+2-a) p a - 2 (x ) p b _i(a;)- 2 (a + 6-l-2f) Pi_i(x) p a+b _ 2 _i(x) ^ + 

V i=a-\-b—n— 2 J 


+ ^n+i S (n + 2-a)p a _ 2 (x)p b _ 2 (;r) - Y (a + 5 - 2i) Pi- 2 {x) p a+b - 2 -i(x) 

v i=a-\-b—n —2 > 


(59) 


Using Eq. m in the case of S n +i, that is for n substituted by n + 1, and Eq. (i54l) . one would instead find, 
Va<6=l,...,n+1: 

a 

Pab{p(y)) = (n + 2 - a) p a -i(y) Pb-iiy) - Y (a + b - 2i) p l - 1 (y) p a +b-i-i{y) = 

i=a-\-b—(n-\- 1) —1 

= (n + 2 - a) (Pa-l(x) + Pa- 2 (x) Xn+l) (Pb-l(x) + Pb- 2 (x) x n+ i) + 
a 

Y ( a + b-2i ) (Pi-i{x) + Pi_ 2 (x) x n+1 ) (p a+b -i-i(x) +p a+b - 2 -i(x) x n+1 ) = 

i=a-\-b—n —2 

= (n + 2 - a) { p a -i{x) p b -i(x) + x n+1 [p a _ x (x) p b - 2 (x) + p a - 2 (x) p b _i(x)] + x 2 n+l p a ~ 2 {x) p b - 2 (x)} + 


Y, (a + b - 2i) pi-i(x) p a +b-i-i(x) - x n +i Y (a + b - 2i) pi-i(x) p a + b - 2 -i{x) + 


i=a-\-b—n —2 


i=a-\-b—n —2 


-x n+1 Y (a + b- 2i) Pi- 2 {x) p a+b -i-i{x) - ; 


i=a+b—n—2 


+. £ 

i—a-\-b—n —2 


(a + b - 2*) p;_ 2 (:r) p a+b _ 2 _j(:r) = 
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a 

= (n + 2 - a)pa-i(x) pb-i(x) - y {a + b - 2i) Pi-\{x) p a +b-\-i{x) + 

i=a-\-b—n— 1 


+ x n +l 



a) [Pa-l(z) Pb-2{%) +Pa- 2{x) Pb-l{ x )] 


a 

{a + b - 2i) Pi_i(x) p a+b - 2 -i(x) + 

i=a-\-b—n —2 


i=a+6— n —1 


(a + b — 2i) pi- 2 {x) p a +b-i-i(x) 


+ 


+ X n+1 < (n + 2 - a) Pa- 2 (x) Pb- 2 {x) - y {a + b-2i) 2 (x) p a +b-2-i(x ) 

v. i=a+b—n—2 

The lower limits i = a + b — n — 2 of the sums in the last expression have sometimes been increases, because when 
there is a factor p a+b -i-i(x) inside the sums, we have the limitation a + b — 1 — i < n, that is, i > a + b — n — 1 , 
otherwise this factor is zero. The terms in this expression for P ab {p(y )) that do not factorize a: n +i and those that 
factorize are the same as those in Eq. (1551) . The terms that factorize must be further transformed to see 
that they also coincide with those in Eq. (1551) . This can be done in the following way: 

a 

(n + 2 — a) [p a -!(x) p b - 2 (x) + p a - 2 {x) p b -i(x)} - y (a + b - 2i) Pi-^x) p a+b - 2 -i(x) + 

i=a-\-b—n— 2 


a 

y (a + b- 2i) Pi- 2 {x) p a+b -i-i{x) = 

i=a-\-b—n— 1 


a— 1 

= {n + 2- a) Pa-i(x) p b - 2 {x) + {n + 2- a) p a -2{x) p b -i{x) - y (a + b - 2i) Pi-\{x) p a +b-2-i{x ) + 

i=a-\-b—n —2 


a— 1 

- {b- a) p a -i{x) p b - 2 {x) - y {a + b-2-2i') pp-^x) p a+b - 2 -i’{x) = 

i'=a-\-b—n —2 


a —1 

= {n + 2-b) p a - i{x) p b - 2 {x) + {n + 2- a) p a - 2 {x) p b -i{x) - y 2{a + b-2i) p t -i{x) p a+b - 2 -i{x) + 

i—a-\-b—n —2 


a— 1 

y ( — 2) Pi- ± {x) p a +b- 2 -i{x) = 

i—a-\-b—n —2 


a —1 

= {n + 2-b) p a -i{x) p b - 2 {x) + {n + 2- a) p a - 2 O) p b -i{x) - y 2{a + b-\-2i) Pi-\{x) p a+b - 2 -i{x ). 

i=a-\-b—n —2 

This expression coincides with the expression in Eq. (1591) inside the curly brackets multiplying x n +i. 

We have so verified that the expression of P ab {p{y )) obtained above, using Eq. (l53l) in the case of .SVi+i coincides 
with the expression of P a b{y) written in Eq. (l59l) . that was obtained from the explicit calculation of P ab {y ) using Eq. 
© and the assumed validity of Eq. (1551) in the case of S n , using the induction principle. This proves that Eq. (l53l) 
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is true for S n for all values of the positive integer n. This holds also for the generating function (l25l) because of the 
equivalence of Eqs. (li?5l) and (15^1) . given the definitions p a = 0, Va < 0 and Va > n, valid for Eq. m- o 


Proof of the generating formula (1261) . 

Theorem [T| tells us how to construct the A-vector Ai det ^ P ^(p) corresponding to the determinant det(P(p)) of the 
P-matrix P(p). 

The positive roots of S n are vectors orthogonal to the reflection hyperplanes, and are for example the ( 2 ) vectors 
e, — ej, 1 < i < j < n. The corresponding linear forms vanishing on the reflecting hyperplanes are then the following: 
Xi — Xj, 1 < i < j < n. Eq. m gives then, for all a = 1 ,..., n: 


X^\p(x)) = 2j2 

gGTl 


S7p a {x) • r g 

lg(x) 


n 

*E 

i<j=i 


S7p a (x) • (ej - ej) 

Xi — Xj 


n 

i<j=1 



f dp a {x) 

V dxi 


dpa(x) \ 

dxj J 


Using the basic invariant polynomials in Eq. (THl) . we see that for a = 1 one has dp g^ = 1, Vi = 1,..., n, so that 
(p(x)) = 0 , while for a > 1 one has: 


that gives, Va = 2 ,..., n: 


dpgjx) 

dxi 


Pa- i(*)U =0 > 


A f et{P)) (p(x)) 


n 

i<j=1 





E 

4 P n — 1 


X, 


—— ( X j Pa- 2 (x)\ x . =x . =0 -Xi Pa - 2 {®)| a 


n 


Xj — Xj 


= -2 E P«-2(*)l*«=* i =o 

*< 2=1 

The last expression is a completely symmetric w-homogeneous polynomial of weight a — 2 in the variables Xi ,..., x n , 
and is formed by terms contained in p a _ 2 {x), so it must be a multiple of the basic invariant polynomial p a _ 2 ( 2 :). To 
find the multiplying factor it is sufficient to count the number of terms, Ci, in the last expression and the number of 
terms, c 2 , in p Q _ 2 (a;). The multiplying factor is then equal to the ratio One easily finds that C\ = ( 2 ) ("- 2 ) anc ^ 
c 2 = ( a ™ 2 )’ so ^at ^ = 5 (n — a + 2)(n — a + 1). One then has: 

X^\p(x)) = -2 i (n -a + 2)(n -a + l)p a - 2 {x) = -(n -a + 2)(n -a+ 1 )p a -i{x ), Va = 2,...,n. 

Eq. (E^l) is then proved. o 


B. Proof of Theorem [3 for A n 
Proof of the generating formula (1321) . 

We adopt here the notation introduced in Section HVl concerning the groups A n and in particular, the definition (l27l) 
of the rotation matrix P n+ 1 , and the derivation of the basic invariant polynomials pi(x'),... ,p n (x'), of A n from the 
basic invariant polynomials Si(x ),..., s n +i(a:), of S n + 1 , as specified by Eqs. (USD (El). 

We know the explicit form of the P-matrix P(s) of S n +±, when the basic invariant polynomials of Sn+i are the 
elementary symmetric polynomials (1281) . This explicit form is given by Eq. (1251) . where the variables p a have to be 
substituted by the variables s a , and n has to be substituted by n + 1. This gives the following expressions of the 
matrix elements of the P-matrix P(s) of SVi.+i : 

min(a,b) 

P a b(s) = [n + 2 - min(a, b)] s a _i s 6 _i - ^ (a + b - 2i) Sj_i s a +b-i-i , Va, b = 1 ,... ,n + 1 , (60) 

2 =max(l,a+b— n— 2) 
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in which so = 1. The matrix P(s) is clearly symmetric. Eq. (1601) is a good point to start with to prove Eq. (1321) . and 
this is the road followed in the present proof. This proof is quite long and will be broken in a few lemmas. 

The P-matrix P(s'), corresponding to the basic invariant polynomials s' a (x') = s 0 (P I [ +1 x'), a = 1,.. ,,n + 1, by 
Theorem [3] item [3], has the same form as P(s), only the variables s a have to be substituted by the variables s ’ a , 
Va = 1,..., n + 1. So, P(s') is also given by Eq. (R)U1) . but with the variables s a substituted by the variables s' a . 

One first proves the following Lemma. 

Lemma 1. The elementary symmetric polynomials Si(x) and s 2 {x) of S n +±, written in Eq. (1281) . are transformed by 
the rotation matrix R„,+ 1 into the following polynomials: 

si(x) ->■ si (a/) = si(Rn + 1 x') = Vn+l x ' n+1 

1 71 

s 2 (x) -+ s' 2 (x') = S 2 {Rl + 1 x') = ^ x'l + i - - J2 X 'i 

i=l 


Proof. Let’s write R for R n+ i in this proof. Let’s first obtain the form of s'^x'). Onehassi(x) = x±+x 2 + . ■ . + x„+i, 
so that, by Theorem [3l itemj^J one has: 

n+1 /n+1 \ n+1 /n+1 \ 

4+) = si(i? T ^) = 5Z E R 1A = £ E R u x 'j 

<=1 \i =1 ) j =1 \ i=1 / 

The sum inside brackets in the last member is the sum of the elements of the j-th row of the matrix R , and by 
definition of R this is 0, except when j = n + 1, in which case it is equal to (n + 1) ^== = \Jn + 1. All this gives 

«!(*') = Vn + l x' n+1 . 

Let’s now obtain the form of s' 2 (x'). One has, by elementary algebraic transformations: 


S 2 (x) = X\X 2 + XiX 3 + ... + x n x n+ i = - 


n+1 

(a^i + X 2 + • • • + #n+ 1) 2 — xj 

i =1 


\ [(si (®)) 2 - INI 2 ] • 


To calculate s' 2 {x') = s 2 (R T x') one has to evaluate the previous expression for x substituted by R T x'. The squared 
norm ||x || 2 = l xf oi x € K n+1 is invariant under orthogonal transformations in K" +1 , that is: ||x || 2 = ||P T x , || 2 = 
||x'|| 2 . Using also the equality s^(x') = Si(P T x') = \Jn + 1 x' n+1 , obtained above, one finds: 


400 = s 2 {R t x') = - [{Vn + l+i) 2 - Ill'll 2 ] = - 


n+1 


(n + 1 ) x ' n+1 


,2 \ ' ,2 

c n+1 - 2 ^ X i 


The Lemma is so completely proved. 


The basic invariant polynomials of A n are obtained by substituting x ra +i = 0 in the polynomials s' a (x'), a = 
2,n + 1. We write x! for the projection in the hyperplane x ' n+1 = 0 of a vector x' S R n+1 and p(x') for the 
restriction to the hyperplane x ' n+1 = 0 of a function p{x') defined in R ™ +1 (typically a polynomial in x'), that is, 
p(x') = p(x') | <+i=0 - 

From the expression of s^ix') recovered in Lemma 1, one obtains 

1 n 

4 +) = ~2 • 

Z i —1 

This expression has not the standard form (|3| one wishes for the quadratic invariant polynomial. To get the standard 
form © it is then necessary a further transformation. The simplest one is to multiply the basic polynomial s' 2 (x') 
by a factor — 2 , for example, by defining s"(x') = s' a (x'), Va = 1 ,..., n + 1 , a ^ 2 , and s'^x') = — 2 s' 2 (x'), but the 
resulting matrix P(s") would then have some fractional coefficients. To avoid fractional coefficients in the matrix 
P{s") it is convenient to multiply all the basic invariant polynomials, not only s' 2 (x by the factor — 2 : 

s a{ x ') = —2 s'(x '), 


Va = l,...,n + 1. 


( 61 ) 
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The P-matrix P(s"), can be determined from P(s') in the way specified by Eq. (fITl) in Theorem [4] The Ja¬ 
cobian matrix J(s') of this transformation is J(s') = —2 7, with I the unit matrix of order n + 1, so one finds 
J(s') P(s') J T (s') = 4 P(s’). The substitution p = p{p’) in Eq. (fT71) . that here becomes s' a = — ^s", Vo = 1,, n+ 1, 
in addition to the substitution of s' a with s", brings a factor —h in front of all variables s' a appearing in the ma¬ 
trix elements of P(s'). In P(s') there are at most quadratic terms in the variables s' 1? ..., s' n+1 , so, after this scale 
transformation the coefficients of the quadratic terms are multiplied by a factor (— those of the liner terms by 
a factor —f, and the constant terms are not multiplied at all. Then, with all these substitutions, all the elements of 
the P-matrix change some coefficients, but P(s") remains without fractions, because of the global factor 4. 

For later use it is convenient to write explicitly the elements of the matrix P(s") corresponding to the basic invariant 
polynomials s"(x') of S n +i (that are obtained from the basic invariant polynomials (l28l) by performing the rotation 
of K . n+1 obtained through the matrix R n + 1 , defined in Eq. m, and the scale transformation (IbTl) '). To obtain the 
matrix P(s"), it is necessary to distinguish the terms containing zero, one or two variables s ' l5 ..., s' n+1 , because they 
are multiplied by different coefficients, so, it is also convenient to distinguish the element Pn(s"), the other elements 
of the first row and column of the P-matrix P(s"), and the rest of the P-matrix P(s"). In this last case, using 
Eq. (1601) . with a, b > 2, a linear term in the variables s^,..., s' n+ i appears only if there is an i for which i — 1 = 0, 
because then one finds the variable Sq = 1. When i = 1, the variable s' a , b _i_i = s ' a+b _ 2 is not equal to 0 only if 
a + b — 2 < n+1, that is when a + b < n + 3. (The only other possibility is when a + b — 1 — i = 1, that is when 
i = a + b — 2 , and the limitations imply that that is possible only if* = a = fe = 2 , but in that case its coefficient 
vanish: a + b — 2i = 0). We can then easily find, from Eq. (1601) (in which the variable s a is substituted by the variable 
s' a , Va = 1,... ,n + 1), using Theorem [4] (that implies the multiplication of all the terms by factors 4 • 1, 4 • (— i), 
4 • as specified above), the following expressions for the elements of the P-matrix P(s"): 

Pn(s") = 4 (n + 1), 

Pla(S') = Palis") = -2(n +2 

Pab{s") = [n + 2 - min(a, b)\ s"_ x s£-i + 2 (a + b - 2 ) s " +b _ 2 + 

min (a,b) 

(a + b-2i) s"_! s" +fe _i_ i , 

2 —max (2, a+6—n—2) 

The term with s ”, b _ 2 is present only ifa + 6 < 7 i + 3, that isifa + & — 2 < ti + 1. To take into account this condition 
it is convenient to define, in addition to Sq = 1, that was required by Eq. (ED, also s" = 0, Va > n + 1. We can 
observe that in Eq. (l62l) the variable Sq cannot be present, because of the limitations, so it is simpler to consider for 
Eq. (l62l) only the constraints s" = 0 ,Va> 7 i + l. 


Vo = 2, ...,7i + l, 

Va, b = 2 ,..., n + 1 , (62) 


The two transformations performed up to now commute. It is possible first to rotate the system of coordinates 
and then to multiply by —2 the basic invariant polynomials or vice versa, first to multiply by —2 the basic invariant 
polynomials and then to rotate the system of coordinates. In both cases one ends with the same expressions both for 
the basic invariants s^a/), and for the matrix P(s ,/ ). Actually, only the multiplication by —2 is sufficient to obtain 
the form (1621) of the P-matrix. 


The basic invariant polynomials of A n are easily found from the basic invariant polynomials a = 1, ..., n + 1, 

by taking x' n+1 = 0. This leaves the n basic invariant polynomials s"(ir'), a = 2, ..., n + 1, of degrees 2,3 .,n, only, 
because s"(x’) = 0. Moreover, the quadratic invariant s 2 (x') has the desired standard form ([3D, equal to the squared 
norm: s 2 (x r ) = || x '\\ 2 . A renaming the basic invariant polynomials is then convenient, like in the following way: 

Pa(x') = s" + i(i'), Va = 1,..., n. (63) 


The n x n matrix P(p) can be obtained through Eqs. ([HD and (HOD from the basic invariant polynomials p a (x'), 
a = 1,..., ti, written in Eq. (ED, but this is not what we want to do here. The same matrix P(p) can also be 
determined from the (n + 1) x [n + 1) matrix P(s"), and this is what we are going to do, but with caution. The first 
idea is to substitute s" = 0 in the matrix P(s") (because s"( x') oc x' n+1 , and when x' n+l = 0, one has s"(x’) = 0), 
to remove its first row and column (because there is no more an invariant of degree 1 , and the first row and column 


of the matrix P{s") 


have no more meaning), and to rename the variables as suggested by Eq. (1S51) : s " +1 = p a , 

s'/=0 


Va = 1,..., n. However, doing that way one generally does not find the right P-matrix P(p), corresponding to the 
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basic invariant polynomials p a (x'), a = 1 that can also be calculated from its definition, that is from Eqs. 

and m- We are now going to see when and how the knowledge of the matrix P(s") can determine the matrix P(p). 

Let’s define a matrix operator [ • ] p that applied to a (n + 1) x (n + 1) matrix P(s"), dependent onn+1 variables 
s'{,... , s" +1 , gives annxn matrix [P(s")] p , dependent on n variables pi,... ,p n , obtained from the matrix P(s") by 
applying the following rules: 

1. make the substitution s" = 0 in the matrix P(s"); 

2 . delete the first row and column of the resulting matrix; 

3. rename the variables like in Eq. CH): s " +1 =p a , Va = 1 ,.. .,n. 

One would desire that [P(s")] p = P(p), with P(p) the P-matrix corresponding (through Eqs. © and (fIU|) l to the 
basic invariant polynomials of A n in Eq. (1631) . However, this is in general not true, because the elements of the matrix 
[P(s")] p can be obtained by making the scalar products of gradients in R n+1 , Vs" +1 W) ■ Vs ' 6 ' +1 (x'), V a, b = 1,..., n, 
by expressing the results in terms of the basic invariant polynomials s"(aV), a = 1 ,... ,n + 1 , and then by doing 
the substitutions s"( x') = 0, s" + 1 (P) = p a , Va = 1 ,...,n, while the elements of the matrix P(p) can be obtained 
by making the scalar products of gradients in R n , Vp 0 (x') • Vpb(ir'), Va, & = 1 by expressing the results 

in terms of the basic invariant polynomials p a (x '), a = 1 ,... ,n, and then by doing the substitutions p a {x') = p a , 
Va = 1,.. . ,n, and these two procedures give generally different results. There are, actually, infinitely many choices 
of a set of n + 1 basic invariant polynomials z a (x') of S n +i, a = 1 ,..., n + 1 , with zi(x') oc x' n+1 , that give rise, by 
taking x ’ n+1 = 0, to the same set of n basic invariant polynomials p a (x') of A n , but that give rise to different matrices 
\P(z)] p . So, only clever choices of the basic invariant polynomials z a {x') are such that P{p) = [P(z)] p , as one desires. 

Let’s prove the following Lemma. 


Lemma 2. A set of basic invariant polynomials of S n + 1 , z a (x'), a = 1, ... ,n + 1, with Z\{x') oc x' n+1 , by taking 
x' n+ i = 0, can determine a set of basic invariant polynomials of A n , p a (x') = z a +i(x'), a = 1 ,...,n, with the 
property that P(p) = [P{z)] p , where P(p) and P{z) are determined, through Eqs. JH]) and (TTU1) . by the sets of basic 
invariant polynomials p a (x'), a = 1 ,n, of A n , and z a (x'), a = 1, . .. ,n + 1, of S n ., respectively, if and only if 


Pla(z) 


= 0, Va = 2,..., n + 1. 

zi—0 


Proof. Let’s use the symbols V(„) and V(„ + i) for the gradients in R” and R n+1 , respectively. Adopting the usual 
notation x’ € R n+1 , and x' = x'\ x , i=0 £ R n , one then has: 


( d 

9 ) 

_ ( d 


\9x[- 

"’dx'J ’ 

(n+1) 

"dx' n+ J 


Consider the matrix element P a b(p), V a, b = 1 , ,n. One has: 


Pab(p(x')) = V( n) p 0 (: r') • V {n )Pb(x') • 

If the equality P{p) = [P(z)] p is true V a, b = 1 ,..., n, the matrix element P a b{p(-%')), after making in it the substitution 
Pa(x') = z a + i{x'), Va = 1,..., n, must be equal to the following matrix element: 


Pa+l,b+l (^'(•^ )) 


,=0 


(V(„+i)2 0+ i(a;') • V(„ + i)Z b+ i(a: , ))| <+i=0 = 


_ , (dz a +i(x' 

— v (n)Za -\-1 ) ' v (n)Zb-\-l\X ) H - I ^ , 

\ Cfx n+1 


dz a +l(x') dz b+1 {x') 


dx' n+1 


x ' n +i=° 


The first term in the last member, with z a +l(x') substituted by p a (x'), Va = 1 ,...,n, is equal to V( n )p a (a: , ) 
V( Jl )Pb(x / ) = P a b(p(x')), and, for the hypothesis P(j>) = [P{z)] p , is also equal to P a+ \^+\{z{x')) . So, if the 

equality P{p) = [P(z)] p is true, one has: 

' dz a+ i(x') dzb+i(x') 


dx 


n+1 


dx: 


n+1 


= o, 


Va, b = 1 ,...,n . 
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In particular, by taking a = b one finds: 

dz a+1 (x') 

f)nr f 
UJj n +1 

From the hypothesis z\{x') oc x' n+1 , one has: 


= 0, Va = l, 


dzi(x') 


V (n-f- 1 ) Zl(x ) — &n+l — C e n+ l , 


(64) 


(65) 


n+1 


where e n +i is the canonical unit vector of R” +1 orthogonal to the hyperplane x' n+1 = 0, and c is a non-zero real 
number. One then obtains, using Eqs. (1^51) and (1M1) . Va = 1,... ,n: 


P 


,a+l(z(x)) ' — iy (n+l)Zl{x ) • V( n+ 1)Z Q +i (lE )) \ x , Q — (c ^ 

X n+1~ u n+L \ u ' L n -\-1 / 


= 0 . 


This condition is equivalent to P\ a {z) = 0, Va = 2,..., n + 1, because z\(x r ) oc x' +1 implies that the condition 

21=0 

x' n+1 = 0 is equivalent to the condition zi(x') = 0. The necessity (only if part) of the statement of the Lemma is 
proved. 

Vice versa, supposing P la (z) = 0, Va = 2 ,...,n + 1, that, because of z\(x') oc x' n+1 , is equivalent to: 

2i=0 


Pl,a+l(z(x')) 


<+l=° 


= 0, Vo = l,...,n, 


one obtains the following conditions: 

{y(n+i) z i{x') ■ V (n+1) 2: Q+ i(a;'))| a .j i 0 = 0, Va = 1,.. ,,n. 

For the hypothesis Z\{x') oc x' n+1 , that implies Eq. (l65l) . the previous equation implies Eq. (l64l) . but then, Va, b = 
1 ,..., n, one has: 


P 


a+1,6+1 


(*(*')) 


= y(n+i)Z a +i{x') ■ V (n+1) z ( , + i(a; , ))| , =0 = =0 = 


= V (n) 2 a+ i{x') • V (ra) 2 b+ i(x') = V^p^x') • V (n) p 6 (x') = P a b(p(x ')), 

with Pa(x') substituted by z a +i(x'), Va = l,...,n, and this implies P(p) = \P(z)] p . The sufficiency (if part) of the 
statement of the Lemma is proved. o 


Now we have to find the transformation sending the basis s'^{x’), a = 1,... ,n + 1, of invariant polynomials of 
S n + 1, to a basis z a {x '), a = 1,..., n + 1, of invariant polynomials of SVi+i that satisfies the hypothesis of Lemma 2, 


that is: P\ a {z) 


z\ —0 


= 0, Va = 2,..., n + 1. To this end, one proves the following Lemma. 


Lemma 3. The basic invariant polynomials z a (x'), a = l,...,n + 1, of SVi+ij that are obtained from the basic 
invariant polynomials s"(x'), a = 1 ,..., n +1 , of SVi+i (that are defined by Eqs. (P251) . (1^1) . and (RJT1) 1. in the following 
way: 


zi(x') = s'l(x'), 

Za(x') = s"(x') + C a s"(x') s"_i(x') , C a = ^ ^ 2 ^ 

2 (n + 1) 


Va = 2,...,n+1, 


( 66 ) 


have the property that [P(z)] p = P(p), where p a { x') = z a + i(x'), Va = 1,..., n, is a basis of invariant polynomials of 
A n , and the matrix P(p) coincides with their P-matrix, as can be calculated from Eqs. (0 and (flOl) . The only other 
possibilities to find basic invariant polynomials z a (x r ), a = 1,..., n + 1, with the property that P(p) = [P(z)] p , is 
obtained by adding terms like s" 2 (x') q a - 2 {s" (x')) to the polynomials z a (x'), defined in Eq. (Ififill . in which a > 2 and 
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q a - 2 (s"(x')) is an arbitrary invariant homogeneous polynomial of degree a — 2 . 


Proof. Because of Lemma 2, to prove that the basic invariant polynomials z a (x'), a = 1, ...,n+ 1, satisfy the 
equation [P(z)] p = P(p), it is sufficient to prove that they satisfy the conditions P\ a {z) = 0, Va = 2 ,..., n + 1. 

zi =0 

As zi(x') = s'/(s') oc x’ n+1 , one has zi(x')\ , =0 = s"(x') | , =0 = Zi(x') = s"(x') = 0, and Pi a (z) = 

n + i n + 1 Zl =0 

Pia{z(x')) , Va = 2,..., n + 1. Using the transformation (l66l) . and the matrix P(s") written in Eq. (l62l) . one 

<+l=° 

finds, Va = 2,..., n + 1: 

Pia{z{x')) ^ = (Vzi(a;') • Vz a ( x'))\ x , n+i=0 = (X7s"(x') • V [s"(x') +c a s'/(a;')s"_ 1 (a;')])| <+i=0 = 

= (Vs'/(a:') • Vs"0r'))| <+i=o + c 0 Vs'/( x') • Vs'/K)) ^ + c a («?(*') Vs'/(a:') • Vs".^')) | x; 0 = 


= Pia(s"(x')) +c a s"_ x (a;') Pn(s"(i')) + 0=-2(n + 2-a) s"_ 1 (x') + c a s"_ 1 (a;')4(n + 1) = 


— ( -2(n + 2 — a) + ———— 4 (n + 1) ) s"_ x (x r ) — 0 . 


The conditions just proved imply Pi a (z) 


z i=0 


2(n+l) 

= 0, Va = 2,..., n+1, and Lemma 2 assures that the equation [P(z)] p = 


P(p) is true. An additional term like s'{ (x') q a - 2 (s"(x')), added to a polynomial z a (x') in Eq. (1M1) . a = 2,..., n + 1, 
implies the following additional term to the gradient Vz a {x'): 2s "(x') q a - 2 {s"(x')) Vs "(x 1 ) + s'{ 2 (x') Vg a - 2 (s"(i')), 
that, however, vanishes when s'{{x') = 0, and hence originates the same matrix [P(z)] p . The same additional 
term s'/ (x')q a - 2 {s"(x')) vanishes when x ' n+1 = 0 and has no influence on the definition of the polynomials p a (x'), 
Va = 1,..., n, too. o 


It is interesting to note that the basis transformation (1M1) was used also in Section (2.4.8) of Ref. [3^], in¬ 
side a proof concerning the P-matrices of the groups of type A n . although the proof in Ref. :39 ! ] was done in 
view to prove the existence of a flat basis (recall the discussion with Eq. (1201) 1 for the groups of type A n , and not 
to find a generating function for the complete structure of the P-matrices for the groups of type A n , as the present one. 


We should now determine the matrix P(z), corresponding to the basic invariant polynomials of S n +i defined 
in Eq. (1661) . because from the knowledge of P(z) it is possible to determine the matrix P(p) from the equation 
P(p) = [P(z)] p . The matrix elements P a b{p ) would then be equal to the matrix elements P a +i^+i(z) , with 


zi—0 


the variables z 0 +i substituted by the variables p a , Va = 1 ,n. The determination of the matrix P(z) is however 
somewhat difficult, and we can simplify the calculation if we determine in one step the matrix P(p) from the matrix 
P(s"), without finding first the matrix P(z). This is done in the following Lemma. 


Lemma 4. The P-matrix P(p) corresponding to the basic invariant polynomials of A n p a (x'), a = 1,..., n, restrictions 
to the hyperplane x' n+l = 0 of the basic invariant polynomials z a (x'), a = 2,..., n + 1, of S n+ 1 , specified in Eq. m, 
and of the basic invariant polynomials s"(x'), a = 2,.. .,n + 1, of S n +i, specified in Eq. m, that is p a {x') = 
z a - fi(x') = s", i(s'), Va = 1 is the following one: 


Pab(p ) = 


i (a,b) - 


n + 1 


min(a,6) — 1 

Pa -1 Pb-1 + 2 (a + b) Pa+b-1 - ^2 (a + 6 - 2i) p a+b _i_i, 

i=max(2,a+b— n— 1) 

Va, b = 1,..., n , (67) 


in which one has to consider po = 0 and p a = 0, V a > n. 


Proof. For Lemma 3, the matrix elements P a b(p), V a,b = 1,..., n, are equal to the matrix elements 
a, b = 1,... ,n, with the variables z a + i substituted by the variables p ai Va = 1 ,,n. 


Pa+l,6-f-l (-^) 


zi=0 
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It is easier to determine the matrix P(p(x')) through Eqs. (0 and (flUl) . instead of using Theorem [4] The explicit 
form of the matrix P(s"), written in Eq. (1621) . will also be used. 

The gradients of the basic invariant polynomials z a+ i(x’), Vo = 1,..., n, are calculated, from Eq. (lC)Sl) . in terms of 
those of s" +1 (x , ) ) Va = 0,...,n, and are as follows, Va = 1,..., n: 

Vz a+1 (x') = V K+i^) + c a + i s"(x') s"(x')\ = Vs" +1 (x') + c 0+ 1 s"(x') Vs" ( x') + c a+ 1 s"{ x') Vs"(x'), 

and when x' n+1 = 0, that imply s"(x') = Z\(x') = 0, they become: 

V z a +i{x')\ x , n+i=0 = S7 sl +l {x')\ x ,^^ 0 + Ca+i s'^x) Vs"(x') \ x ' n+1 =o > Vo=l,...,n. 

We know that [P(z)) p = P(jp), so we know that: 

(V4+l(l') • V2i,+ i(x')) 

Then, Va, b = 1,...,n, one has: 


= (V«a+i(® , )| x ; +1=0 ) • (V2 6 + i(x')| <+i= o) • 


a+ 1 , 6+1 \ 


{z{x')) = (Vz a +i(x') • Vzfe+i^')) 

x' , . —0 V / 


= (V S a+l ( X ')\ x ' n+1 =0 + C o+l S 'a{x') Vs"(x') |^ +1=0 ) ' ( Vs b +1 (%') L ; +1 =0 + C 0 -l S b( x ') Vs l ( x ')lx' n+1 =o) = 

= P a+1:b+1 (s"(x')) + C b+ 1 S b (x') P 1>a+1 (s"(x')) + C a+ 1 s"(x') A , 6+1 (s" (x')) + C a +1 C b+1 s"(x') S b (x') Pl,l(s"(x')) = 


s7=0 


— Pa+l,6+l(s ,/ ) + C b +1 s" -Pl 5 a+l(s") + C a +i s" Pl,6+l(s") + C a +i C b +i s" S b Pl,l(s") 


sy — 0 


( 68 ) 


where in the last member the dependence on x' has been understood, and the condition x' n _ t _ 1 = 0 has been replaced 
by the equivalent condition s” = 0. 

The terms inside the square brackets can be simplified by substituting the definitions of the coefficients c 0 , given in 
Eq. (1661) . and the explicit form of the matrix elements of the first row and column of the matrix P(s"), given in Eq. 
m- In this way the content of the square brackets becomes as follows: 

n + 1 - b „// S r //x , n + 1 - a // S ( , ( n + 1 - a )( n + 1 - b ) n n u < 

S b P h a+l(s ) + —, - ——S a P 1 , b+1 {s ) + - — - —2 - S a S b P hl {s ) 


■■■ Jls i =0 [2(n + l) 

n + 1 — b 


2(n+l) 


4 (n + l) 2 


s7=0 


2 (n + 1) 


„ , n+1 — a ,. ... (n + 1 — a)(n + 1 — b) 

s b ( — 2(n + 1 — a) s a ) + ———— s a (—2(n + 1 — b) s b ) H ————^- s a s b 4 (n + 1) 


2 (n+1) 


4 (n + l) 2 


s '-/=0 


(n + 1 — a) (n + 1 — 6) 

l s a s b ) ls'.'=0 ) 


n+1 

expression that is equal to 0 if a = 1 or b = 1, because s" = 0. 

To simplify the explicit expression of the matrix element P a +i,6+i(s") 


s'/=0 


, that appears in the first term of the last 


member of Eq. (ESI), one has to take into account the third line of Eq. (1621) and the conditions valid for Eq. (1621) . 
that is: s" = 0, Va > n + 1. One so finds, Va, b = 1,..., n: 

min(a,b)+l 

Pa+I,b +1 (s") = I [n + 1 - min(a, b )] s" s" + 2 (a + b) s" +b - ^ (a + b + 2 - 2 i) s-'.j. s"+6+i-i 

s/=0 \ 

i=max(2, a-\-b—n) ) s //_q 


min(a,6)+l 

= [n + 1 - min(a, b)] (s" s b )| s ,, =0 + 2 (a + b) s" +6 - ^ (a + fo + 2-2f) (s"_! s" +b +i-i) | s » =0 = 

i=max(3, a+b—n) 
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= [n + 1 - min (a, b )] (s" s b )| s ,, =0 + 2 (a + 6 ) s" +& - [a + b - 2 min (a, 6 )] (s" in(a6) s" ax (a,b)) 


, =o 


min(a,fe) 

F, (a + + 2 — 2 i) s"_ 1 s" +b+1 _ i = 

i=max(3,a+6— n) 


= [n + 1 


min(a,6) — 1 

max(a, 6 )] (s" s b )L"=o + 2 ( a + b ) s "+b ~ Y ( a + b 

i'=max(2,a+6—n—1) 


20 4 


where in the third expression the lower limit of the sum has been increased from 2 to 3 because the condition s'{ = 0 
implies s"_i = 0, if i = 2, and where in the last expression one one could consider s" in ^ a b j s" ax ( a b ) = s" s b . We can 
then write for Eq. (EHD the following expression, where, as before, the dependence on x' is understood: 


zi—0 


P a+ i,b+i{z) = (n+ 1 - max(a, b) - 


(n + 1 — a)(n + 1 — b) 


^ ( s a S 6)ls"=0 + 2 (a + b) s " + 6 + 


min(a,6) — 1 

(a + b - 2 i) s" s" + b-i = 

i=max(2,a+b— n— 1) 


min(a, b) — 


a b 

n + 1 




,'=o + 2 


(a + &) 


s a+b — 


min(a,b) — 1 

Y (a + b — 2i) s" s" +b -i 

i=max(2,a+6—n—1) 


Vo, b = 1, 


In the last member, one has to substitute the polynomials Va = 1 , ,n + 1, with their expressions in terms 

of the polynomials z a {x'), Va = 1,... ,n + 1, that can be obtained from Eq. m, but for s"(x') = 0 this is easy, 
because Eq. (El gives z”(x') = 0 and s"(x') = z a (x'), Va = 2,..., n + 1. To find the matrix P(p ) one has, finally, to 
substitute p a in place of z a +± (that is of s" +1 ), Va = 1,... ,n. One so finds: 


Pab{p) = 


min(a, b) — 


a b 


n + 1 


min(a,6) — 1 

Pa -1 Pb-1+2 {a+b) Pa+b-i- Y, (a+b-2i) Pi -1 Pa+b-i-i 

i=max(2,a+6— n— 1) 


Va, b = 1, 


where the conditions z\ = 0 (corresponding to x' n+1 = 0) and z a = 0, V a > n +1, become po = 0 and p a = 0, V a > n. 
The proof of the Lemma is complete. o 


The matrix P(p) in Eq. (l67l) has fractional coefficients, because of the denominator n + 1. If one wishes to have 
a P -matrix with only integer coefficients one can make a scale transformation of the invariant polynomials p a {x'), 
a = 2,... ,n, avoiding to modify the quadratic invariant pi{x') that has its standard form given in Eq. Ej). This is 
done in the following Lemma. 

Lemma 5. The basic invariant polynomials of A n p' a (x'), a = 1,..., n, defined by the scale transformation 

Pa(x') = ( n + l )^ Pa(x '), Va = 1,..., n, (69) 

in terms of the basic invariant polynomials of A n p a (x r ). a = 1,..., n, given in Lemma 4, are such that p'i(x'), has the 
standard form of Eq. 0 : Pi( x ') = ||a/|| 2 , and the corresponding P-matrix P(p') can be obtained with the following 
formula: 

Pab(p') = [(n + l) mm(a,b) - ab}p' a _ 1 p ' b _ 1 + 2 (a + b) p ’ a+b _ 1 + 

min(a,6) — 1 

- O+l) Y (.a + b- 2 i)p’ i _ x p’ a+b _ x _ i , 

i=max(2,a+6-n-l) 


Va, b = 1 ,..., n 


(70) 
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where one has to consider p' 0 = 0 and p' a = 0, V a > n. The matrix P(p') has then only integer coefficients. 

Proof. We first note that for a = 1, Eq. (1551) gives p'\{x') = pi(x'), hence p[(x') = ||x'|| 2 . To obtain the matrix P{p') 
we use Theorem 3) The Jacobian matrix of the transformation (1551) is diagonal: 

Jab(p) = (n + 1)~ Sab, Ma,b= 

that implies J T (p) = J{p). The substitution p = p(p') in Eq. 1171) . using Eq. (1551) . becomes the following one: 

Pa = (n + l)^ p' a , Va = 


Using Eq. (1571) . Va, b = 1,..., n, Eq. (1171) gives: 

Pab{p') = Jac(p) Pcd(p) Jj b {p) = ((n + 1)“ S ac Pcd(p) (n + l)“ S db \ 

P=P(p') V ' 


P-P(P') 


= (n+ 1) +2 P ab (p ) 


p=p{p') 


= {n + 1) ; 


1 (a, 6) - 


n + 1 


{n + 1) ' 2 ' p' a _ x (n + 1) 1 <2 11 p' b _ x + 2 (a + b) [n + 1)* ^ 


min(a,6) — 1 


E 


(a + 6 - 2i) (n + 1 ) ' ^ (n + 1 ) ' ( “ +2 ' ° ^ +6 _i_i > = 


z=max(2,a+£>— n— 1) 


= (n + l)“ 2 < [(n + 1) min(a,6) - ab] (n + 1) 2 p ' a _ 1 p ' bl + 2 (a + 6) {n + 1) 2 Pa +b _i+ 


min(a,b) —1 


-(^ + 1 )" 


(a + 6 - 2*) Pa+b-i-* / = 


i=max(2,a+6—n—1) 


min(a,6) — 1 

= \(n + 1) min (a, b) - a b\ p' a _ x p' b _ x + 2 (a + b) p' a+b _ x - (n + 1) ^ (a + 6 - 2 i) p\_ x p' a+b _ x _i ■ 

i=max(2,a+b—n—1) 

This is the P-matrix in Eq. Con and is clear that all the coefficients of the monomials in the variables p' a , a = 1,..., n, 
that appear in it, are integer numbers. o 


The P-matrix in Eq. m is the P-matrix corresponding to the basic invariant polynomials p' a (x'), a = 1,..., n, of 
A n , and is the one also reported in Eq. (l32l) . Let’s review how, in this proof, the polynomials p' a (x'), a = 1,... ,n, 
are obtained from the basic invariant polynomials s a (x), a = 1,..., n + 1, of S„+i, given in Eq. (1751) . 

1. We performed a rotation of the space R” +1 by the matrix R n+x given in Eq. m- x’ = Rx. The basic 
polynomials of S„+i were transformed into s' a (x') = s a (Rn + ix), Va = 1 ,...,n + 1. 

2. We performed the scale transformation by the factor —2 as written in Eq. m- The basic polynomials of S n +i 
were transformed into s"(x') = —2 s' a (x'), Va = 1 ,..., n + 1. 


3. We performed the basis transformation written in Eq. (1551) . 
into zi(x!) = s'/O'), z a (x') = s"(x') + s'/O') 


The basic polynomials of S n +i were transformed 
Va = 2 ,...,n + 1. 


4. We took the restrictions to the hyperplane x' n+1 = 0 of the basic invariants z a {x'), Va = l,...,n + 1, of 
S n+ i. One would obtain the same polynomials by taking the restrictions to the hyperplane x' n+1 = 0 of the 
basic invariant polynomials s"(x'), Va = l,...,n + 1, of S n+ 1 (although [P(z)] p ^ [P(s")] p ). One so finds 
s'/O') = s'/OOIx' =o = 0, and the restrictions to the hyperplane x' n+l = 0 of the non-linear polynomials 
s"(x'), a = 2,..., n + 1, form a basis of invariant polynomials of A n . These basic invariant polynomials of A n , 
after renaming and renumbering, are specified in Eq. (1551) : p a {x') = s" + i(x') = s" +1 (x')| a; , =Q , Va = 1,...,n. 
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5. We performed the scale transformation reported in Eq. (1691) . in which the basic invariant polynomials of A n 
were transformed into p' a {x') = (n + l )^ 1 Pa(x'), Va = l,...,n. 

Regarding the construction of the basic invariant polynomials of A n p' a (x'), a = 1 ,..., n, the relevant transformations 
are those in items 1., 2., 4., and 5. The order is not important, except for the restriction to the hyperplane x' n+1 = 0 
(item 4.) that must be performed after the rotation (item 1.). For example, one obtains the same basic invariant 
polynomials following the operations in the following order: item 2., item 5, item 1., item 4., or in the following one: 
item 1., item 2., item 5., item 4., that is the one reported in Section HVl The combination of the scale transformations 
in items 2. and 5., if done before the renaming of item 4, results in the following scale transformation: 

t a (x') = -2 (n + l)~ s' a (x '), a = 1 ,... ,n + 1 , 

where the index a in Eq. (1691) has been changed into a — 1 because the scale transformation is done before the renaming 
of item 4.. The scale transformation has been extended also to s\ (x r ) without any problem, because the restriction to 
the hyperplane x (j , 1 = 0 implies in any case that t\(x') = 0. This scale transformation is that one reported in Eq. m- 

Without the restriction to the hyperplane x' n+ i = 0 specified in item 4. (and with or without the rotation specified 
in item 1 .), one obtains a set of basic invariant polynomials of S n + 1, but both the invariant polynomials and the 
corresponding P-matrices depend on the order these transformations are performed. As we have seen the operator 
[ • ] p , when applied to the P-matrix corresponding to the transformed invariant polynomials of Sn+i, gives a P-matrix 
for A n , only if the hypothesis of Lemma 2 are satisfied. 

We started this proof from the P-matrix of S^+i in Eq. (IflOl) . The final matrix reported in Eq. Col), is the P-matrix 
corresponding to the basic invariants p' a (x'), a = 1 ,..., n, of A„, obtained in this proof and coincides with the matrix 
reported in Eq. m, that is the P-matrix corresponding to the basic invariants p a (x'), a = 1,..., n of A n obtained 
in Eqs. (128)) - (I3T1) . This completes the proof of the first part of Theorem [7] o 


Proof of the generating formula (1331) . 

The basis of invariant polynomials p a {x'), a = 1,... ,n, of A n , to which it corresponds the matrix P{p) of A n that 
can be generated by Eq. (l32l) . was derived from the basis of invariant polynomials s a (x ), a = 1,..., n + 1, of S n + 1 , 
that is written in Eq. (l28l) . through the transformations defined in Eqs. (l29l) . (l30l) and (1311) . We shall now determine 
the A-vector corresponding to the determinant det(P(p)) of the P-matrix of A n , using TheoremlU starting from the 
A-vector corresponding to the determinant det(P(s)) of the P-matrix of S n +i, that is written in Eq. (l26l) . but in which 
n must be substituted by n + 1, and the variables p a must be substituted with the variables s Q , Va = 1,..., n + 1, 
that is from the following A-vector of det(P(s)): 


A (det(P)) (s) = (A (det(P)) (s)j . . . 5 Ai d6t ( P »( S ) 

in which one has to consider s_i = 0 and sq = 1. 


— (n — a + 3)(n — a + 2) s a - 2 , o = 1,...,n+1, (71) 


With the basis transformation written in Eq. 


has, by Theorem [3j item [3], P(s') 



consequent to the rotation produced by the matrix R n + 1 , one 
and Theorem [5l item[l] asserts that 


A (det(P))( a /) = A (det(P))( s ) 


so that A^ det ^^(s') is given by Eq. (ITTl) . with s substituted by s'. 

To the basis transformation of the basic invariant polynomials of S'n+i written in Eq. (1301) . it corresponds the 
(n + 1) x (n + 1) diagonal jacobian matrix J(s'), whose elements are the following: 

{Af 

Jab{s') = = -2 y/(n+ l) a ~ 2 Sab , a,b = l,...,n + l, 


with S a b the Kronecker delta. 

The determinant of the P-matrix of S n +i P(t), corresponding to the basic invariant polynomials ti(x'),... ,tn+l(x'), 
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satisfies the boundary equation (fl2l) with the A-vector Aa det ^ P ^(f) that can be determined from A^ det ^^(s') using 
Theorem^ item[2] This calculation gives, Va = 1,..., n + 1: 


A(det(P))( t ) = j Qh(s /) 


s'=s'(t) 


= -2 ^(n + l) 0 - 2 S ab Ar t(P)) (s'; 


'=?(*) 


= -2 V(n + 1 )“- 2 Ai de ‘( p » (s') ^ = 2 V(n+ l ) a " 2 (n - a + 3) (n - a + 2) s' 

s'=s'(i) 

The inverse transformation of the transformation in Eq. (l30l) is the following one: 

1 


s'=s'(i) 


s'a(t) = -- yj{n + l ) 2 “i„, a = 1 ,..., n + 1 , 


that gives S a_ 2 | s ' = p( t ) = s'a- 2 (t) = —\^/(n + l ) 4 a t a - 2 ■ In the formula for X^ et( ~ P ^ (t), we have however to 
consider separately the cases a = 1 and a = 2 from those with a > 2 , because when a = 1 in the formula one finds 
s' a _ 2 = s'^ = 0 , and when a = 2 one finds s ' a _ 2 = s' 0 = 1 , so in both these cases no substitution has to be performed 

because no variables s' a , a = 1,. .., n + 1 are present in the expression for Ai det ^ P ^ ( t ). One finds then: 

A( de t ( p)) (t) = 0 5 A (det ( p)) ( t) = 2 n ( n + , 


and Va = 3,..., n + 1, one lias: 

A (det(P))(£) _ 2 ^/(n+ l) a ~ 2 (n - a + 3)(n - a + 2) i/(n + l) 4_a t a _ 2 ^ = — (n + l)(n - a + 3)(ra - a + 2) t a _ 2 . 

The last transformation to do is that one in Eq. ED, that can be performed in two steps. In the first one one has 

to take x ' n+1 = 0, that means ti(x') = 0, and implies A ( dct ( p ))^ = 0. In the second and last step one has to rename 
the basic polynomials and the variables in the following way: t a + i(x') = p a {x '), a = 1 because t i(x') is no 

longer an invariant polynomial, so that it does no longer interest, A^^ P ^(t) = Ai det ^ P ^(p), Va = 1 and this 

gives the following final form of the A-vector of det(P(p)), where its first component has been dropped: 

A ( d et(P)) (rt = (2n (n + 1), A^ det(P)) (p), ..., X^Hp)) , 

A^ det(P)) (p) = -(n + l)(n - a + 2 )(n - a + 1 )p a -2 , Va = 2 ,... ,n, 

in which one has to consider po = 0. This is Eq. EH), so the proof is complete. Using p 0 = 0 one sees that 
A( det(p)) (p) = 0. o 


C. Proof of Theorem [8] for B n 
Proof of the generating formula (1351) . 

The proof can be obtained by induction on the rank n. The lowest possible value for n is n = 2 (because of the 
isomorphism B\ — A\), so we first prove that Eq. (1351) is true for n = 2. 

The basic invariant polynomials of P 2 , from Eq. ED), are the following: pi(x) = x\ +x 2 , p 2 (x) = xfx 2 , their gradi¬ 
ents are the following: Vpi(x) = (2xi, 2x 2 ), Vp 2 (x) = (2xix|, 2x 2 , x 2 ), and the elements of the matrix P(x), calculated 
with Eqs. El) and (TTU1) . are the following: Pn(x) = 4(x 2 + x 2 ) = 4pi(x), Pl 2 (x) = P 2 i(x) = 4(xfx 2 +xfx 2 ) = 8 p 2 (x), 
P 22 (x) = 4{x\x\ + xfx 2 ) = 4pi(x)p 2 (x). Using Eq. (1551) one finds the following expression for the P-matrix elements 
of B n : Pn(p) = 4pi, P i 2 (p) = P 2 i(p) = 8 p 2 , P 22 (p) = 4pip 2 . We see then that for n = 2 Eq. ED is true. 

Using the induction principle, we suppose that Eq. (l35l) is true for a general n > 2, and prove in all generality that 
it is true also for n substituted by n + 1. This implies then that Eq. (l35l) is true for all n > 2. 
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We use the symbols x and y for general vectors of 1" and R n+1 , that is: x = (xi, X 2 , ■ ■ ■, x n ), and 
y = (xi, X 21 ■ ■ •, x n , x n +i). B n acts in R" and its basic invariant polynomials are p a {x), a = 1 defined 

in Eq. (j34]), and B n+ 1 acts in M" +1 and its basic invariant polynomials are p a (y), a = 1,... ,n + 1, also defined in 
Eq. (1551) . but with n substituted by n + 1. 

The basic invariant polynomials of B n +± can be expressed in terms of those of B n in the following way: 

Pi(y) = Pi(x)+xl +1 , 

P2{y) = P 2 (x) + Pi(x) xl +1 , 

P3 (y) = P3 (x) + P 2 (x) xl +1 , 

Pn(y) = Pu{x) +Pn-l(x)xl +1 , 

Pn+l(y) = Pn(x)xl +1 , 

By defining po(x) = 1 and p n + i(x) = 0, one has the following general formula: 

Pa{y) =Pa{x) + Pa-l{x) X 2 n+l V a = 1, 2, . . . , 71 + 1. (72) 

Let’s use the “nabla” operators defined in Eq. (O: 


03 

03 

d o) 

( d 

d 

d 

9 \ 

V^Xl ’ dx2 

''dx n ' ) ’ 

|co 

1 

> 

’ 9x2 ’ 

dx n 

dxn~\- 1 y 


One has then, Va = 1,..., n + 1, for Eqs. (1721) and (1551) : 

V v Pa(y) = Vy [p a {x) + Pa- i(x) x^ +1 ] = V x p a (x) + x ^ +1 V x p a -i(x) + p a - i(x) 2x„+i e n+1 , (73) 

where e n +i is the canonical unit vector of K” +1 . 

We can now calculate a general element of the matrix P(y) of B n+ 1, using the basic invariant polynomials p a {y ), 
a = 1 ,...,n + 1, of -B n +i, and Eqs. EJ) and GD- For all a, b = 1,..., n + 1 we have: 

Pab(y) = v v p a (y) ■ VyPb(y) = 

(x) T X n _i_i ^ x Pa— 1 (x) 4“ Pa— 1 (x) 2 X^-j-i e^-l-l j • [V x pb(x) T X n _j_^ ^xPb—1 (x) T Pb—1 (x) 2 X n _|_i — 

= v x p a {x) ■ V x p b (x ) + x ^ +1 V x p a _i(x) ■ S7 X p b (x) + 0 + x £ +1 V x p a (x) ■ S7 x p b - l(x) + x ^ +1 V a .p 0 _i(x) • V x p 6 _i(x) + 


+ 0 + 0 + 0 + 4 x^ +1 Po -1 (x) Pb-i (x) = 

= Pab(x) +x^ +1 [Pa- 1 , fc(x) + Pa, 6-1 (x) + 4p a _i (x) p b -l (x)] + X^ +1 P a -l,6-l(x) ■ 

Of course, in the preceding expression, P c d(x) = 0, if one of the indices c, d (that in the above expression can be equal 
to a, b, a — 1, b — 1) is equal to 0 or to n + 1, because for B n the allowed indices run from 1 to n, and the gradients of 
Po{x) = 1, and p ra +i(x) = 0 are null vectors. By the induction hypothesis, Eq. (1551) is true for B n , so we can use it in 
the preceding equation. We obtain then, Va, b — 1,..., n + 1: 

min(a,6) —1 

Pab{y) = ^2 4(a + &- 1 - 2i) Pi(x) p a+ b-l-i(x) + 

i=max(0,a+6-n-l) 

{ min(a—1,6) —1 

Y. {a + b - 2 - 2i) Pi(x) p a+ b- 2 -i(x) + 

2 =max( 0 ,a+ 6 — n— 2) 


min(a,6—1) —1 

+ y (a + b - 2 - 2i) Pi(x) p a+b -2-i(x) + p a -i(x)p b -i(x) f + 

i=max(0,a+6— n— 2) 
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min(a—1,6—1) —1 

+ 4 x^ +1 Y (a + b - 3 - 2i) pi(x) p a +b- 3 -i(x) ■ 

z=max(0,a+6— n— 3) 

The term inside the curly brackets multiplying x‘^ L+1 can be simplified. One notes that both sums contain equal 
terms up to the value i = min(o — 1, b — 1) — 1 = min(a, b) — 2. Moreover the first sum contains one more term if 
min(a — 1,6—1) ^ min(a — 1, 6), that is if a > 6, that corresponds to i = b — 1 and the second sum contains one 
more term if min(a — 1,6—1)^ min (a, 6—1), that is if a < 6, that corresponds to i = a — 1. These two terms never 
appear simultaneously and one can write that the additional term appears for i = min(a, 6) — 1. One can then write 
the term inside the curly brackets multiplying x„ +1 as follows: 

{ min(a,6)—2 

Y 2 (a + b - 2 - 2i) Pi(x) p a+b - 2 -i{x) + 

i=max(0,a+6— n— 2) 


+ (a + 6 — 2 — 2(min(a, 6) - 1)) p min ( 0 ,b)- i(x) p a +b- 2 -(min(a,b)-i)( x ) + p a -i{x) p b -i(x)} 

One can see from this expression that if a = 6 there is no additional term originating from the two sums, as it should 
be. The terms in the second line of the above expression, if a < 6 can be written in the following way: 

(a + 6 — 2 — 2(a — 1)) p a -i(x) P a +b-2-(a-i)(x) + p a -i{x)Pb-i{x) = 


= (-a + 6) p a - iO) pb-i(x) + p a -i{x)p b -i{x) = (6- a + 1) p a -i(x) p b -i(x) 


(74) 


and if a > b can be written in the following way, obtained by interchanging a and 6: (a — 6 + 1) p a -i(x) pb-i(x). In 
all cases, one can then write the term inside the curly brackets multiplying as follows: 

{ min(a,6)— 2 

Y, 2 (a + b - 2 - 2i) Pi(x) p a+b - 2 -i(x) + {\b - cl\ + 1) p a -i(x) Pb-i{x) 

i=max(0,a+6— n— 2) 

In the sum multiplying 4x^ +1 , it is convenient to change the summation variable in i' = * + 1. One so finds: 

min(a—1,6—1) —1 min(a,6) —1 

Yj (a + 6 — 3 — 2i) Pi(x) p a+b -3-i(x) = Y (a + 6 - 1 - 2i') Pi>-i{x) p a+b - 2 -i'(x). 

i=max(0,a+6-n-3) i'=max(l,a+6— n— 2) 

All this gives: 

min(a,6) —1 

Pab(y) = Y 4 (a + 6 — 1 — 2i) pi(x) p a + b -i-i(x) + 

i=max(0,a+6— n— 1) 


+ 4 a;; 


n+1 


min(a,6)—2 


2 (a + 6 - 2 - 2 i) p t (x) p a+b - 2 -i{x ) + (|a - 6| + 1) p a - i(x) p b -i(x) 


i=max(0,a+6— n— 2) 


min(a,6) —1 

+ 4x* +1 Y (a + 6 -1 - 2i) Pi-i{x) p a+b - 2 -i(x) . (75) 

z=max(l,a+6— n— 2) 

This expression coincides exactly with that one one would find from Eq. (1551) in the case of B n+ 1. In fact, Eq. (1551) 
gives, Va<6 = l,...,n + 1: 

min(a,6) — 1 

Pab{p{y)) = Y 4(a + 6- 1 - 2 i) Pi(y) p a+6 _i_i(?/) = 

i=max(0,a+6— (n+1) — 1) 
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min(a,6) —1 


^ 4(a + 6-l-2i) (pi(x)+pi- 1 (x)x 2 n+1 ) (p a+b -i-i{x) + p a+h -2-i{x) x 2 n+1 ) 


i=max(0,a+6-n-2) 


min(a,6) — 1 

Y\ 4 (a + b - l - 2i) pi(x) p a +b-i-i(x) + 

i=max(0,a+6— n— 1) 


■4xl + i 


min(a,b) — 1 

T. (a + b - 1 - 2i) pi(x) p 0 +b _ 2 -i(a;) + 

2 =max( 0 ,a+&—n—2) 


min(a,6) — 1 

+ Y (a + b - 1 - 2i) p,;_i(a;) p a+b _i_i(:r) 

i=max(l,a+6-n-l) 


min(a,6) — 1 

+ 4x* +1 Y (a + 6 — 1 — 2*) Pi_i{x) p a +b-2-i(x) ■ 

i=max(l,a+6-n-2) 

In the first and third sums the lower limit of the sum has been changed from max(0, a+b—n— 2) to max(0, a+b—n— 1), 
because the factor p a+ b-i-i(x) vanish when i = a + b—n — 2, and in the two last sums the lower limit of the sum 
has been changed from max( 0 , a + b — n — 2 ) (or from max( 0 , a + b — n — 1 ) ) to max(l, a + b — n — 2 ) (or to 
max(l, a + b — n — 1 ) ), because the factor p,;_ \{x) vanish when * = 0 . 

The first and the last sums coincide exactly with the first and the last sums in Eq. m- The expression inside the 
curly brackets needs to be transformed somewhat to see that it is identical to the expression inside the curly brackets 
in Eq. (USD- The road to follow is indicated by the fact one wants to obtain the sum that appears inside the curly 
brackets in Eq. m- One has: 


min(a,6) — 1 min(a,b) — 1 

Y (a + b- 1 - 2i) pi(x) p a +b- 2 -i{x) + Y (a + b - 1 - 2i) Pi-i{x) p a+b -i-i(x) = 

i=max(0,a+6-n-2) i=max(l,a+6— n— 1) 


min(a,b) — 1 min(a,b)—2 

(a + b-l-2i) pi(x) p a+b - 2 -i{x)+ Y (a + b-l-2{i'+ 1)) pi*(x) p a+b -i-(i'+i){x) = 

z=max(0,a+6—n—2) z / =max(0,a+6—n—2) 


min(a,6) — 1 

Y (°+ & 

i=max(0,a+6-n-2) 


min(a,b)—2 

1 - 2 i) pi(x) p a +b-2-i{x) + Y (a + b - 3 - 2i) pi(x) p a+b -2-i(x) = 

2 =max( 0 ,a+&— n— 2) 


min(a,b) — 1 min(a,6) — 1 

Y (a + b - 2 - 2i) pi(x) p a +b-2-i(x) + Y Pi{x) p a +b-2-i{x) + 

2 =max( 0 ,a+ 6 — n— 2) i=max(0,a+fo— n— 2) 


min(a,6)—2 min(a,fe)— 2 

+ Y (a + b - 2 - 2i) pi(x) p a +b- 2 -i(x) - Y Pi{x) p a+b -2-i(x) = 

2 —max(0,a+6—n—2) i=max(0,a+6- n— 2) 


min(a,b)—2 

Y 2 (a + b — 2 — 2i) pi(x) p a +b-2-i{x ) + 

i=max(0,a+6— n— 2) 
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+ 


a + b — 2 — 2(min(a, b) — 1) 


+ 1 


Pmin(a,b) — 1 (x) Pa+b— 2— (min(a,b) — 1) (x) 


min(a,6)— 2 

E 2 (a + b - 2 - 2i) pi(x) p a +b- 2 -i(x ) + (|a - b| + 1) p a -i(x) Pb-i(x). 

i=max(0,a+6—n—2) 


The simplification in the last line is obtained by distinguishing the cases a > b, a = b and a < b. This expression is 
identical to the expression inside the curly brackets in Eq. C5D- We have so proved that the whole expression for the 
matrix element P a b{y) obtained from the application of Eq. (1351) for the group B n+ 1 coincides with Eq. urn that is 
obtained from the induction principle and the assumed validity of Eq. (1551) for the group B n . This proves that Eq. 
© is valid for B n for all values of n. o 

Proof of the generating formulas (13611 . (13711 and (138(1 . 

The groups of type B n have two different root systems (and as crystallographic groups have two different root lengths 
in the ratio y/2). Theorem [2] then asserts that there are two different active factors s(p) and l(p) of det(P(p)) that 
are the re-homogeneous polynomials in pi,... ,p n , that represent the invariant polynomials product of the squares of 
the linear forms l r (x) = r ■ x, corresponding to the n positive short roots and to the 2Q) positive long roots of B n , 
respectively. Theorem [2] also tells us how to construct the corresponding A-vectors. 

Let’s start with the short roots. As recalled in the beginning of Section |Vj a set of n positive short roots are e*, 
i = 1,..., n. The linear forms vanishing in the reflection hyperplanes of the short roots are the variables themselves: 
Xi, i = 1,... ,n, and the polynomial s(p(x)), obtained as specified in Theorem [2j coincides with the highest degree 
basic invariant polynomial p n (x): 


s(p{x)) = ]\ x i — p n (x). 

2=1 

Its A-vector can be obtained from the first of Eqs. m ■ For all a = 1,..., n, one has: 


aW( P (x)) = 2 y 

reK+,s 


VPa(x) 

l r (x) 


r 




V Pa(x) ■ e t 
Xi 



dpgjx ) 
dxi 


The basic invariant polynomial p a {x), written in Eq. (13dl) . has (") terms, and the number of those containing Xi 
(actually xf) is ("~ 1 ). One then finds: 


l_ dp a (x) 
Xi dxi 


2 Xi 
Xi 


E 

d*i —i)C(l,2,...,v..,n) 


2 2 
X il X *2 


.. x 2 

X la-1 


= 2 


E 

v(*i ,i2,...,ia-i)C(l,2,...,i,...,n) 


2 2 2 
X il X l2 ’ ' ' X i a -1 


where the sum is over all the combinations of a — 1 integers (i\, % 2 , ■ ■ ■, i a - 1 ) among the first n integers, excluding the 
integer i (the hat in i means exclusion, and writing the symbol ■ ■ ■ ,i a - 1 ) we automatically suppose i\ < %2 < 

... < i a - 1 ). There are ("ij) such combinations. One then has: 


^a ] (p(x)) = 4 E 

i=l 

The second member is a completely symmetric polynomial with terms that are products of a — 1 different variables 
squared, among xi,...,x n . It must then be a multiple of the basic invariant polynomial p a _i(x). To find its 
coefficient one can count the number of terms in the expression above and in the expression of p a - i(x). The ratio of 
these two numbers gives the coefficient of p a _i(x) one is searching in. The expression at the second member above has 
ci = <-) terms, while the explicit expression of p a - i(x), given in Eq. has C 2 = ( Q Zi) terms. The coefficient 

of p a -i(x) one is searching in is then the following one: 

ci n(n—1)! (a — 1)! (n — a + 1)! 

C 2 (a — 1)! (n — a)! n! 




-i)C(l,2,. 


2 2 2 
X H X i 2 -" X i a . i 


,,n) 


(n — a + 1), 
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and one finally obtains: 

Ki S) (p(x )) = 4(n- a+ l)p a -i(x ), 


that concludes the proof of Eq. (13(1 . 

There is a much simpler proof of Eq. (1361) if one makes use of the explicit form of the P-matrix given in Eq. (1351) . One 
starts by verifying that p n is an active polynomial. To do this one has to verify if p n satisfies the boundary equation 
m with a proper A-vector A ^ Pn \p). Let’s take a(p) = p n in Eq. m- The first member of Eq. m then reduces to 
the following one: 


~ dp n ^ ^ 

y Pabip) TT^ = p ab(p) <>bn = Pan(j> ) = 4 (n - Cl + 1) p a -l p n , d = 1, . . . , 71 , 


6= 1 


dp b 


6=1 


If one defines the A-vector A ^ Prl \p) = (A i Pn ^(p),..., A n n \p)) such that 

Ai Pn) (p) = 4(ra- a+ l)p a -i , a = l,...,n, 

one sees that p n is active because it satisfies the boundary equation (fTIT) . Theorem 4.1 of Ref. 
establishes then that p n is a factor of det(P(p)). 


item iv), 


Let’s consider now the long roots. As recalled in the beginning of Section El a set of 2Q) positive long roots are 
± 6j. 1 < * < j < n. The linear forms vanishing in the reflection hyperplanes of the long roots are then: Xj ± Xj, 
1 < * < j < n, and the polynomial l(jp(x)), obtained as specified by Theorem [2], is the following one: 

n n 

l(p(x)) = (Xi + Xj)(Xi - Xj) = 

i<j=l i<j=l 

Its A-vector can be obtained from the second of Eqs. m- For all a = 1,..., n, one has: 


A«(p(x))=2 £ 

r£TZ+, 1 


VPajx) 

l r (x) 


r 


2 Vp a (x) 


E 

i<j=l 


Ci + Cj 
Xi + Xj 


Ci Cj 
Xi — Xj 


2 Vp a (x) 


E 

i<j =1 


2 Xi Ci 2 Xj Cj 
(Xi + Xj)(Xi - Xj) 


n 

= 4 E 

i<j —1 



(„ dp a (x) dp a (x)\ 
V r * dx.i Xj dxj ) ' 


(76) 


The basic polynomial p a (x ), written in Eq. (l34l) . has (”) terms, and the number of those containing Xi (actually xf) 
is . One then finds: 


dp a (x) 
dxi 


= 2 xj 


du.*2,. 


E 

o —i)C(l,2,. 


2 2 
X h X i2 


where the sum is over all the combinations of a — 1 integers (*i, i%,..., i a -i) among the first n integers, excluding the 
integer i. There are ("~ x ) such combinations. It is convenient to distinguish the terms containing the factor x'j from 
those that do not. We can then write: 


dpajx) = 2 
dxi 



2 2 2 2 
X i X j X i\ X i2 


E 



The expression for Xj is similar, with i and j interchanged. When making the difference of the two expressions 

so found, all the terms of the first sums cancel, because of the symmetry in i and j, and one finds: 


dpg(x) _ dp a (x) 

dxi X ' 3 dxj 


2 ( x i - x 2 i) 




E 

)C(1,2,. 


2 2 
X h X i 2 


., 2 ,. 


.,n) 
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Inserting this result in Eq. urn one obtains: 

n 

Ai°(p(*)) = 8 

i<j =1 

The second member is a completely symmetric polynomial with terms that are products of a — 1 different variables 
squared, among aq,...,x„. It must then be a multiple of the basic invariant polynomial p a _i(x). To find its 
coefficient one can count the number of terms in the expression above and in the expression of p a -i(x). The ratio of 
these two numbers gives the coefficient of p a -i{x) one is searching in. The expression at the second member above 
has Ci = ( 2 ) ("~i) terms, while the expression of p Q _i(x), given in Eq. dnu, has C 2 = ( terms. The coefficient of 
Pa-i(x) one is searching in is then the following one: 

ci n(n—l) (n— 2)! (a — 1)! (n — a + 1)! (n — a+l)(n — a) 

C 2 2 (a — 1)! (n — a — 1)! n\ 2 

and one obtains: 

A { a\p( x )) = 8—— a+ ^ ^ — —Pa- i(x) =A(n-a + l)(n-a)p a - 1 (x). 

Eq. (1371) is so proved. 

From Eq. m one then has: 

y(d et (P))( p ) = xl s \ p ) + ( p ) =4(n- a + l)p a _! + 4 (n - a + l)(n - a)p a _i = 4 (n - a + l) 2 p a -i ■ 

This concludes the proof of Eq. (l38l) . o 


v(*i>*2|...,*o-i)C(1,2,.. 


2 2 
X h X i 2 




D. Proof of Theorem [9l for D n 


Proof of the generating formula 1)40(1 . 

Comparing Eq. (13911 . giving the basic invariant polynomials of D n , with Eq. (1341) . giving the basic invariant poly¬ 
nomials of B n , one sees that, for a < n the basic invariant polynomials of D n are identical with the basic invariant 
polynomials of B n , and the last basic invariant polynomial p n (x) of B n is just the square of the last basic invariant 
polynomial p n (x) of D n . It is convenient in this proof to write p^ix) and p^{x) to distinguish the basic invariant 
polynomial p n (x) of B n and D n , respectively. One has thus (p®{x)) 2 = Pn{x). No confusion arises for the other 
basic invariant polynomials, so p a (x ), when a < n, needs no specification. This implies that for a,b < n, the matrix 
elements of the matrix P(x ), obtained through Eq. (J6J), when calculated with the basic invariants of D n or of B n , is 
the same, and when one has to use Eq. m to express these invariant polynomials in terms of the basic invariant 
polynomials, the expressions obtained for D n differ from those of B n only for the use of (p^ 3 ) 2 in place of p^. Then 
the matrix elements P^(p) of the P-matrix of D n are the same of the matrix elements P^(p) of the P-matrix of B n , 
except for the substitution —> (p^ 3 ) 2 , that is, the variable p n of B n must be substituted with the square of the 
variable p n of D n . Using Eq. (l35l) and the substitution p n —>■ p 2 , the first line of Eq. (l40l) follows immediately. 

To prove the second and third lines of Eq. (l40l) it is possible to calculate the matrix elements of the matrix P D (p) of D n 
using the explicit form of the matrix elements of the matrix P B (p) of B n , written in Eq. (EH). As p B {x) = {jPn O)) 2 , 
one has X/p B (x) = V(p^(x)) 2 = 2p B (x) Vp^( x). Then, Va = 1 ,n— 1: 


Pan( X ) = Vp a (a:) ■ Vp% (x) 


1 

2 pS(x) 


\7p a (x) ■ Vpf (x) 


1 

2 Pn{x) 


Pan(p ( x )) 


1 

2 Pn(x) 


4 (n— 


a + 1) Pa-l(x) Pn {x) = 


Pn{x) 


2 (n - a + 1) Pa-i(x) (p n (x)) =2 (n-a+ 1)p Q _i(x)p„ (x) 


The second line of Eq. (HOI) is thus proved. To prove the third line, we have: 

1 1 


P nn(x) = Vp n ( x ) ' ^Pn = 


4 (Pn (x))' 


■ Vp„ (x) ■ Vp„ (x) = 


4 Pn(x) 


Pnn(P(x)) = 


4 Pn(x) 


4Pn-l(x) Pn( x ) =p n -l{x) 
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The third line of Eq. (l40l) is thus proved and (1401) is then completely proved. 


o 


Proof of the generating formula ()41|) . 

Theorem [T| tells us how to construct the A-vector A \F\p) corresponding to the determinant det(P(p)) = D(p) of the 
P-matrix P(p). 

As recalled in the beginning of Section ED a set of 2(™) positive roots of D n is for example formed by the vectors 
e.i ± e.j, 1 < i < j < n, and the linear forms vanishing in the corresponding reflection hyperplanes are Xi ± Xj, 
1 < i < j < n. Eq. m gives then, for all a = 1 ,,n: 


^a D) (p(x)) = 2 

ren+ 


V Pa(x ) 
l r { x) 


r 


2 S7p a {x) 


E 

i<j =1 


Xi + Xj 


€i Gj 
Xi — Xj 


2 S7p a {x) 


E 

i<j=l 


2 G^ 2 Xj Gj 

(Xi + Xj)(Xi - Xj) 


n 

= * E 

i<j—l 



(„ dp a (x) __ dp a (x)\ 
V Xi dXi Xj dXj J ' 


For a = 1,..., n — 1, the basic invariant polynomials of D n are the same as those of B n , so the previous calculation 
has already been done in the proof of Eq. (EH). Using the result there obtained for A a\p(x)), corresponding to the 
long roots of B n , we can write: 

A] [ D \p) = 4 (n - o+ l)(n - a)p a ~ 1 , Va = l,...,n-1. 

This is equal to Eq. m, but only for a = 1,..., n — 1. It remains to calculate the case a = n. One has p n (x) = 
xix 2 ,---,x n , that implies: Vp n (x)= • • •, Pn{x). Then: 


*L D) (p(x ))= 4 E 


_ „2 
i<j=1 » J 


9p„(a:) 9p„(a;) 

--- Xn 


dxi 


dxj 


= 4 


~2 _ ~2 
i<j=l » J 


— Pn(») - —Pn(x) ) = 0. 


Eq. dSTj) is so completely proved. 


o 


E. Proof of Theorem 1101 

Let’s consider the generating Eq. (13511 . relative to the group B n , in which, for the symmetry, we can suppose a < b. 
The variable p n , representing the basic invariant polynomial of the highest degree 2n, can be found in the element 
P a b{p ) only in p a+ b_i_i when i = a+b—n— 1. From the lower limit of the sum in Eq. (1351) we know that i = a+b—n— 1 
is possible only ifa + 6 — n — 1 > 0, that isifa + 6>n+l. This tells us that p n can be found not above the auxiliary 
diagonal of the P-matrix, only. (This is obvious also from the calculation of the weights of the matrix elements of the 
P-matrix: w(P a b ( P )) = d a +db~ 2 = 2a + 25 — 2 = 2(a + 5) — 2, and above the auxiliary diagonal one has a + 5 < n +1 
that gives w(P a b(p)) < 2 n = d n ). The coefficient of p n in P a b{p ), supposing a + 5>n + l, is obtained from Eq. (l35l) 
with i = a + b — n—1, and is the following one: 4 (a + b — 1 — 2(a + 5 — n — l))p 0+ t,_„_i = 4 (2 n — a — b+1) Pa+b-n- 1 - 
This result depends on a and 5 only through the sum a + 5, and this means that all the elements of the P-matrix that 
are in a same up going diagonal (that is parallel to the auxiliary diagonal), contain p n with the same coefficient, the 
same of the element in the last row, that one corresponding to 5 = n. Eq. (1551) gives P na (p) = 4 (n — a+ l)p a -ip n , so 
the coefficient of p n in all elements of the up going diagonal starting at P na is 4 (n — a + 1) p a - 1 , V a = 1,..., n. Using 
the definition of the Hankel matrices Ha" 1 (p) , given just after Eq. (1421) . we can say that the terms in the matrix P{p) 
containing p n form the Hankel matrix 4 p n H^\p). 

Consider now the matrix P(p) — 4p n H„ n \p). It does not contain p n and its greatest weight variable is p n -i- The 
variable p n _i, when not multiplied by p n , is contained in the element P a b(p) only inp a+ b_i_i when i = a+b—n. From 
the lower limit of the sum in Eq. (1351) we know that i = a + 5 — n is possible only ifO < a + b — n< min(a, 5) — 1 = a — 1, 
and we obtain the inequalities 5 < n — 1 and a + b > n. This tells us that p n ~ 1 , when not multiplied by p n , can be 
found only not above the diagonal of the P-matrix connecting P n -i.i (p) with Pi jn _i(p), and not in the last row (and 
column) of the P-matrix. The coefficient of p n -i in P a b{p ), supposing a + b > n and 5 < n — 1, is obtained from Eq. 
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(l35l) with i = a + b — n, and is the following one: 4 (a + b — 1 — 2(a + b — n))p a +b-n = 4 (2 n — a — b — l)p a+ &_„. 
This result depends on a and b only through the sum a + 6 , and this means that all the elements of the matrix 
P(p) — 4 p n Hr' 1 ' 1 (p) that are in a same up going diagonal, contain p n -i with the same coefficient, the same of the 
element in the second-last row, that one corresponding to b = n — 1 in Eq. (1551) . From Eq. (1551) one sees that 
Pn-i,a{p), Va = 1 ,..., n — 1 , contains the term 4 (n — a)p a -ip n -i, so the coefficient of p n -i in all elements of the up 
going diagonal starting at P„_i iQ (p) is 4 (n — a) p a -i- Using the definition of the Hankel matrix B^"\(p), we can say 
that the terms in the matrix P(j>) — 4p n H^\p) containing p n _ i, but not p ni form the Hankel matrix 4p n _i h2\(p)- 

We can continue this way finding where p n -2 is contained in the matrix P(p) — 4 p n H^(p) — 4p„_i H^ l } 1 (p) and 
again for p n - 3 , p n - 4 , and so on, and the formula (1551) easily follows. o 


F. Proof of Theorem 1111 

Theorem m states that Eq. (1551) is equivalent to Eq. (1551) . A comparison of Eqs. m and m is sufficient to 
understand that the P- matrix of D n in the upper left (n — 1) x (n — 1) block is equal to the P-matrix of B n except 
for the substitution of p n with pj ). Using Eq. (1451) . one finds the first two terms at the second member of Eq. (1551) . A 
comparison of Eq. ®5J), for b = n, that implies i = a — 1, and the second line of Eq. pop . is sufficient to understand 
that the elements in the last row and column of the P-matrix of D n , except the (n,n) element, are one half of the 
corresponding elements of the P-matrix of B n . Using Eq. (1431) for a = n, this gives the third term at the second 
member of Eq. (1551) . The third line of Eq. (150l) is represented by the fourth term at the second member of Eq. (1441) . 
Eq. (1551) and Theorem fill are so proved. o 


G. Proof of Theorem 1121 

The P-matrix corresponding to the basic invariant polynomials of A n defined in Eqs. (1281) (1311) can be constructed 
using Eq. (1551) . Let’s consider the three terms in the expression of P a b(p ) given by Eq. (1551) . By considering all 
possible values of a, b = 1 ,,n, they define three n x n matrices. The matrix corresponding to the first term in Eq. 
(155]) corresponds to the symmetric matrix (p) defined in Eq. (l45l) , and corresponds to the first term in the second 
member of Eq. (1461) . The matrix corresponding to the second term in Eq. (1551) is a Hankel matrix because it depends 

only on the sum a + b, and corresponds to the matrix 2 ]Co=i( a + 1 )Pa Ya^ in Eq. (l46l) . We have then just to prove 
that the matrix T(p) corresponding to the third term in Eq. (1321) . that one with the summation symbol, corresponds 

to the matrix — (n + 1) P a Ka n \p) in Eq. (l46l) . From Eq. (1551) . we have: 


min(a,b) — 1 

T a b(p) = — (n + 1) ^2 (a + b - 2i) pi-i pa+b-i-i, Va, b = 1 ,..., n , 

i=max(2,a+6- n— 1) 


in which one has to consider p o = 0. We can simplify the formula defining the matrix T{p). For the symmetry of the 
matrix T(p) we can suppose a < b. Moreover, we can define p a = 0, Va < 0 and Va > n, and can then write: 

a— 1 

Tabip) = -(n + 1 ) ^ (a + 6 - 2 i) pj_i p a+ 6 _i_j, Va, b = 1 ,..., n. 

i=2 

We can now change the lower limit of the sum in the following way: 

a— 1 

T a b{p ) = ~{n + 1) ^2 (« + b - 2i) pt- 1 p a +b-i-i , Va, b = 1,..., n , 

i=i 

because for i = 1 there is a factor p 0 that is equal to 0. If a = 1, T a b(p) = 0, V 6 = 1,..., n, so T(p) has its first row 
(and column) with elements equal to 0. (Actually, also the second row (and column) of T(p) is zero, but this does 
not matter here). We can define a matrix T'(p) that has for its first n—1 rows the last n — 1 rows of the matrix T(p) 
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(that are so shifted one position towards the top). We can write the formula for the matrix elements of T'(p) just by 
replacing the row label a of the formula for T(p) with a + 1. One then has: 

a 

T’ ab (p) = T a+ltb (p) = — (n + 1) ^2 (a + b+ 1 - 2i) Pi-i p a +b-i, Va= 1,... ,n-1, V6 = 1,... ,n. 

2=1 

Forgetting T{jp) we can allow the index a run up to n. Then: 

a 

T'atip) = -(n + 1) ^ (a + 6 + 1 - 2*) pi_ x p a+b -i , Va, b = 1,..., n , 

2=1 

and changing the summation index from i to i' = i — 1, one can write: 

a— 1 

T ab(p ) = -(n + 1) ^ (a + 6 - 1 - 2 i') p v p a+b -i-i >, Va, b = 1,..., n. 
i'= 0 

With the same definitions for the variables p a , that is po = 0 and p a = 0, Va < 0 and Va > n, the P-matrix of B n , 
in Eq. (1351) . can be written in the following way: 

a— 1 

Pabip) = ^2 4 (a + b - 1 - 2i) pi p a+b -i-i , Va, b = 1,..., n . 

2=0 

Comparing with the expression of T'{jp) we see that the only difference, in addition to po = 0 instead of po = 1, is the 
coefficient: — (n + 1) in the case of T'{p) and 4 in the case of the P-matrix of B n . We can then use Theorem [TUI and 
Eq. (TI51) and write: 


T\p) = -{n + 1 )J2 Pa H a n) (P ) - 

a=1 

but we have to consider po = 0 instead of po = 1. We can now obtain T(p) just by shifting one position to the bottom 
the first n — 1 rows of T'(p), and setting all the elements in the first row of T(p) equal to zero. For this it is then 
convenient to define the matrices Ki n \p) that are obtained from the matrices Ha U \p ) just by shifting one position 
towards the bottom their first n — 1 rows, and by substituting p^ = 0 in place of po = 1. We have then: 

n 

Tip) = + 1 )J2 Pa K a n) (P ) > 

a =1 

and this completes the proof of Theorem [12] o 


XI. APPENDIX 

In this Appendix I list alternative expressions of the generating formulas (l25l) . (l32l) . (l35l) and (l40l) in which the 
symmetry of the resulting matrices is hidden, but that can, however, be of some interest. 

Formula equivalent to Eq. G3, for the groups of type S n : 

a— 1 

Pabip) = in + 1 - b) Pa -1 Pb -1 + ^2{a-b-2i) p a -i-i Pb-i+i, Va, b = 1,..., n , 

i= 1 

in which one has to consider po = 1 and p a = 0, V a > n. 

Formula equivalent to Eq. ([32]), for the groups of type A n : 

a—2 

P ab {p ) =a(n + l-b) Pa-i Pb-i + 2 (a + b) p a+b -1 + (n + 1) ^ (a - b - 2i) p a -i-i Pb-i+i , Va, b = 1,..., n , 

2=1 


57 


in which one has to consider po =Q and p a = 0, V a > n. 

Formula equivalent to Eq. (j35]), for the groups of type B n : 

a 

Pab{p ) = 4 ^2 (b - a - 1 + 2i) p a -i p b -i+i , Va, b = 1 ,..., n , 

i= 1 

in which one has to consider po = 1 and p a = 0, V a > n. 

Formula equivalent to Eq. (l40l) . for the groups of type D n : 

Pab{p ) = = 4 J2i =1 {b-a- 1 + 2 i) Pa _ iPb _ 1+i | pn ^ p 2 

Pan{p) = Pna{p) = 2 (n - a + 1) p a -l pn 

P‘n,n (p) Pn — 1 - 

where one has to consider po = 1 and p a = 0, V a > n. 

To obtain the formulas written in this appendix from those in Eqs. (E5l) . (15^1) . (l35l) and SO]), one starts to consider 
a < b, for the symmetry of the P-matrices obtained by Eqs. S2D, S3 and SOD , define p a = 0 V a < 0 and 

V a > n, in such a way to drop all min and max in Eqs. dm S3- S3 and S3- and then change the summation 
index from i to i' = a — i. The formulas written above generate symmetric matrices (a formal proof is not short so I 
do not reproduce here), so it is not necessary to consider the limitation a < b, but it is anyway convenient because 
if one considers a < b one finds a smaller number of terms originating from the sums. If a > b many of these terms 
cancel each others. 


Va, b = 1,..., n — 1, 
Va = 1,..., n — 1, 
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